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Abstract:  A  dynamic  simulation  modeling  approach  to  describing  carbon- 
flow-based,  ecophysiological  processes  and  biomass  dynamics  of  fresh¬ 
water  submersed  aquatic  plant  species  has  been  developed.  The  models 
describe  major,  carbon-flow-based  ecophysiological  processes  and  bio¬ 
mass  dynamics  of  four  common  freshwater  species  and  how  these  are  in¬ 
fluenced  by  factors  such  as  light,  temperature,  current  velocity,  dissolved 
inorganic  carbon  availability,  oxygen  concentration,  and  human  influences 
such  as  management  measures  (changes  in  turbidity,  mechanical  harvest¬ 
ing,  grazing,  flooding).  The  model  species  are  Vallisneria  americana 
[model  VALLA],  Potamogeton pectinatus  [POTAM],  Hydrilla  verticillata 
[HYDRIL],  and  Myriophyllum  spicatum  [MILFO].  These  plant  species  are 
similar  in  growth  strategy  but  differ  significantly  in  morphology  and 
physiology.  The  same  modeling  approach  was  followed  for  all  species,  with 
species-characteristic  morphology  and  physiology  incorporated  in  four 
separate  models.  V.  americana  and  P.  pectinatus  are  considered  as  desir¬ 
able,  and  H.  verticillata  and  M.  spicatum  as  invasive  species  in  parts  of 
the  United  States.  In  the  models,  phenology  is  tied  to  day-degree  sum. 

This  enables  simulations  for  different  sites  and  climates  facilitating  model 
operation.  The  models  contain  unique  descriptions  of  (l)  species- 
characteristic  vertical  distribution  of  shoot  biomass  in  the  water  column; 
(2)  recalculation  procedures  of  vertical  distribution  with  daily  changes  in 
water  level;  (3)  ibidem  with  daily  removal  of  shoot  biomass  at  various  lev¬ 
els  within  the  water  column;  (4)  species-characteristic  epiphytic  light  in¬ 
terception;  (5)  species-characteristic  effects  of  current  velocity  on  photo¬ 
synthesis;  (6)  removal  of  periodic  shoot  and  tuber/root  crown  biomass; 
and  (7)  relationships  of  plant  process  parameters  with  site-specific  climate 
by  linkage  with  formatted  weather  files  and  calculation  of  latitude-specific 
effects  on  day  length.  Sensitivity  analysis  showed  that  the  models  are  very 
sensitive  to  changes  in  process  parameters  influencing  carbon  flow.  A  good 
fit  was  found  between  simulated  and  measured  biomass  in  the  field.  The 
models  can  be  used  to  simulate  the  daily  changes  in  plant  processes,  bio¬ 
mass,  and  persistence  under  various  site-specific  conditions  and  manage¬ 
ment  scenarios  over  a  1-  to  5-year  period  of  four  submersed  aquatic  vege¬ 
tation  species,  in  climates  varying  from  temperate  to  tropical. 
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Preface 

This  report  summarizes  the  expanded  development  of  aquatic  plant 
growth  models  that  can  be  used  as  stand-alone  versions  or  linked  with  a 
variety  of  hydrologic  and  sediment  modeling  systems.  This  development 
effort  was  performed  by  the  U.S.  Army  Corps  of  Engineers  (USACE)  Engi¬ 
neer  Research  and  Development  Center  (ERDC),  Environmental  Labora¬ 
tory  (EL),  Vicksburg,  MS.  Funding  for  this  expansion  effort  was  provided 
under  the  System-Wide  Water  Resources  Program  (SWWRP).  Dr.  Steven 
L.  Ashby  was  Program  Manager  of  SWWRP.  Previous  funding  for  the  de¬ 
velopment  of  the  versions  1.0  of  these  models  was  provided  by  the  USACE 
Aquatic  Plant  Control  Research  Program  (APCRP);  and  for  the  refinement 
of  VALLA  and  POTAM  for  use  in  the  Upper  Mississippi  River  System  was 
provided  by  the  USACE  Rock  Island  District. 

Special  thanks  go  to  Drs.  C.  F.  Cerco,  Water  Quality  and  Contaminant 
Modeling  Branch,  and  D.  M.  Soballe  (Environmental  Processes  Branch, 
USACE)  who  reviewed  and  provided  comments  on  an  earlier  version  of 
this  report. 

This  study  was  performed  by  Dr.  Elly  P.  H.  Best,  Principal  Investigator,  of 
the  Environmental  Risk  Assessment  Branch  (ERAB),  in  cooperation  with 
William  A.  Boyd,  who  assisted  with  programming  and  model  runs,  of  the 
Ecosystem  Processes  and  Engineering  Branch,  both  of  the  Ecological 
Processes  and  Engineering  Division  (EPED).  The  study  was  conducted 
under  the  general  supervision  of  Dr.  Robert  P.  Jones,  Chief,  ERAB; 

Richard  E.  Price,  Chief,  EPED;  Dr.  Mike  Passmore,  Acting  Deputy 
Director,  and  Dr.  Beth  C.  Fleming,  Director  of  EL. 

COL  Richard  B.  Jenkins  was  Commander  and  Executive  Director  of  ERDC. 
Dr.  James  R.  Houston  was  Director. 
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1  Introduction 

Submersed  aquatic  vegetation  (SAV)  may  play  important  roles  in  aquatic 
ecosystems.  Roles  attributed  to  desirable  species  are  stabilization  of  sedi¬ 
ment,  amelioration  of  transparency  and  regulation  of  nutrient  availability 
in  the  water  column,  and  serving  as  habitat  and  food  source  for  inverte¬ 
brates,  fish,  and  waterfowl.  Many  SAV  communities  in  freshwater  and  ma¬ 
rine  environments  have  experienced  dramatic  losses  during  the  past  three 
to  five  decades.  Most  commonly  declines  are  attributed  to  decreases  in  wa¬ 
ter  transparency  due  to  anthropogenic  influences,  but  often  they  have 
been  attributed  to  other  factors  such  as  high  water  levels,  extended  draw¬ 
downs,  changes  in  current  velocity,  and  epiphyte  shading,  or  to  combina¬ 
tions  of  factors.  However,  once  the  vegetation  is  lost  from  a  given  locale, 
increased  sediment  resuspension  and  current  velocity  may  place  signifi¬ 
cant  constraints  on  plant  recolonization  of  that  site.  Roles 
attributed  to  nuisance  or  invasive  SAV  species,  however,  are  excessive 
biomass  production  interfering  with  human  utilization  of  freshwater  re¬ 
sources  or  displacing  desirable  indigenous  communities. 

The  degree  to  which  SAV  influences  the  ecosystem  is  proportional  to  plant 
mass  and  depends  on  plant  species  and  physical  and  chemical  factors. 
Therefore,  predictions  of  the  environmental  impact  of  management  meas¬ 
ures  concerning  the  aquatic  system  in  which  SAV  grows  should  be  based 
on  accurate  estimates  of  plant  (l)  species,  mass,  and  pertinent  physiologi¬ 
cal  properties;  (2)  contribution  to  various  food  chains  (live  and  dead  bio¬ 
mass);  (3)  contribution  to  biogeochemical  cycling  (carbon,  nitrogen,  phos¬ 
phorus)  and  oxygen  regime;  and  (4)  contribution  to  habitat  availability. 
Simulation  models  that  include  descriptions  of  SAV  responses  in  terms  of 
biomass  dynamics  to  changes  in  physical  and  chemical  factors  in  various 
climates  can  be  useful  tools  for  water  resource  managers  because  they  can 
be  used  to  evaluate  key  environmental  conditions  in  which  SAV  would 
persist  or  produce  excessive  biomass,  with  ensuing  consequences  for  the 
systems  in  which  they  grow,  either  affected  or  not  affected  by  management 
scenarios  (Carr  et  al.  1997;  Best  et  al.  2001).  Although  the  number  of 
simulation  models  for  production  of  monotypic  SAV  is  increasing  (e.g., 
Titus  et  al.  1975;  Best  1981;  Collins  and  Wlosinski  1985;  Best  and  Jacobs 
1990;  Hootsmans  1994;  Scheffer  et  al.  1993;  Herb  and  Stefan  2003),  it  is 
still  relatively  low  compared  with  that  for  terrestrial  vegetation. 
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This  technical  report  describes  aquatic  vegetation  models  that  simulate 
major,  carbon  flow-based  ecophysiological  processes,  biomass  dynamics 
and  persistence  of  four  common  freshwater  species  and  how  these  are  in¬ 
fluenced  by  factors  such  as  light,  temperature,  current  velocity,  dissolved 
inorganic  carbon  (DIC)  availability,  oxygen  concentration,  and  human  in¬ 
fluences  such  as  management  measures  (changes  in  turbidity,  mechanical 
harvesting,  grazing,  flooding).  The  model  species  are  Vallisneria  ameri- 
cana  (American  wildcelery)-VALLA,  Potamogeton  pectinatus  (sago 
pondweed)-POTAM,  Hydrilla  verticillata  (dioecious  hy d r i  1 1  a) - H Y D RI L , 
and  Myriophyllum  spicatum  (Eurasian  watermilfoil)-MILFO.  These  plant 
species  are  similar  in  growth  strategy,  but  they  differ  significantly  in  mor¬ 
phology  and  physiology.  Therefore,  the  same  modeling  approach  was  fol¬ 
lowed  for  all  species,  with  species-characteristic  morphology  and  physiol¬ 
ogy  incorporated  in  four  separate  models.  V.  americana  and  P.  pectinatus 1 
are  considered  as  desirable,  and  H.  verticillata  and  M.  spicatum  as  inva¬ 
sive  species  in  parts  of  the  United  States. 

The  models  can  be  used  to  simulate  the  daily  changes  in  plant  processes 
and  biomass  under  various  site-specific  conditions  and  management  sce¬ 
narios  over  a  l-  to  5-year  period  of  four  SAV  species,  in  climates  varying 
from  temperate  to  tropical. 


1  P.  pectinatus  has  relatively  recently  moved  from  the  Potamogetonaceae  (Voss  1972)  into  the 
Stukeniaceae  (Crow  and  Hellquist  2000),  and  its  current  taxonomic  name  is  S tukenia  pectinata.  The 
authors  use  the  taxonomic  name  commonly  cited  up  to  2000  in  this  paper,  since  all  literature  refer¬ 
ences  pertain  to  the  formerly  used  name  of  P.  pectinatus  for  this  plant. 
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2  Modeling  Concepts 

Each  model  simulates  growth  of  a  typical  SAV  community.  In  the  models 
growth  is  considered  as  the  plant  dry  matter  accumulation  including  sub¬ 
terranean  tubers  in  tuber-forming  plants  or  rhizome/root  crown  systems 
in  M.  spicatum,  under  ample  supply  of  nitrogen  and  phosphorus,  in  a 
pest-,  disease-,  and  competitor- free  environment  under  the  prevailing 
weather  conditions.  At  least  one  plant  cohort  (a  cohort  is  a  plant  group  ex¬ 
hibiting  the  same  phenological  cycle)  waxes  and  wanes  per  season  in  dif¬ 
ferent  climatological  regions,  varying  from  temperate  to  tropical.  The  rate 
of  dry  matter  accumulation  is  a  function  of  irradiance,  temperature,  car¬ 
bon  dioxide  (C02)  availability,  and  plant  characteristics.  The  rate  of  C02 
assimilation  (photosynthesis)  of  the  SAV  community  depends  on  the  radi¬ 
ant  energy  absorbed  by  the  canopy,  which  is  a  function  of  incoming  radia¬ 
tion,  reflection  at  the  water  surface,  attenuation  by  the  water  column,  the 
epiphytes,  plant  shoots,  and  associated  leaf  area.  The  daily  rate  of  gross 
C02  assimilation  of  the  community  is  calculated  from  the  absorbed  radia¬ 
tion,  the  photosynthetic  characteristics  of  individual  shoot  tips  and  the 
pH-determined  C02  availability.  Calculations  are  executed  in  a  set  of  sub¬ 
routines  added  to  the  model.  Part  of  the  carbohydrates  produced  is  used  to 
maintain  the  existing  biomass.  The  remaining  carbohydrates  are  con¬ 
verted  into  structural  dry  matter  (plant  organs).  In  the  conversion  process, 
part  of  the  weight  is  lost  in  respiration.  The  dry  matter  produced  is  parti¬ 
tioned  among  the  various  plant  organs  using  partitioning  factors,  defined 
as  a  function  of  the  phenological  cycle  of  the  community.  The  dry  weights 
(DW)  of  the  plant  organs  are  obtained  by  integration  of  their  growth  rates 
over  time.  The  plant  winters  through  tubers  (rhizome/root  crown  for 
M.  spicatum )  in  the  sediment  without  or  with  biomass  present.  All  calcula¬ 
tions  are  performed  on  a  square  meter  basis.  Since  environmental  factors 
and  plant  growth  characteristics  vary  with  depth,  in  each  model  the  water 
column  and  associated  growth-related  processes  have  been  partitioned  in 
o.io-m  depth  classes. 

Seed  formation  has  not  been  included  in  the  models  because  its  role  in 
maintaining  existing  SAV  communities  in  a  temperate  climate  is  minimal. 
Dispersal  and  colonization  of  new  habitats  by  seeds  are  recognized,  impor¬ 
tant  characteristics  of  SAV.  The  latter  processes,  however,  are  better  de¬ 
scribed  using  other  modeling  approaches  as  discussed  by  Scheffer  (1991), 
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i.e.,  based  on  logistic  regression  or  on  descriptions  of  population  dynamics 
varying  in  time  and  space. 
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3  Description  of  Models 

Model  approach,  variables,  and  implementation 

Model  approach 

The  models  follow  a  mechanistic  approach  explaining  plant  growth  on  the 
basis  of  the  underlying  processes,  such  as  phenology,  C02  assimilation, 
and  respiration,  as  influenced  by  environmental  conditions.  They  are 
carbon-flow-based  mass  balance  models.  This  type  of  model  follows  the 
state-variable  approach,  in  that  it  is  based  on  the  assumption  that  the  state 
of  each  system  can  be  quantified  at  any  moment  and  that  changes  in  the 
state  can  be  described  by  mathematical  equations.  In  this  type  of  model, 
state,  rate,  and  driving  variables  are  distinguished.  State  variables  are 
quantities  such  as  biomass.  Driving  variables  characterize  the  effect  of  en¬ 
vironment  on  the  system  at  its  boundaries,  such  as  climate  and  environ¬ 
mental  conditions.  Each  state  variable  is  associated  with  rate  variables  that 
characterize  its  rate  of  change  at  a  certain  instant,  as  a  result  of  specific 
processes.  These  variables  represent  flows  of  material  between  state  vari¬ 
ables,  the  values  of  which  are  calculated  from  the  state  and  driving  vari¬ 
ables  according  to  knowledge  of  the  physical,  chemical,  and  biological 
processes  involved.  After  all  rate  variables  are  calculated,  these  values  are 
then  used  to  calculate  the  state  variables  according  to  the  following 
scheme:  state  variable  at  time  t  +  At  equals  state  variable  at  time  t  plus  the 
rate  at  time  t  multiplied  by  At.  This  procedure,  called  numerical  integra¬ 
tion,  gives  the  new  values  of  the  state  variables,  from  which  the  calculation 
of  rate  variables  is  repeated.  To  avoid  instabilities,  the  time  interval  At 
must  be  small  enough  so  that  the  rates  do  not  change  materially  within 
this  period.  This  is  generally  the  case  when  the  time  interval  of  integration 
is  smaller  than  one-tenth  of  the  time  coefficient  or  response  time.  This 
characteristic  time  of  a  system  is  equal  to  the  inverse  of  the  most  rapid 
relative  rate  of  change  of  one  of  its  state  variables.  The  smaller  the  time 
coefficient,  the  smaller  the  time  interval  of  integration  (Rabbinge  and 
De  Wit  1989). 

The  predictive  ability  of  mechanistic  models  does  not  always  equal  expec¬ 
tations.  It  should  be  realized,  however,  that  each  parameter  estimate  and 
process  formulation  has  its  own  uncertainty,  and  that  uncertainties  in  pa¬ 
rameter  estimates  may  accumulate  in  the  prediction  of  the  final  yield.  The 


ERDC/EL  TR-07-14 


6 


primary  aim  of  this  model  approach  is  to  increase  insight  in  the  system 
studied  by  quantitatively  integrating  the  current  knowledge  in  a  dynamic 
simulation  model.  By  studying  the  behavior  of  such  a  model,  better  insight 
in  the  real  system  is  gained. 

Variables 

A  list  of  variables  used  in  the  models  discussed  in  this  report  is  presented 
in  Table  l. 

Implementation 

Each  model  is  implemented  as  a  FORTRAN77  program.  For  numerical  in¬ 
tegration,  the  Runge-Kutta  technique  is  used,  which  allows  employing  a 
variable  time-step.  The  program  integrates  the  equations  once  per  day  in 
the  main  subroutine  (MODEL),  once  per  second  in  the  subroutines  calcu¬ 
lating  day  length  and  instantaneous  irradiance  (ASTRO)  and  instantane¬ 
ous  gross  assimilation  (ASSIM),  and  at  three  times  of  the  day  in  the  sub¬ 
routine  calculating  daily  total  gross  assimilation  (TOTASS;  Gaussian 
integration).  Instantaneous  gross  assimilation  is  calculated  per  second  and 
converted  to  hourly  rates  within  ASSIM. 

Model  organization  is  similar  to  that  used  for  agricultural  crops 
(SUCROSi;  Goudriaan  et  al.  1992).  Each  model  runs  within  a  FORTRAN 
Simulation  Environment  (FSE)  shell,  Version  2.1,  to  enable  easy  handling 
of  input  and  output  files  and  rapid  visualization  of  the  simulation  results 
(Van  Kraalingen  1995).  In  the  model  input  file,  default  values  on  plant 
characteristics  and  environmental  conditions  are  provided.  Users  may  en¬ 
ter  site-specific  data  into  this  input  file.  In  the  weather  file,  data  on  irradi¬ 
ance,  temperature,  and  latitude/longitude  are  provided.  Users  may  enter 
site-specific  data  into  this  weather  file.  Two  types  of  output  files  are  gener¬ 
ated,  i.e.,  a  binary  file  and  an  ASCII  text  file.  The  binary  file  is  used  for 
visualization  of  the  results  within  the  shell.  The  ASCII  file  can  be  imported 
in  various  software  packages  and  serve  as  the  basis  for  the  preparation  of 
publishable  tables  and  graphics. 

Each  model  can  be  executed  on  IBM  PC-ATs  and  compatibles  as  a 
stand-alone  version.  Because  of  its  language  and  simple  structure,  the 
model  will  generally  be  compatible  with  ecosystem  models  that  accept 
FORTRAN.  Executable  model  versions  are  available  via  the  internet  at 
http://el.erdc.usace.mil/products.cfm?Topic=model&Type=aquatic.  The  organization  of  the 
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models  and  their  subroutines  in  combination  with  the  FSE  shell  is  illus¬ 
trated  in  Figure  l. 


Table  1.  Variable  listing  of  the  models. 


Name 

Explanation 

Dimension 

AH(i) 

Absolute  height  of  vegetation  on  top  of  stratum  1,  measured 
from  the  plant  top 

m 

AMAX 

Actual  CO2  assimilation  rate  at  light  saturation  for  shoots 

g  CO2  g1  DW  h1 

AMTMP 

Daytime  temperature  effect  on  AMX  (relative) 

- 

AMTMPT 

Table  of  AMX  as  function  of  DVS 

AMX 

Potential  CO2  assimilation  rate  at  light  saturation  for  shoots 

g  CO2  g1  DW  h1 

ASRQ 

Assimilate  requirement  for  plant  dry  matter  production 

g  CH2O  g 1  DW 

ASRQSO 

Assimilate  requirement  for  storage  component  production 

g  CH2O  g 1  DW 

ATMTR 

Atmospheric  transmission  coefficient 

- 

COSLD 

Intermediate  variable  in  calculating  solar  height 

- 

CRIFAC 

Critical  shoot  weight  per  0.1-m  depth  layer 

g  DW  per  (0.1-m  depth 
layer)  plant*1 

CRIGWT 

Critical  weight  per  0.1-m  depth  layer 

g  DW  per  (0.1-m  depth 
layer)  nr2 

CRRIZ* 

Critical  rhizome/root  crown  system  weight 

g  DW  m*2 

CVT 

Conversion  factor  for  translocated  dry  matter  into  CH2O 

- 

DAVTMP 

Daily  average  temperature 

°C 

DAY 

Day  number  (January  1=1) 

d 

DAYEM 

First  Julian  day  number 

d 

DAYL 

Day  length 

h 

DDELAY 

Integer  value  of  DELAY 

- 

DDTMP 

Daily  average  daytime  temperature 

°C 

DEC 

Declination  of  the  sun 

radians 

DELAY 

Lag  period  chosen  to  relate  water  temperature  to  air  temp.,  in 
cases  where  water  temp,  has  not  been  measured 

d 

DEPTH 

Water  depth 

m 

DLV 

Death  rate  of  leaves 

g  DW  rrr2  d*1 

DMPC(i) 

Dry  matter  allocation  to  each  plant  layer  (relative) 

- 

DMPCT 

Table  to  read  DMPC(i)  as  function  of  depth  layer  (relative) 

- 

DPTT 

Table  to  read  water  depth  as  a  function  of  day  no. 

m,  d 

DRT 

Death  rate  of  roots 

g  DW  rrr2  d*1 
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DSINB 

Integral  of  SINB  over  the  day 

s.d'1 

DSINBE 

Daily  total  of  effective  solar  height 

s.d'1 

DSO 

Daily  extra-terrestrial  radiation 

J  m  2  d'1 

DST 

Death  rate  of  stems 

g  DW  rrr2  d'1 

DTEFF 

Daily  effective  temperature 

°C 

DTGA 

Daily  total  gross  CO2  assimilation  of  the  vegetation 

g  CO2  rrr2  d1 

DTR 

Measured  daily  total  global  radiation 

J  nr2  d-1 

DVR 

Development  rate  as  function  of  temperature  sum 

d1 

DVRRT 

Table  of  post-anthesis  development  rate  as  function  of  tem¬ 
perature  sum 

Q. 

0 

O 

DVRVT 

Table  of  pre-anthesis  development  rate  as  function  of  tem¬ 
perature  sum 

CL 

O 

O 

DVS 

Development  phase  of  the  plant 

- 

EE 

Initial  light  use  efficiency  for  shoots 

g  CO2  J1 

EPISHD 

Table  to  read  fraction  of  irradiation  shaded  by  epiphytes  as 
function  of  DVS 

EPHSWT 

Relative  epiphyte  shading  switch  (O=off,  l=on) 

- 

EPHY 

Table  to  read  relative  epiphyte  shading  as  function  of  DVS 

FGROS 

Instantaneous  CO2  assimilation  rate  of  the  vegetation 

g  CO2  nr2  h-1 

FGL 

Instantaneous  CO2  assimilation  rate  per  vegetation  layer 

g  CO2  rrr2  h  1 

FL 

Leaf  dry  matter  allocation  to  each  layer  of  shoot  (relative) 

- 

FLT 

Table  to  read  FL  as  function  of  DVS 

FLV 

Fraction  of  total  dry  matter  increase  allocated  to  leaves 

- 

FLVT 

Table  to  read  FLV  as  function  of  DVS 

- 

FRDIF 

Diffuse  radiation  as  a  fraction  of  total  solar  radiation 

- 

FRT 

Fraction  of  total  dry  matter  increase  allocated  to  roots 

- 

FRTT 

Table  to  read  FRT  as  function  of  DVS 

FST 

Fraction  of  total  dry  matter  increase  allocated  to  stems 

- 

FSTT 

Table  to  read  FST  as  function  of  DVS 

GLV 

Dry  matter  growth  rate  of  leaves 

g  DW  nr2  d-1 

GPHOT 

Daily  total  gross  assimilation  rate  of  the  vegetation 

g  CH2O  nr2  d1 

GRT 

Dry  matter  growth  rate  of  roots 

g  DW  m-2  d-1 

GST 

Dry  matter  growth  rate  of  stems 

g  DW  m-2  d-1 

GTW 

Dry  matter  growth  rate  of  the  vegetation  (plant  excluding  tu¬ 
bers  or  rhizomes) 

g  DW  nr2  d1 

HAR 

harvesting  switch  (O=off,  l=on) 

- 
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HARDAY 

Harvesting  day  number 

d 

HARDEP 

Harvesting  depth  (measured  from  water  surface) 

m 

HIG(i) 

Height  on  top  of  stratum  1  (measured  from  water  surface) 

m 

HOUR 

Selected  hour  during  the  day 

h 

1 

Counter  in  DO  LOOP 

- 

IABS(i) 

Total  irradiance  absorbed  per  depth  layer  i 

J  nr 2  S'1 

IABSL(i) 

Total  irradiance  absorbed  by  shoots  per  depth  layer  i 

J  nr 2  S'1 

IDAY 

Integer  equivalent  of  variable  DAY 

d 

INTUB 

Tuber  size 

g  DW  tuber1 

IREMOB 

Initial  value  remobilization 

g  CH2O  m  2 

IRS 

Total  irradiance  just  under  the  water  surface 

J  m2  S'1 

IRZ(i) 

Total  photosynthetically  active  part  of  irradiance  on  top  of 
depth  layer  i 

J  m 2  s1 

IWGRIZ* 

Initial  rhizome  weight 

g  DW  nr2 

IWLVD 

Initial  dry  matter  of  dead  leaves 

g  DW  nr2 

IWLVG 

Initial  dry  matter  of  green  (live)  leaves 

g  DW  nr2 

IWRTD 

Initial  dry  matter  of  dead  roots 

g  DW  nr2 

IWRTG 

Initial  dry  matter  of  green  (live)  roots 

g  DW  nr2 

IWSTD 

Initial  dry  matter  of  dead  stems 

g  DW  m-2 

IWSTG 

Initial  dry  matter  of  green  (live)  stems 

g  DW  ITT2 

K 

Plant  species  specific  light  extinction  coefficient 

m2  g 1 DW 

KCOUNT 

Counter  used  to  calculate  number  of  consecutive  days  in 
which  sprouts  have  a  negative  net  photosynthesis 

- 

KT 

Table  to  read  K  as  function  of  DVS 

m2  g 1  DW,  - 

L 

Water  type  specific  light  extinction  coefficient 

m1 

LAT 

Latitude  of  the  site 

degrees 

LT 

Table  to  read  L  as  function  of  day  number 

d,  m-1 

MAINRT* 

Maintenance  respiration  rate  rhizome/root  crown  system 

g  CH2O  m  2  d1 

MAINT 

Maintenance  respiration  rate  vegetation 

g  CH2O  nr2  d1 

MAI  NTS 

Maintenance  respiration  rate  vegetation  at  reference  tempera¬ 
ture 

g  CH2O  nr2  d1 

NDTUB 

Dormant  tuber  density  (identical  to  NT) 

dormant  tubers  m  2 

NGLV 

Net  growth  rate  of  leaves 

g  DW  m-2  d-1 

NGRT 

Net  growth  rate  of  roots 

g  DW  m-2  d1 

NGST 

Net  growth  rate  of  stems 

g  DW  m-2  d1 
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NGTUB 

Sprouting  tuber  number 

spr.  tubers  rrr2 

NINTUB 

Tuber  number  concurrently  initiated  per  plant 

cone.  in.  tubers  plant1 

NNTUB 

New  tuber  number 

new  tubers  rrr2 

NPL 

Plant  density 

plants  rrr2 

NTM 

Tuber  density  measured  (field  site) 

tubers  rrr2 

NTMT 

Table  to  read  NTM  as  function  of  day  number 

tubers  rrr2,  d 

NTUBD 

Dead  tuber  number 

dead  tubers  rrr2 

NTUBPD 

Dead  tuber  number  previous  day 

dead  p.d.tubers  rrr2 

NUL 

Zero  (0) 

- 

PAR 

Instantaneous  flux  of  photosynthetically  active  radiation 

J  rrr2  S'1 

PARDIF 

Instantaneous  flux  of  diffuse  PAR 

J  rrr2  S'1 

PARDIR 

Instantaneous  flux  of  direct  PAR 

J  rrr2  S'1 

PI 

Ratio  of  circumference  to  diameter  of  circle 

- 

RAD 

Factor  to  convert  degrees  to  radians 

radians  degree1 

RC 

Reflection  coefficient  of  irradiance  at  water  surface  (relative) 

- 

RCSHST 

Relation  coefficient  tuber  weight-stem  length 

m  g 1  DW 

RDR 

Relative  death  rate  of  leaves  (on  DW  basis) 

d1 

RDRIZ* 

Relative  rhizome  death  rate 

d1 

RDRT 

Table  to  read  RDR  as  function  of  DAVTMP 

d1,  °C 

RDS 

Relative  death  rate  of  stems  and  roots  (on  DW  basis) 

d1 

RDST 

Table  to  read  RDS  as  function  of  DAVTMP 

d1  °C 

RDTU 

Relative  tuber  death  rate  (on  number  basis) 

d1 

RE  DAM 

Reduction  factor  to  relate  AMX  to  pH  and  oxygen  levels  of  the 
water  (relative) 

- 

REDAM1 

Table  to  read  factor  relating  AMX  to  water  current  velocity 
(relative) 

-,  cm  S'1 

REDF(i) 

Relative  reduction  factor  for  AMX  to  account  for  senescence 
plant  parts  over  vertical  axis  of  vegetation  (relative) 

- 

REMOB 

Remobilization  rate  of  carbohydrates 

gCFDO  nr2  d1 

ROC 

Relative  conversion  rate  of  tuber  into  plant  material 

gCFhOg1  DWd-1 

RTR 

Maximum  relative  tuber  growth  rate  at  20  °C 

g  DW  tuber1  d'1 

RTRL 

Relative  tuber  growth  rate  at  ambient  temperature 

g  DW  tuber1  d'1 

SC 

Solar  constant  corrected  for  varying  distance  sun-earth 

J  rrr2  s-1 

SC(i) 

Shoot  dry  matter  in  depth  layer  i 

g  DW  rrr2 

SHTBIO 

Shoot  biomass;  one  term  for  sum  WLV  +  WST 

g  DW  rrr2 
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SINB 

Sine  of  solar  elevation 

- 

SINLD 

Intermediate  variable  in  calculating  solar  declination 

- 

STEMLE 

Stem  length 

m 

SURFAC 

Expression  of  warning  that  plant  canopy  is  not  at  water  surface 
and  tuber  class  has  died 

- 

SSURPR 

Integer  value  of  SURPER 

- 

SURPER 

Survival  period  for  sprouts  without  net  photosynthesis 

d 

TBASE 

Base  temperature  for  juvenile  plant  growth 

°C 

TEFF 

Factor  accounting  for  temperature  effect  on  maintenance  res¬ 
piration,  remobilization,  maximum  relative  tuber  growth  and 
death  rates 

TEFFT 

Table  to  read  TEFF  as  function  of  temperature  (Q10  of  2,  up  to 
45  °C) 

-,  °C 

TGRIZ* 

Total  live  dry  weight  rhizome/root  crown  system  of  previous 
day 

g  DW  m-2 

TGW 

Total  live  plant  dry  weight  (excluding  tubers) 

g  DW  m-2 

TGWM 

Total  live  plant  dry  weight  measured  (field  site) 

g  DW  m-2 

TGWMT 

Table  to  read  TGWM  as  function  of  day  number 

g  DW  m-2,  d 

TL 

Thickness  per  depth  layer 

0.1  m 

TMAX 

Daily  maximum  temperature 

°C 

TMIN 

Daily  minimum  temperature 

°C 

TMPSUM 

Temperature  sum  after  1  January 

°C 

TRAFAC* 

Translocation 

- 

TRANS 

Translocation  rate  of  carbohydrates 

gCFDO  rn-2  d1 

TREMOB 

Total  remobilization 

g  DW  m-2 

TW 

Total  live  +  dead  plant  dry  weight  (excluding  tubers) 

g  DW  m-2 

TWCTUB 

Total  critical  dry  weight  of  new  tubers 

g  DW  m-2 

TWGRIZ* 

Total  live  dry  weight  rhizome/root  crown  system  of  current  day 

g  DW  m-2 

TWGTUB 

Total  dry  weight  of  sprouting  tubers 

g  DW  m-2 

TWLD(1,2,3)* 

Total  dry  weight  of  dead  leaves 

g  DW  m-2 

TWLG(1,2,3)* 

Total  dry  weight  of  live  leaves 

g  DW  m-2 

TWLVD 

Total  dry  weight  of  dead  leaves 

g  DW  m-2 

TWLVG 

Total  dry  weight  of  live  leaves 

g  DW  m-2 

TWNTUB 

Total  dry  weight  of  new  tubers 

g  DW  m-2 

TWRD(1,2,3)* 

Total  dry  weight  of  dead  roots 

g  DW  m-2 

TWRG(1,2,3)* 

Total  dry  weight  of  live  roots 

g  DW  m-2 
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TWRIZD* 

Total  weight  dead  rhizome/root  crown  system 

g  DW  nr2 

TWRTD 

Total  dry  weight  of  dead  roots 

g  DW  nr2 

TWRTG 

Total  dry  weight  of  live  roots 

g  DW  nr2 

TWSD(1,2,3)* 

Total  dry  weight  of  dead  stems 

g  DW  nr2 

TWSG(1,2,3)* 

Total  dry  weight  of  live  stems 

g  DW  nr2 

TWSTD 

Total  dry  weight  of  dead  stems 

g  DW  m-2 

TWSTG 

Total  dry  weight  of  live  stems 

g  DW  m-2 

TWTUB 

Total  dry  weight  of  tubers 

g  DW  m-2 

VELSWT 

Switch  to  relate  AMX  to  current  velocity  (O=off,  l=on) 

- 

WLV 

Dry  weight  of  leaves  (live  +  dead) 

g  DW  m-2 

WRT 

Dry  weight  of  roots  (live  +  dead) 

g  DW  m-2 

WST 

Dry  weight  of  stems  (live  +  dead) 

g  DW  m-2 

WTMP 

Daily  water  temperature 

°C 

WTMPT 

Table  to  read  WTMP  as  function  of  day  number 

°C,  d 

WVEL 

Water  type  specific  current  velocity 

cm  s1 

YRNUM 

Year  number  simulation  (1-5) 

- 

^Variables  unique  for  MILFO. 


Figure  1.  Relational  diagram  illustrating  the  organization  of  each  model  and  its  subroutines 
in  combination  with  the  FSE  shell. 
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Model  processes  and  parameter  values 

Model  processes 

Morphology,  phenological  cycle,  and  development 

Morphology  and  phenological  cycle.  All  four  plant  species  are  rooted  into 
the  sediment.  Summer  habit,  reproduction  strategy,  flowering  period,  and 
frequency  may  vary  strongly  with  plant  species.  Flowering  usually  coin¬ 
cides  with  peak  biomass  and  is  followed  immediately  by  sloughing. 

Phenology  and  development.  The  phenology  of  a  plant  community,  for 
which  development  phase  can  be  used  as  a  measure,  quantifies  physiologi¬ 
cal  age  and  is  related  to  the  plants’  morphological  appearance. 

Development  phase  can  not  be  expressed  simply  as  chronological  age  be¬ 
cause  several  environmental  factors  such  as  temperature  and  stress  (e.g., 
nutrients,  grazing)  can  speed  up  or  reduce  the  rate  of  phenological  devel¬ 
opment.  Contrary  to  what  is  suggested  by  intuition,  the  rate  of  plant 
growth  per  se  has  no  effect  on  phenological  development,  as  long  as  the 
growth  rate  is  not  very  low  (Penning  de  Vries  et  al.  1989b,  and  citations 
therein).  The  concept  of  development  phase  is  used  to  characterize  the 
whole  plant  community;  it  is  not  appropriate  for  individual  organs. 

The  response  of  developmental  rate  to  temperature  in  the  current  models 
is  in  accordance  with  the  degree-day  hypothesis  (Thornley  and  Johnson 
1990a).  The  idea  is  as  follows.  The  mean  temperature  T)  for  each  day  i  is 

measured,  and  a  sum  h  is  formed  according  to  the  following  equation: 

h  =  i(%-Tc)  (1) 

1=1 


which  includes  only  those  terms  where  T{  is  above  some  threshold  value 

Tc.  When  h  reaches  a  particular  value,  this  signifies  that  a  phase  in  devel¬ 
opment  is  complete,  and  this  is  generally  associated  with  a  biological  event 
that  occurs  over  a  short  period  of  time  and  is  readily  observed.  The  day- 
degree  sum  h  essentially  integrates  some  underlying  temperature- 
dependent  processes.  For  each  plant  species  there  are  various  phases  in 
the  development  of  the  plant,  and  the  temperature  sum  is  found  to  have  a 
certain  value  for  the  successful  completion  of  each  phase.  The  temperature 
threshold  Tc  may  be  different  for  each  of  these  phases.  The  approach  is 
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based  on  the  notion  of  a  developmental  rate,  whose  response  to  tempera¬ 
ture  is  approximately  linear  over  a  restricted  temperature  range.  Compari¬ 
son  with  actual  temperature  responses  found  in  agricultural  crops  suggests 
that  this  is  reasonable,  and  the  method  works  well  in  practice.  It  is  implic¬ 
itly  assumed  that  the  organ  possesses  a  developmental  clock  that  is  pro¬ 
ceeding  at  the  rate  kd .  In  general,  it  is  to  be  expected  that  the  development 
rate  kd  may  depend  on  a  number  of  quantities.  This  can  be  represented  by 
the  following  equation: 

kd=f{V,P,E )  (2) 

in  which/ represents  some  function  of  the  state  variables  V ,  parameters  P, 
and  environmental  quantities  E.  The  temperature-sum  rule  works  because 
the  most  important  environmental  variable  is  temperature,  and  the  re¬ 
sponse  to  temperature  is  approximately  linear. 

The  phenological  cycle  in  each  model  is  described  based  on  published  field 
data  collected  in  a  temperate  climate,  and  development  rates  have  been 
generated  by  scaling  the  development  phase  to  the  pertinent  climatological 
day-degree  sum.  In  cases  where  climatological  data  were  lacking  for  the 
year  when  the  field  data  were  collected,  climatological  data  for  the  same  or 
a  nearby  site  of  a  different  year  were  used,  after  verifying  climatological 
conditions  did  not  deviate  from  the  usual  at  that  site.  Three  plant  species, 
i.e.,  V.  americana,  P.  pectinatus,  and  H.  verticillata,  are  modeled  as  if  they 
propagate  vegetatively  by  tubers  with  and  without  plant  shoots  present, 
because  this  is  the  predominant  wintering  mechanism  supported  by  the 
calibration  data  that  were  collected  in  a  temperate  climate.  Tubers  are 
relatively  small,  dormant  organs  that  develop  on  most  stolons  under  spe¬ 
cial  day  length  and  temperature  conditions  in  autumn,  and  that  grow  into 
the  sediment.  Tubers  are  composed  of  a  small  amount  of  dividing  tissue 
surrounded  by  several  fleshy  leaves.  The  parent  plants  senesce  and  disin¬ 
tegrate  at  the  end  of  the  growth  season,  and  only  the  tubers  hibernate 
within  the  sediment  until  their  emergence  the  following  spring,  which 
completes  the  annual  growth  cycle.  Tubers  are  depleted  and  disintegrate 
during  the  summer  in  which  they  were  formed  in  V.  americana  and 
P.  pectinatus,  but  may  survive  several  years  in  H.  verticillata  (Van  and 
Steward  1990).  The  fourth  plant  species,  M.  spicatum,  is  modeled  as  if  it 
propagates  vegetatively  by  rhizome/root  crowns  with  and  without  plant 
shoots  present. 


ERDC/EL  TR-07-14 


15 


Development  phase  (DVS)  is  a  state  variable  in  the  models.  The  develop¬ 
ment  phase  is  dimensionless,  and  its  value  gradually  increases  within  a 
growing  season.  The  development  rate  has  the  dimension  d1.  The  multiple 
of  rate  and  time  period  yields  an  increment  in  phase.  In  the  models,  the 
temperature  affecting  plant  development  can  be  chosen  as  equal  to  the 
daily  average  air  temperature  at  the  height  of  the  growing  point  of  the 
shoots,  with  a  user-defined  lag-period  to  correct  for  deviations  in  tempera¬ 
ture  of  the  water  body  in  which  the  SAV  grows  compared  with  air  tempera¬ 
tures  (7  days  is  nominal).  It  is  more  accurate  to  use  water  temperatures  for 
this  purpose,  but  since  water  temperatures  are  not  always  available  for  the 
site  for  which  the  user  wants  to  run  the  model,  the  models  can  be  run  us¬ 
ing  either  value. 

The  rate  of  phenological  development  can  be  affected  by  temperature  dif¬ 
ferently  in  the  vegetative  phase  and  in  the  reproductive  phase.  These  dif¬ 
ferences  indicate  that  the  physiological  process  of  development  may  not  be 
the  same  before  and  after  flowering. 

Plant  density  and  maximum  biomass 

In  each  model,  calculations  are  performed  on  a  square  meter  basis  assum¬ 
ing  a  typical  plant  density.  Typical  plant  density  indicates  in  this  case,  that 
it  is  possible  that  at  some  point  in  time  under  natural  conditions  plant 
densities  may  be  different,  but  that  an  established,  monotypic  SAV  opti¬ 
mizes  at  this  typical  number  of  plants  per  square  meter.  Lower  densities 
may  occur  in  the  establishment  phase,  when  some  plants  may  not  yet  have 
neighbors  and  become  relatively  large,  whereas  higher  plant  densities  may 
occur  early  in  the  season  when  more  sprouts  emerge  but  the  plantlets  are 
subsequently  thinned  to  the  typical  number  of  plants  per  square  meter  by 
self-shading  of  the  vegetation.  In  tuber-forming  species,  this  implies  also 
that  the  number  of  sprouting  tubers  in  the  tuber  bank  equals  the  typical 
number  of  plants  per  square  meter,  while  the  remaining  tubers  continue  to 
senesce.  However,  at  tuber  bank  densities  lower  than  the  typical  number 
per  square  meter,  the  number  of  sprouting  tubers  is  recalculated  and  set 
equal  to  the  actual  tuber  bank  density.  If  wintering  plants  are  present, 
plant  biomass  is  redistributed  over  the  typical  number  of  plants  per  square 
meter. 

Seasonal  plant  biomass  maxima  may  vary  considerably  over  time  and 
space.  In  the  models  maximum  biomass  is  calculated  as  explained  in 
“Light,  photosynthesis,  maintenance,  growth,  and  assimilate  partitioning” 
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up  to  a  level  that  equals  the  highest  published  biomass  value.  Higher  bio¬ 
mass  values  are  considered  as  unrealistic  and  attributed  to  selective  sam¬ 
pling  or  erroneous  extrapolations  of  biomass  measured  on  small  surface 
areas  to  per  square  meter  area. 

Wintering,  sprouting  in  spring,  and  growth  of  sprouts  to  water  surface 

In  tuber-forming  SAV,  tubers  in  the  tuber  bank  within  the  sediment  are 
the  main  storage  organs  for  carbohydrates  during  wintering  in  a  temperate 
climate.  Tuber  bank  weight  is  calculated  by  multiplication  of  tuber  bank 
density  and  typical  tuber  weight.  Tuber  bank  density  may  decrease  by  tu¬ 
ber  death  with  a  typical  tuber  death  rate,  and  through  grazing  by  water- 
fowl  and  other  animals.  Tuber  bank  weight  may  decrease  by  maintenance, 
tuber  death  and  by  the  sprouting  of  tubers,  which  transform  into  plants. 
Both  tuber  bank  weight  and  density  may  increase  by  the  formation  of  new 
tubers.  Tubers  lie  dormant  if  not  disturbed  (Van  and  Steward  1990),  and  it 
is,  therefore,  to  be  expected  that  maintenance  processes  proceed  at  a  very 
low  level  of  activity.  Tuber  death  rate,  maintenance,  and  growth  increase 
with  increasing  temperature  through  a  relative,  effective  temperature 
function  ( TEFF ),  being  a  Q10  of  2  with  a  20  °C  reference  temperature. 

A  certain  minimum  weight  is  required  for  a  tuber  to  sprout.  The  number  of 
tubers  that  sprout  is  set  equal  to  the  typical  plant  density.  Sprouting  of  tu¬ 
bers  is  a  function  of  development  phase  (e.g.,  for  V.  americana  a  DVS  of 
0.291).  Plantlets  remobilize  the  tuber  carbohydrate  reserves  by  conversion 
into  plant  material,  and  elongate.  The  plant  material  is  allocated  to  leaves, 
stems  and  roots  following  a  fixed  biomass  partitioning  pattern  (see  “Light, 
photosynthesis,  maintenance,  growth,  and  assimilate  partitioning”). 

These  plantlets  may  survive  at  the  plant  height  they  can  maximally  reach, 
provided  the  balance  of  the  carbohydrate  influx  from  tuber  carbohydrate 
reserves,  efflux  from  tuber  growth  respiration,  and  influx  from  plantlet 
photosynthesis  is  positive.  Sprouts/plantlets  die  if  they  have  a  negative  net 
assimilation  rate  over  a  user-defined  number  of  days,  and  the  program 
stops  with  a  warning  “KCOUNT.”  After  the  death  of  one  tuber  class,  the 
next  tuber  class  may  sprout,  provided  more  tubers  are  available  in  the  tu¬ 
ber  bank  and  DVS  is  lower  than  required  for  flowering  (DVS  1).  Up  to 
three  tuber  classes  may  sprout  annually.  The  program  can  resume  running 
for  the  same  year  after  the  user’s  pressing  “ENTER”  provided  the  proper 
conditions  are  met. 
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A  relational  diagram  illustrating  the  wintering  and  sprouting  tubers  of 
tuber-forming  plant  species  is  shown  in  Figure  2. 


TEMPERATURE  LIGHT 


Figure  2.  Relational  diagram  illustrating  wintering  and  sprouting  of  tubers  in  tuber-forming  SAV. 
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TWGTUB 

=  NPL  x  INTUB 

(3) 

NDTUB 

=  NDTUB  -  (NTUBD  -  NTUBPD) 

(3.a) 

NTUBD 

=  RDTU  x  NDTUB  x  TEFF 

(3.a.l) 

IF  (DVS.  GE.  0.291.  AND.  DVS  .LT.  1.)  THEN 
IF  (TWGTUB  .LE.  (0.01  x  NPL  x  INTUB))  NGTUB  =  0.0 


NGTUB 

=  NPL 

(3.1) 

TWGTUB 

=  INTGRL  (TWGTUB,  -REMOB,  DELT) 

(3.2) 

REMOB 

=  TWGTUB  x  ROC  x  TEFF 

(3.3) 

In  M.  spicatum,  rhizomes/root  crowns  in  the  sediment  are  the  main  stor¬ 
age  organs  for  carbohydrates  during  wintering  in  a  temperate  climate.  An 
initial  biomass  equal  to  the  lowest  biomass  published  (g  DW  nr2,  meas¬ 
ured)  is  used  to  enable  the  model  plants  to  sprout.  Rhizome  weight  may 
decrease  by  rhizome  death  with  a  typical  rhizome  death  rate,  and  through 
grazing  by  waterfowl  and  other  animals.  Rhizome  weight  may  also  de¬ 
crease  by  maintenance  and  by  the  sprouting  of  the  growing  points,  which 
transform  into  plants.  Rhizome  weight  may  increase  by  downward  trans¬ 
location  from  shoots. 

The  number  of  plants  that  sprout  is  set  equal  to  the  typical  plant  density. 
Sprouting  of  growing  points  on  the  rhizome  is  a  function  of  development 
phase.  It  occurs  between  DVS  0.375  and  the  flowering  period  for  cohort  1, 
between  DVS  1.001  and  the  flowering  period  for  cohort  2.  In  a  tropical 
climate,  a  third  cohort  may  sprout  between  DVS  2.001  and  the  flowering 
period  for  cohort  3. 

Plantlets  remobilize  the  rhizome  carbohydrate  reserves  by  conversion  into 
plant  material  and  elongate.  The  plant  material  is  allocated  to  leaves, 
stems,  and  roots  following  a  fixed  biomass  partitioning  pattern  (see  “Light, 
photosynthesis,  maintenance,  growth,  and  assimilate  partitioning”). 
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These  plantlets  may  survive  at  the  plant  height  they  can  maximally  reach, 
provided  the  balance  of  the  carbohydrate  influx  from  tuber  carbohydrate 
reserves  depletion,  efflux  from  tuber  growth  respiration,  and  influx  from 
plantlet  photosynthesis  is  positive.  Sprouting  plantlets  die  if  they  have  a 
negative  net  assimilation  rate  over  a  user-defined  number  of  days,  and  the 
program  stops  with  a  warning  “Vegetation  is  dying.”  After  initiation  of 
flowering  of  plant  cohort  l,  the  next  plant  cohort  2  may  sprout,  provided 
more  rhizome  mass  is  available  and  DVS  is  lower  than  required  for  flower¬ 
ing  of  cohort  2  (DVS  1.630).  Up  to  three  plant  cohorts  may  sprout  annu¬ 
ally.  The  program  can  resume  running  for  the  same  year  after  the  user’s 
pressing  “ENTER”  provided  the  proper  conditions  are  met. 


TWGRIZ  =  INTGRL  [  TWGRIZ, 


^ -TRANS  1  +  REMOBi  +  MAINRT  ' 

v  1-242  y 


+  ( RDRIZ  x  TGRIZ )  , 


BELT' 


(3.4) 


IF  (DVS.  GE.  0.376.  AND.  DVS.  LT.  1.0)  THEN 
IF  (TWGRIZ.  GT.  CRRIZ )  THEN 


REMOBi 

=  ROC  x  TWGRIZ 

(3.5) 

TWRIZD 

=  INTGRL  ( TWRIZD ,  RDRIZ,  DELT) 

(3.6) 

where  1.242  equals  assimilate  requirement  for  rhizome  dry  matter  produc¬ 
tion  (g  CH20  g  DW1). 

A  relational  diagram  illustrating  the  wintering  and  sprouting  rhizomes/ 
root  crowns  of  M.  spicatum  is  shown  in  Figure  3. 

Light,  photosynthesis,  maintenance,  growth,  and  assimilate  partitioning 

Light.  Light  availability  is  an  important  factor  controlling  the  distribution 
and  abundance  of  SAV.  In  aquatic  systems  a  small  part  of  the  irradiance 
can  be  reflected  by  the  water  surface,  and  further  attenuation  occurs  by 
water  and  its  suspended  solids,  by  epiphytes,  and  by  SAV  itself.  Measured 
daily  total  irradiance  (wavelength  300-3,000  nm)  is  used  as  input  in  the 
model.  Only  half  of  the  irradiance  reaching  the  water  surface  is  considered 
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Figure  3.  Relational  diagram  illustrating  wintering  and  sprouting  of  the  rhizome/root  crown 
of  M.  spicatum. 


to  be  photosynthetically  active  and  is,  therefore,  used  as  a  base  for  the  cal¬ 
culation  of  C02  assimilation.  Six  percent  of  the  irradiance  is  reflected  by 
the  water  surface  (Golterman  1975).  The  subsurface  irradiance  is  attenu¬ 
ated  by  dissolved  substances  and  particles  within  the  water  column  result¬ 
ing  in  a  site-  and  season-specific  water  extinction  coefficient.  The  remain¬ 
ing  radiation  may  be  further  reduced  by  epiphyte  shading.  Epiphyte 
shading  can  be  activated  and  deactivated  in  the  model  input  file.  The  verti¬ 
cal  profiles  of  the  radiation  within  the  SAV  layers  are  characterized  also. 
The  absorbed  irradiance  for  each  horizontal  SAV  layer  is  derived  from 
these  profiles.  The  SAV  light  extinction  coefficient  K  is  plant  species- 
characteristic  and  assumed  to  be  constant  throughout  the  year.  The  in- 
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coming  irradiance  is  attenuated  by  the  shoots,  part  of  which  is  absorbed  by 
the  photosynthetic  plant  organs,  i.e.  the  leaves. 


IRZi+1=IRZixe 


(- TLxL-KxSC; ) 


(4) 


ilRZ.-IRZ^ixSC^K  . 

IABS .  =  - — - - - - x  1.0  -  EPISHD) 

1  ( KxSC+TLxL )  v  ’ 


(4.1) 


IABSL-  =  IABS-  x  FL 


(4.2) 


Photosynthesis.  Instantaneous  rates  of  gross  assimilation  are  calculated 
from  the  absorbed  light  energy  and  the  photosynthesis  light  response  of 
individual  shoots,  here  used  synonymously  to  leaves. 

The  photosynthesis-light  response  of  leaves  is  described  by  the  following 
exponential  function 


FGL  =  SCixAMAXx 


i-exp 


' -EE  x  IABS  {x  3600^ 
v  AM  AX  x  SCt  , 


(4.3) 


In  the  photosynthesis-light  response  equation,  the  value  of  potential  pho¬ 
tosynthetic  activity  at  light  saturation  (AMX)  is  species-characteristic  and 
the  initial  light-use  efficiency  (EE)  typical  for  C3  plants.  AMX  is  affected  by 
temperature  via  a  fitted,  relative  function,  AMTMPT,  accounting  for  the 
measured  effect  of  daytime  temperature,  and  enabling  the  calculation  of 
the  actual  photosynthesis  rate  (AMAX).  AMAX  is  also  affected  by  current 
velocity  via  a  fitted,  relative,  function,  REDAMi,  accounting  for  the  meas¬ 
ured  effect  of  current  velocity  on  AMX.  Senescence  and  daily  changes  in 
pH  and  oxygen  concentrations  may  affect  AMX.  Substituting  the  appro¬ 
priate  value  for  the  absorbed  photosynthetically  active  radiation  yields  the 
assimilation  rate  for  each  specific  shoot  layer.  The  instantaneous  rate  of 
gross  assimilation  over  the  height  of  the  vegetation  is  calculated  by  relat¬ 
ing  the  assimilation  rate  per  layer  to  the  species-characteristic  biomass 
distribution  and  by  subsequent  integration  of  all  vegetation  layers.  The 
daily  gross  assimilation  rate  is  calculated  by  using  the  Gaussian  integra¬ 
tion  method.  This  method  specifies  the  discrete  points  at  which  the  value 
of  the  function  to  be  integrated  has  to  be  calculated,  and  the  weighting  fac¬ 
tors  that  must  be  applied  to  these  values  to  attain  minimum  deviation 
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from  the  analytical  solution.  A  three-point  method  performs  well  for  calcu¬ 
lating  daily  total  assimilation  (Goudriaan  1986;  Spitters  1986). 

Maintenance,  growth,  and  assimilate  partitioning .  A  portion  of  the  car¬ 
bohydrates  formed  is  respired  to  provide  energy  for  maintaining  the  exist¬ 
ing  plant  components.  The  maintenance  costs  increase  with  metabolic  ac¬ 
tivity,  probably  because  of  higher  enzyme  turnover  and  higher  transport 
costs  (Penning  de  Vries  1975).  The  maintenance  cost  can  be  estimated 
from  the  chemical  composition  of  the  plant.  Typical  maintenance  coeffi¬ 
cients  for  various  plant  organs  have  been  derived,  based  on  numerous 
chemical  determinations  in  agricultural  crops.  They  typically  range  from 
0.010  to  0.016  g  CH20  g  AFDW1  d1  (Penning  de  Vries  and  Van  Laar 
1982b).  These  coefficients  have  been  used  in  the  current  models.  Mainte¬ 
nance  respiration  increases  with  increasing  temperature  through  a  relative 
effective  temperature  function  (TEFF),  as  described  under  “Wintering, 
sprouting  in  spring,  and  growth  of  sprouts  to  water  surface.” 

Equations  describing  maintenance  costs  are: 

MAINTS  =  0.016  x  TWLG  +  0.010  x  TWSG  +  0.015  x  TWRG  (4.4) 

MAINT  =  MAINTS  x  TEFF  (4.4a) 

The  remainder  of  the  carbohydrates  is  available  for  conversion  into  struc¬ 
tural  plant  material,  resulting  in  growth.  In  this  conversion  process  of  the 
glucose  molecule,  C02  and  water  are  released.  The  assimilate  requirement 
to  produce  one  unit  weight  of  any  particular  plant  organ  can  be  calculated 
from  its  chemical  composition  and  the  assimilate  requirements  of  the 
various  chemical  components.  Typical  values  are  1.46  g  CH20  g  DW'1  for 
leaves,  1.51  for  stems,  and  1.44  for  roots  (Penning  de  Vries  and  Van  Laar 
1982b;  Penning  de  Vries  et  al.  1989a),  confirmed  by  Griffin  (1994).  At 
higher  temperatures  the  conversion  processes  are  accelerated,  but  the 
pathways  are  identical.  The  recently  determined  construction  costs  for 
several  submersed  plant  species,  using  a  different  method  (Williams  et  al. 
1987),  are  generally  lower  and  range  from  0.99  to  1.11  (Spencer  et  al. 

1997).  However,  some  of  the  latter  plants  are  relatively  poor  in  nitrogen 
and  transport  costs  have  not  been  included,  both  factors  which  may  have 
contributed  to  the  lower  cost  found. 
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In  the  models  the  construction  costs  typical  for  agricultural  plants  have 
been  used  because  construction  costs  calculated  for  M.  spicatum  leaves 
with  an  average  chemical  composition  were  similar  to  those  in  agricultural 
plants,  i.e.,  1.465  CH20  g  DW1,  and  those  of  stems  and  roots  were  pre¬ 
sumed  to  be  similar  also. 


The  following  equation  describes  growth: 

[(REMOB  x  CVT )  +  GPHOT  -  TRANS  -  MAINT ] 
~~  ASRQ 


(4.5) 


Assimilate  partitioning  is  the  process  by  which  assimilates  available  for 
growth  are  partitioned  over  leaves,  stems,  roots,  and/or  storage  organs.  It 
depends  on  physiological  age.  Assimilate  partitioning  patterns  in  most 
SAV  species  are  unknown.  However,  partitioning  of  the  biomass  resulting 
from  this  process  in  full-grown  plants  has  been  published,  and  these  val¬ 
ues  have  been  used  for  calibration  in  the  models.  The  species- 
characteristic  distribution  of  shoot  biomass  is  allocated  over  the  vertical 
axis  via  a  dry  matter  partitioning  coefficient  function  (DMPC),  which  al¬ 
lows  changes  in  shoot  biomass  distribution  with  changes  in  water  depth 
and  differs  significantly  between  species. 


The  following  equation  describes  biomass  partitioning  over  plant  organs: 


GLV  =  FLV  x  GTW 

(4.6) 

GRT  =  FRT  x  GTW 

(4.7) 

GST  =  FST  x  GTW 

(4.8) 

A  relational  diagram  illustrating  photosynthesis,  respiration,  and  biomass 
formation  is  shown  in  Figure  4. 
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TEMPERATURE  LIGHT  C02  AVAILABILITY  PHOTOSYNTHETIC 

EFFICIENCY 


Figure  4.  Relational  diagram  illustrating  photosynthesis,  respiration,  and  biomass  formation 
in  tuber-forming  SAV. 


Flowering,  translocation,  senescence,  and  formation  of  wintering  organs 

The  flowering  event  affects  subsequent  metabolic  activity  of  the  vegeta¬ 
tion.  The  timing  of  flowering  is,  therefore,  extremely  important  for  physio¬ 
logical  activity  and  biomass  formation,  whereas  the  actual  investment  of 
dry  matter  in  flowers  and  seeds  is  minor.  In  the  models,  a  considerable 
part  of  net  production  is  translocated  downwards  after  flowering  in  the 
form  of  carbohydrates  considered  as  equivalent  to  starch  (Gijzen  1985).  In 
tuber-forming  species  the  translocation  products  are  largely  transported  to 
the  tubers,  in  M.  spicatum  to  the  rhizome/root  crowns.  The  remainder  of 
net  production  is  allocated  according  to  the  species-characteristic  biomass 
partitioning  pattern.  Senescence  sets  in  resulting  in  loss  of  particulate 
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plant  material  via  a  species-characteristic  relative  death  rate,  which  in¬ 
creases  with  an  increase  in  temperature  through  the  previously  mentioned 
relative,  effective  temperature  function  (TEFF). 

In  tuber-forming  plants,  tubers  are  induced  after  flowering,  i.e.,  DVS>i.o, 
and  formed  under  short  day  conditions  and  within  a  certain  temperature 
range.  The  combination  of  tuber  size  and  concurrently  initiated  tuber 
number  is  species-characteristic  (see  “Wintering,  sprouting  in  spring,  and 
growth  of  sprouts  to  water  surface”).  Once  initiated,  a  tuber  class  grows 
with  a  typical  actual  relative  tuber  growth  rate  (RTRL)  from  translocated 
material  until  a  preselected  critical  weight  of  the  new  tuber  class 
(TWCTUB)  is  reached.  Temperature  influences  the  maximum  relative 
growth  rate  (RTR)  of  tubers  through  the  TEFF  function.  Once  finished,  a 
tuber  class  is  added  to  the  dormant  tuber  bank,  and  the  plant  starts  form¬ 
ing  a  new  tuber  class.  Tuber  initiation  continues  as  long  as  environmental 
conditions  permit,  and  tubers  are  formed  as  long  as  plants  provide  assimi¬ 
lates  to  fill  them. 

The  following  equations  describe  induction  and  formation  of  new  tubers 
(illustrated  for  V.  americana;  DVS  and  temperature  range  differ  between 


species). 

IF  (REMOB  .  EQ.  o.o)  THEN 

IF  (DVS.  GT.  l.o.  AND.  DAYL.  LT.  14.7)  THEN 

IF  (DDTMP  .GT.  5.0  .AND.  DDTMP  .LT.  25.0)  THEN 
IF  (TGW .GT  .0.1 )  THEN 

TRANS  =  AMAXi  [o.,  (RTRL x l./CVT) x(GPHOT  -  MAINT )]  (5) 

NNTUB  =  NPL  x  NINTUB  (5.1) 

TWNTUB  =  INTGRL  ( TWNTUB,  TRANS,  DEBT )  (5.2) 

IF  (TWNTUB  .GE.  TWCTUB)  THEN 
NDTUB  =  NDTUB  +  ( NPL  x  NINTUB )  (5.3) 
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A  relational  diagram  illustrating  translocation  and  senescence  following 
anthesis  in  tuber-forming  plants  is  shown  in  Figure  5. 


Figure  5.  Relational  diagram  illustrating  translocation  and  senescence  following  anthesis 
in  tuber-forming  SAV. 


In  M.  spicatum,  rhizomes  increase  in  mass  through  translocation  after 
flowering,  i.e.,  DVS>i.o.  Once  started,  the  rhizomes  grow  with  a  typical 
relative  translocation  rate  (TRAP AC)  from  translocated  material  that 
decreases  exponentially  to  o  with  concomitantly  decreasing  biomass  of 
the  translocating  plant  cohort  and  increasing  biomass  of  the  successive 
growing  plant  cohort.  New  plant  cohorts  can  only  sprout  until  a  DVS  of 
2.001  is  reached. 
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Translocation  is  described  by  the  following  equation: 

TRANSi  =  CVT  x  GPHOT 

(TWLG2  +  TWSG2  +  TWRG2)  ]  (5_/ 

X  (TWLGi  +  TWSGi  +  TWRGi  +  TWLG2  +  TWSG2  +  TWRG2)  x  TRAFAC 

Parameter  values 

The  relatively  simple  simulation  model  type  described  in  this  paper  in¬ 
cludes  parameter  values  that  can  be  defined  with  varying  certainty.  Most 
parameters  have  been  calculated  or  estimated  from  published  literature 
(Table  2).  The  choice  of  parameter  values  is  summarized  below.  The  pa¬ 
rameter  values  are  presented  in  four  species-  and  model-characteristic 
tables:  Table  2a  pertains  to  V.  americana,  model  VALLA,  Table  2b  to 
P.  pectinatus,  model  POTAM,  Table  2c  to  H.  verticillata,  model  HYDRIL, 
and  Table  2d  to  M.  spicatum,  model  MILFO. 


Table  2a.  Parameter  values  used  in  VALLA. 


Parameter 

Name 

Value 

Reference 

Morphology,  phenological  cycle,  and  development 

Fraction  of  total  dry  matter  increase  allocated  to 
leaves 

FLV(T) 

0.718 

1,  2 

Fraction  of  total  dry  matter  increase  allocated  to 
stems 

FST(T) 

0.159 

1,  2 

Fraction  of  total  dry  matter  increase  allocated  to 
roots 

FRT(T) 

0.123 

1,  2 

Development  rate  as  function  of  temperature* 

DVRV 

0.015  d1 

Calibrated 

DVR  prior  to  flowering  DVRV,  subsequently  DVRR 

DVRR 

0.040  d1 

Calibrated 

Plant  density  and  maximum  plant  biomass 

Plant  density 

NPL 

30  nr2 

1 

Maximum  plant  biomass 

496  g  DW  m-2 

2 

Wintering,  sprouting,  and  growth  of  sprouts  to  water  surface 

Dormant  tuber  density 

NDTUB 

233  m 2 

3 

Tuber  size 

INTUB 

0.090  g  DW  tuber 1 

3,  4 

Rel.  tuber  death  rate  (on  number  basis) 

RDTU 

0.018  d1 

1 

Rel.  conversion  rate  of  tuber  into  plant  material 

ROC 

0.0576  g  CH2O. 
g1  DW  d-1 

5 

Relation  coefficient  tuber  weight-stem  length 

RCSHST 

12  m  g1  DW 

5,  6 

Critical  shoot  weight  per  0.1-m  depth  layer 

CRIFAC 

0.0091g  DW  (0.1-m 
depth  layer)1  plant1 

3,4 

Survival  period  for  sprouts  without  net  photosynthe¬ 
sis 

SURPER 

23  d 

7,  8 
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Light,  photosynthesis,  maintenance,  growth,  and  assimilate  partitioning 

Plant  species  specific  light  extinction  coefficient 

K(T) 

0.0235m2  g1  DW 

9 

Fraction  of  irradiation  shaded  by  epiphytes* 

EPISHD 

0-0.43 

10 

Potential  CO2  assimilation  rate  at  light  saturation  for 
shoots 

AMX 

0.0165  g  CO2. 
g-1  DW  h1 

9 

Initial  light  use  efficiency  for  shoots 

EE 

0.000011  gC02  J1 

11 

Rel.  reduction  factor  for  AMX  to  account  for  senes¬ 
cence  plant  parts 

REDF(T) 

1.0 

User  def. 

Reduction  factor  to  relate  AMX  to  water  pH 

RE  DAM 

1 

Daytime  temperature  effect  on  AMX  as  function  of 
DVS* 

AMTMP(T) 

0-1 

12 

Rel.  factor  to  relate  AMX  to  water  current  velocity* 

REDAM1 

0-1 

10 

Dry  matter  allocation  to  each  plant  layer* 

DMPC 

0-1 

9 

Flowering,  translocation,  senescence,  and  formation  of  wintering  organs 

Max.  relative  tuber  growth  rate  (part  of  net  photosyn¬ 
thetic  rate) 

RTR 

0.247  d1 

4, 12,  13 

Conversion  factor  for  translocated  dry  matter  into 

CH2O 

CVT 

1.05 

11 

Tuber  number  concurrently  initiated  per  plant 

NINTUB 

5.5  plant1 

13 

Total  critical  dry  weight  of  new  tubers 

TWCTUB 

14.85  g  DW  m-2 

1,  3,  13 

Rel.  death  rate  of  leaves  (on  DW  basis;  Q10=2) 

RDR(T) 

0.021  d1 

1 

Rel.  death  rate  of  stems  and  roots  (on  DW  basis; 
Q10=2) 

RDS(T) 

0.021  d1 

1 

Site  information 

Water  depth  (field  site) 

DPT(T) 

1.4  m 

User  def. 

Daily  water  temperature  (field  site) 

WTMP(T) 

-,  °C 

User  def. 

Water  type  specific  light  extinction  coefficient  (field 
site) 

L(T) 

0.43-0.80  m  1 

1 

Water  type  specific  current  velocity  (field  site) 

WVEL 

0-100  cm  S'1 

User  def. 

Total  live  dry  weight  measured  (field  site) 

TGWM(T) 

-,  g  DM  rrr2 

User  def. 

Tuber  density  measured  (field  site) 

NTM(T) 

233  m-2 

1 

Harvesting 

Harvesting 

HAR 

0  or  1 

User  def. 

Harvesting  day  number 

HARDAY 

1-365 

User  def. 

Harvesting  depth  (measured  from  water  surface;  1-5 
m) 

HARDEP 

0.1  m<DEPTH 

User  def. 

*  Calibration  function. 

1.  Titus  and  Stephens  1983.  2.  Haller  1974.  3.  Korschgen  and  Green  1988.  4.  Korschgen  et  al.  1997.  5.  Bowes  et  al.  1979. 
6.  Best  and  Boyd  1996.  7.  Titus  and  Adams  1979b.  8.  Best  and  Boyd  2003c.  9.  Titus  and  Adams  1979a.  10.  Bestetal. 
2005.  11.  Penning  de  Vries  and  Van  Laar  1982  a,  b.  12.  Don nermeyer, 1982.  13.  Donnermeyer  and  Smart  1985. 
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Table  2b.  Parameter  values  used  in  POTAM. 


Parameter 

Name 

Value 

Reference 

Morphology,  phenological  cycle,  and  development 

Fraction  of  total  dry  matter  increase  allocated  to 
leaves 

FLV(T) 

0.731 

1,  2 

Fraction  of  total  dry  matter  increase  allocated  to 
stems 

FST(T) 

0.183 

1,  2 

Fraction  of  total  dry  matter  increase  allocated  to 
roots 

FRT(T) 

0.086 

1 

Development  rate  as  function  of  temperature* 

DVRV 

0.015  d1 

Calibrated 

DVR  prior  to  flowering  DVRV,  subsequently  DVRR 

DVRR 

0.040  d1 

Calibrated 

Plant  density  and  maximum  plant  biomass 

Plant  density 

NPL 

30  m2 

1,  3 

Maximum  plant  biomass 

1,952  g  DW  m  2 

4 

Wintering,  sprouting,  and  growth  of  sprouts  to  water  surface 

Dormant  tuber  density 

NDTUB 

240  m 2 

1 

Tuber  size 

INTUB 

0.083  g  DW  tuber1 

1 

Rel.  tuber  death  rate  (on  number  basis) 

RDTU 

0.026  d1 

5 

Rel.  conversion  rate  of  tuber  into  plant  material 

ROC 

0.0576  g  CH2O. 
g 1  DW  d-1 

6 

Relation  coefficient  tuber  weight-stem  length 

RCSHST 

12  m  g1 DW 

6,  7,  8 

Critical  shoot  weight  per  0.1-m  depth  layer 

CRIFAC 

0.0076  g  DW  (0.1-m 
depth  layer)-1  plant1 

7,  8 

Survival  period  for  sprouts  without  net  photo¬ 
synthesis 

SURPER 

27  d 

1 

Light,  photosynthesis,  maintenance,  growth  and  assimilate  partitioning 

Plant  species  specific  light  extinction  coefficient 

K(T) 

0.095  m2  g1  DW 

1 

Fraction  of  irradiance  shaded  by  epiphytes* 

EPISHD 

0-1.0 

9 

Potential  CO2  assimilation  rate  at  light  saturation  for 
shoots 

AMX 

0.019  g  C02. 
g1  DW  h  1 

10 

Initial  light  use  efficiency  for  shoots 

EE 

0.000011  gC02  J'1 

11 

Rel.  reduction  factor  for  AMX  to  account  for  senes¬ 
cence  plant  parts 

REDF(T) 

1.0 

User  def. 

Reduction  factor  to  relate  AMX  to  water  pH 

RE  DAM 

1 

1 

Daytime  temperature  effect  on  AMX  as  function  of 
DVS* 

AMTMP(T) 

0-1 

1 

Rel.  factor  to  relate  AMX  to  water  current  velocity* 

REDAM1 

0-1 

9 

Dry  matter  allocation  to  each  plant  layer* 

DMPC 

0-1 

1 

Flowering,  translocation,  senescence,  and  formation  of  wintering  organs 

Max.  relative  tuber  growth  rate  (part  of  net  photosyn¬ 
thetic  rate) 

RTR 

0.19 

1,  12 

Conversion  factor  for  translocated  dry  matter  into 

CH2O 

CVT 

1.05 

11 

Tuber  number  concurrently  initiated  per  plant 

NINTUB 

8  plant1 

1,  8 
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Total  critical  dry  weight  of  new  tubers 

TWCTUB 

19.92  g  DW  m  2 

1,  3 

Relative  death  rate  of  leaves  (on  DW  basis;  Q10  = 2 ) 

RDR(T) 

0.047  d1 

1 

Relative  death  rate  of  stems  and  roots  (on  DW  basis; 
Q10=2) 

RDS(T) 

0.047  d1 

1 

Site  information 

Water  depth  (field  site) 

DPT(T) 

1.3  m 

User  def. 

Daily  water  temperature  (field  site) 

WTMP(T) 

°C 

User  def. 

Water  type  specific  light  extinction  coefficient  (field 
site) 

L(T) 

1.07  m  1 

1 

Water  type  specific  current  velocity  (field  site) 

WVEL 

0-100  cm  s1 

User  def. 

Total  live  dry  weight  measured  (field  site) 

TGWM(T) 

-,  g  DM  nr2 

User  def. 

Tuber  density  measured  (field  site) 

NTM(T) 

400  m  2 

3 

Harvesting 

Harvesting 

HAR 

0  or  1 

User  def. 

Harvesting  day  number 

HARDAY 

1-365 

User  def. 

Harvesting  depth  (measured  from  water  surface; 

1-5  m) 

HARDEP 

0.1  m<DEPTH 

User  def. 

*  Calibration  function. 

I.  Best  and  Boyd  2003c.  2.  Sher-Kaul  etal.  1995.  3.  Van  Wijk  1989.  4.  Howard-Williams  1978.  5.  Van  Wijk  1988. 

6.  Best  and  Boyd  1996.  7.  Spencer  1987.  8.  Spencer  and  Anderson  1987.  9.  Best  etal.  2005.  10.  Van  der  Bijl  et  ai.  1989. 

II.  Penning  de  Vries  and  Van  Laar  1982a,  b.  12.  Van  Wijk  et  al.  1988. 


Table  2c.  Parameter  values  used  in  HYDRIL. 


Parameter 

Name 

Value 

Reference 

Morphology,  phenological  cycle  and  development 

Fraction  of  total  dry  matter  increase  allocated  to 
leaves 

FLV(T) 

0.34 

1,  2,  3 

Fraction  of  total  dry  matter  increase  allocated  to 
stems 

FST(T) 

0.60 

1,  2,  3 

Fraction  of  total  dry  matter  increase  allocated  to 
roots 

FRT(T) 

0.06 

1,  2,  3 

Development  rate  as  function  of  temperature* 

DVRV 

0.012  d1 

Calibrated 

DVR  prior  to  flowering  DVRV,  subsequently  DVRR 

DVRR 

0.012  d1 

Calibrated 

Plant  density  and  maximum  plant  biomass 

Plant  density 

NPL 

35  m-2 

4,  5 

Maximum  plant  biomass 

900  g  DW  m  2 

4 

Wintering,  sprouting,  and  growth  of  sprouts  to  water  surface 

Dormant  tuber  density 

NT 

500  m  2 

4 

Tuber  size 

INTUB 

0.1  g  DW  tuber1 

6 

Rel.  tuber  death  rate  (on  number  basis) 

RDTU 

0.36  d1 

7 

Rel.  conversion  rate  of  tuber  into  plant  material 

ROC 

0.0576  g  CH20. 
g 1  DW  d-1 

8 

Relation  coefficient  tuber  weight-stem  length 

RCSHST 

12  m  g  DW-1 

3,  8 
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Light,  photosynthesis,  maintenance,  growth  and  assimilate  partitioning 

Plant  species  specific  light  extinction  coefficient 

K 

0.01  m2  g 1  DW 

4,  9 

Fraction  of  irradiation  shaded  by  epiphytes* 

EPISHD 

0-0.43 

10 

Potential  CO2  assimilation  rate  at  light  saturation  for 
shoots 

AMX 

0.0158  g  CO2. 
g1  DW  h  1 

4,  11 

Initial  light  use  efficiency  for  shoots 

EE 

0.000011  gC02  J1 

12 

Rel.  reduction  factor  for  AMX  to  account  for  senes¬ 
cence  plant  parts 

REDF(T) 

1.0 

User  def. 

Reduction  factor  to  relate  AMX  to  water  pH 

RE  DAM 

0.581 

4,  13 

Daytime  temperature  effect  on  AMX  as  function  of 
DVS* 

AMTMP(T) 

0-1 

13 

Rel.  factor  to  relate  AMX  to  water  current  velocity* 

REDAM1 

0-1 

10 

Dry  matter  allocation  to  each  plant  layer* 

DMPC 

0-1 

4,  14 

Flowering,  translocation,  senescence,  and  formation  of  wintering  organs 

Max.  relative  tuber  growth  rate  (part  of  net  photosyn¬ 
thetic  rate) 

RTR 

0.4  d1 

15 

Conversion  factor  for  translocated  dry  matter  into 

CH2O 

CVT 

1.1 

12 

Tuber  number  concurrently  initiated  per  plant 

NINTUB 

7  plant1 

4 

Total  critical  weight  new  tubers 

TWCTUB 

24.5  g  DW  nr2 

4,  6 

Rel.  death  rate  of  leaves  (on  DW  basis) 

RDR(T) 

0.033  d1 

8 

Rel.  death  rate  of  stems  and  roots  (on  DW  basis) 

RDS(T) 

0.033  d1 

8 

Site  information 

Water  depth  (field  site) 

DPT(T) 

1.0  m 

User  def. 

Daily  water  temperature  (field  site) 

WTMP(T) 

-,  °C 

User  def. 

Water  type  specific  light  extinction  coefficient  (field 
site) 

L(T) 

0.83  m  1 

1,4 

Water  type  specific  current  velocity  (field  site) 

WVEL 

0-100  cm  s1 

User  def. 

Total  live  dry  weight  measured  (field  site) 

TGWM(T) 

-,  g  DW  nr2 

User  def. 

Tuber  bank  density  measured  (field  site) 

NTM(T) 

58-510  m  2 

8 

Harvesting 

Harvesting 

HAR 

0  or  1 

User  def. 

Harvesting  day  number 

HARDAY 

1-365 

User  def. 

Harvesting  depth  (measured  from  water  surface; 

1-5  m) 

HARDEP 

0.1  m<DEPTH 

User  def. 

*  Calibration  function. 


1.  Haller  and  Sutton  1975.  2.  Van  et  al.  1978b.  3.  Van  der  Zweerde  1981.  4.  Bowes  et  al.  1979.  5.  Barko  and  Smart  1981. 
6.  Van  etal.  1977.  7.  Sutton,  pers.  comm.  1995.  8.  Bowes  etal.  1977.  9.  Ikusima  1970.  10.  Best  etal.  2005.  11.  Vanetal. 
1978a.  12.  Penning  de  Vries  and  Van  Laar  1982a,  b.  13.  Van  et  al.  1976.  14.  Ambasht  and  Ram  1976. 15.  Haller  et  al.  1976. 
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Table  2d.  Parameter  values  used  in  MILFO. 


Parameter 

Name 

Value 

Reference 

Morphology,  phenological  cycle  and  development 

Fraction  of  total  dry  matter  increase  allocated  to  leaves 

FLV(T) 

0.47 

1 

Fraction  of  total  dry  matter  increase  allocated  to  roots 

FRT(T) 

0.06 

1 

Fraction  of  total  dry  matter  increase  allocated  to  stems 

FST(T) 

0.47 

1 

Development  rate  as  function  of  temperature* 

DVRV 

0.015  d1 

Calibrated 

DVR  prior  to  flowering  DVRV,  subsequently  DVRR 

DVRR 

0.022  d1 

Calibrated 

Maximum  plant  density  and  plant  biomass 

Plant  density 

NPL 

11  nr2 

2 

Maximum  plant  biomass 

2,283  g  DW  m  2 

2 

Wintering  and  sprouting  of  rhizomes/root  crowns,  and  growth  of  sprouts  to  water  surface 

Initial  rhizome  weight 

IWGRIZ 

50  g  DW  nr2 

3 

Rel.  rhizome  death  rate 

RDRIZ 

0.00042  d1 

4 

Rel.  conversion  rate  of  rhizome  weight  into  plant  material 

ROC 

0.0576  gCH20. 
g1  DW  d-1 

5 

Relation  coefficient  rhizome/ root  crown  weight-stem  length 

RCSHST 

12  m  g  DW-1 

5,  6 

Critical  rhizome  weight 

CRRIZ 

12  g  DW  m-2 

7 

Light,  photosynthesis,  maintenance,  growth  and  assimilate  partitioning 

Plant  species  specific  light  extinction  coefficient 

K(T) 

0.006  m2  g1  DW 

8 

Fraction  of  irradiation  shaded  by  epiphytes* 

EPISHD 

0-0.43 

9 

Potential  CO2  assimilation  rate  at  light  saturation  for  shoot 
tips 

AMX 

0.0165  g  C02. 
g-i  DW  h1 

10,  11 

Initial  light  use  efficiency  for  shoot  tips 

EE 

0.000011  g  CO2 

J'1 

12 

Rel.  reduction  factor  for  AMX  to  account  for  senescence 
plant  parts 

REDF(T) 

1.0 

User  def. 

Reduction  factor  to  relate  AMX  to  water  pH 

REDAM 

0.5 

13,  14 

Daytime  temperature  effect  on  AMX  as  function  of  DVS 

AMTMP(T) 

0-1 

8,  15 

Rel.  factor  to  relate  AMX  to  water  current  velocity* 

REDAM1 

0-0.43 

9 

Dry  matter  allocation  to  each  plant  layer 

DMPC(T) 

0-1 

1 

Flowering,  translocation,  senescence,  and  formation  of  wintering  organs 

Translocation  (part  of  net  photosynthetic  rate) 

TRAFAC 

0.35 

5 

Conversion  factor  for  translocated  dry  matter  into  CH2O 

CVT 

1.05 

12 

Rel.  death  rate  of  leaves  (on  DW  basis;  Q10=2) 

RDR(T) 

0.042  d1 

10 

Rel.  death  rate  of  stems  and  roots  (on  DW  basis;  Q10=2) 

RDS(T) 

0.042  d1 

10 

Site  information 

Water  depth  (field  site) 

DPTT(T) 

1.5  m 

User  def. 

Daily  water  temperature  (field  site) 

WTMP(T) 

-,  °C 

User  def. 

Water  type  specific  light  extinction  coefficient  (field  site) 

L(T) 

1.15-2.00  m1 

13 

Water  type  specific  current  velocity  (field  site) 

WVEL 

0-100  cm  s1 

User  def. 

Total  live  dry  weight  measured  (field  site) 

TGWM(T) 

-,  g  DW  rrr2 

User  def. 
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Harvesting 

Harvesting 

HAR 

0  or  1 

User  def. 

Harvesting  day  number 

HARDAY 

1-365 

User  def. 

Harvesting  depth  (measured  from  water  surface  in  m) 

HARDEP 

0.1  m<DEPTH 

User  def. 

*  Calibration  function. 

I.  Adams  et  al.  1974.  2.  Budd  et  al.  1995.  3.  Smith  and  Adams  1986.  4.  Vogt  et  al.  1991.  5.  Bowes  et  al.  1979.  6.  Van  der 
Zweerde  1981.  7.  Madsen  1997.  8.  Titus  and  Adams  1979a.  9.  Bestetal.  2005.  10.  Adams  and  McCracken  1974. 

II.  Van  et  al.  1976.  12.  Penning  de  Vries  and  Van  Laar.  1982a,  b.  13.  Lee  and  Kluesener  1972.  14.  Titus  and  Stone  1982. 
15.  Stanley  and  Nailor  1972.  16.  Kooman  1995. 


Morphology,  phenological  cycle,  and  development 

Morphology  and  phenological  cycle. 

•  V.  americana.  V.  americana  plants  have  a  short  stem  axis,  bearing  a 
variable  number  of  stolons.  The  plant  forms  a  rosette  on  top  of  the 
sediment  (“meadow  formation”)  with  up  to  1.2-m-long,  linear  tape¬ 
shaped  leaves,  submersed  or  floating  at  the  water  surface  (Lowden 
1982;  Titus  and  Stephens  1983;  Korschgen  and  Green  1988).  Typical 
partitioning  of  biomass— 71.8  percent  in  leaves,  15.9  percent  in  stems, 
and  12.3  percent  in  roots— is  based  on  summer  habit  of  full-grown 
plants  (Titus  and  Stephens  1983;  Haller  1974;  FLVT,  FSTT,  FRTT). 

This  species  propagates  largely  vegetatively  by  tubers  in  North  Amer¬ 
ica.  The  production  of  viable  seeds  is  low  and  seedling  establishment  in 
the  field  is  rare  in  the  Northern  U.S.  (Choudhuri  1966;  Titus  and 
Stephens  1983;  Kimber  et  al.  1995).  Plants  flower  once  a  year,  from 
late-June  to  August  in  the  northern  hemisphere  (Donnermeyer  1982; 
Titus  and  Stephens  1983;  Donnermeyer  and  Smart  1985). 

•  P.  pectinatus.  P.  pectinatus  plants  have  long,  flexible  stems  with  nar¬ 
row,  mostly  filiform  leaves  near  or  at  the  water  surface  (“canopy  forma¬ 
tion”)  and  emergent  flowers  (Den  Hartog  1982).  Typical  partitioning  of 
biomass— 73.1  percent  in  leaves,  18.3  percent  in  stems,  and  8.6  percent 
in  roots— is  based  on  summer  habit  of  full-grown  plants  (Best  and  Boyd 
2003c;  Sher-Kaul  et  al.  1995;  FLVT,  FSTT,  FRTT ).  In  shallow  waters, 
this  species  usually  behaves  as  an  annual  because  it  is  sensitive  to  frost 
and  decomposes  rapidly  in  water  (Lohammar  1938;  Lapirov  and 
Petukhova  1985).  However,  in  deeper  waters,  green  perennial  shoots 
can  be  found  (Hammer  and  Heseltine  1988).  This  species  propagates 
largely  by  tubers.  Regeneration  from  rhizomes,  turions  or  seeds  may 
occur  also.  Plants  flower  once  a  year,  from  late-June  to  August  in  the 
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northern  hemisphere  (Yeo  1965).  In  mild  climates  flowering  may  ex¬ 
tend  up  to  five  months  (Gupta  1968;  Ramirez  and  San  Martin  1984). 

•  H.  verticillata.  H.  verticillata  plants  have  long,  flexible,  stems,  that 
branch  close  to  the  water  surface  (“canopy  formation”)  and  are  sur¬ 
rounded  by  whorls  of  leaves.  Typical  partitioning  of  biomass— 34  per¬ 
cent  in  leaves,  60  percent  in  stems,  and  6  percent  in  roots— is  based  on 
summer  habit  of  full-grown  plants  (Haller  and  Sutton  1975;  Van  et  al. 
1978b;  Van  der  Zweerde  1981;  FLVT,  FSTT,  FRTT).  H.  verticillata  bio¬ 
types  regenerate  largely  from  tubers,  but  may  have  stem  fragmenta¬ 
tion,  horizontal  shoot  formation  in  the  sediment,  and  turion  formation 
on  the  shoots  in  addition  to  this  vegetative  reproduction  mechanism. 
Plants  flower  once  a  year,  i.e.,  from  early  June  to  the  end  of  August  in 
near-subtropical  Florida  (Bowes  et  al.  1979). 

•  M.  spicatum.  M.  spicatum  plants  have  long,  flexible  stems,  that 
branch  close  to  the  water  surface  (“canopy  formation”)  and  are  sur¬ 
rounded  by  whorls  of  leaves.  Typical  distribution  of  biomass— 47  per¬ 
cent  in  leaves,  47  percent  in  stems,  and  6  percent  in  roots— is  based  on 
summer  habit  of  full-grown  plants  (Adams  et  al.  1974;  FLVT,  FSTT, 
FRTT).  This  species  is  able  to  propagate  in  an  evergreen  condition,  by 
vegetative  fragmentation,  and  by  seeds.  It  most  frequently  winters  in 
an  evergreen  form,  as  root  crowns/rhizome  system  (Grace  and  Wetzel 
1978;  Madsen  et  al.  1988;  Madsen  1997)  and  may  maintain  consider¬ 
able  winter  biomass  (Stanley  et  al.  1976).  Plants  flower  from  June  to 
November  in  the  northern  hemisphere;  one  (Aiken  et  al.  1979;  Carpen¬ 
ter  1980),  two  (Nichols  1971;  Lind  and  Cottam  1969)  and  three  (Grace 
and  Wetzel  1978)  flowering  periods  per  year  have  been  reported.  The 
timing  of  flowering  is  “less  predictable”  in  more  southern  areas  (Grace 
and  Wetzel  1978),  whereas  it  may  last  the  whole  growth  season  in  the 
tropics  (Zutschi  and  Vass  1973). 

Phenology  and  development.  Plant  development  rates  were  derived  from 

field  data,  typical  for  temperate  regions.  They  are  all  in  the  same  order  of 

magnitude  (Tables  3a,  b,  c,  d,  e). 

•  V.  americana :  0.015  d_1  prior  to  the  flowering  period  and  0.040  d1 
subsequently,  at  a  reference  temperature  of  30  °C  and  a  temperature 
threshold  of  3°C. 

•  P.  pectinatus:  0.015  d1  prior  to  the  flowering  period  and  0.040  d1  sub¬ 
sequently. 
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Table  3a.  Relationship  between  development  phase  (DVS)  of  V.  americana,  day  of  year,  at  a  reference 
temperature  of  30  °C  and  3  °C  day-degree  sum  in  a  temperate  climate  (DVRVT=  0.015;  DVRRT=  0.040). 


Developmental  Phase 

Day  Number 

3  °C  Day-Degree 
Sum 

Description 

DVS  Value 

First  Julian  day  number  *  tuber  sprouting  and 
initiation  elongation 

0  ■*  0.291 

0  -*  105 

1  *  270 

Tuber  sprouting  and  initial  elongation  ->  Leaf 
expansion 

0.292  -*  0.875 

106  ■*  180 

271  -*1215 

Leaf  expansion  ->  floral  initiation  and  anthesis 

0.876  -*1.000 

181  -*191 

1216  -*  1415 

Floral  initiation  and  anthesis  ->  induction  of  tu¬ 
ber  formation,  tuber  formation  and  senescence 

1.001  -*  2.000 

192  ■*  227 

1416  -*  2072 

Tuber  formation  and  senescence  -*  senesced 

2.001  ->  4.008 

228  -*  365 

2073  ■*  3167 

Senesced 

4.008 

365 

3167 

Note:  Calibration  was  on  field  data  on  biomass  and  water  transparency  from  Chenango  Lake,  NY,  1978  (Titus  and  Stephens 
1983)  and  climatological  data  from  Binghamton  (air  temperatures)  and  Ithaca  (irradiance),  NY,  1978. 


Table  3b.  Relationship  between  development  phase  (DVS)  of  P.  pectinatus,  day  of  year,  at  a  reference 
temperature  of  30  °C  and  3  °C  day-degree  sum  in  a  temperate  climate  (DVRVT=  0.015;  DVRRT=  0.040). 


Developmental  Phase 

Day  Number 

3  °C  Day-Degree 
Sum 

Description 

DVS  Value 

First  Julian  day  number  *  tuber  sprouting  and 
initiation  elongation 

0  -*  0.210 

0  -*  77 

1  ■*  193 

Tuber  sprouting  and  initial  elongation  ■*  Leaf 
expansion 

0.211  -*  0.929 

78  ■*  187 

194  ■*  1301 

Leaf  expansion  ■*  floral  initiation  and  anthesis 

0.930  -*  1.000 

188  ■*  195 

1302  ■*  1434 

Floral  initiation  and  anthesis  *  induction  of  tu¬ 
ber  formation,  tuber  formation  and  senescence 

1.001  -*  2.000 

196  ■*  233 

1435  *  2077 

Tuber  formation  and  senescence  -*  senesced 

2.001  ■*  4.033 

234  -*  365 

2078  -*  3193 

Senesced 

4.033 

365 

3193 

Note:  Calibration  was  on  field  data  on  biomass  and  water  transparency  from  the  Western  Canal  near  Zandvoort,  NL,  1987 
(Best  and  Boyd  2003c)  and  climatological  data  from  De  Bilt,  NL,  1987. 
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Table  3c.  Relationship  between  development  phase  (DVS)  of  H.  verticillata,  day  of  year,  at  a  reference 
temperature  of  30  °C  and  3  °C  day-degree  sum  in  a  temperate  climate  (DVRVT=  0.012;  DVRRT=  0.012). 


Developmental  Phase 

Day  Number 

3  °C  Day-Degree 
Sum 

Description 

DVS  value 

First  Julian  day  number  ->  tuber  sprouting  and 
initiation  elongation 

0  -»  0.326 

0  ^  74 

1  ->  596 

Tuber  sprouting  and  initial  elongation  ->  Leaf 
expansion 

0.327  ->  0.360 

75  -»  78 

597  -»  673 

Leaf  expansion  ->  floral  initiation  and  anthesis 

0.361  ->  1.000 

79  ->  150 

674  2063 

Floral  initiation  and  anthesis  ->  induction  of  tu¬ 
ber  formation 

1.001  2.000 

151  ->  241 

2064  ->  4279 

Induction  tuber  formation  ->  tuber  formation  and 

senescence 

2.001  ->  2.200 

242  -»  260 

4280  4744 

Tuber  formation  and  senescence  ->  senesced 

2.201  ->  2.979 

261  -»  365 

4745  6357 

Senesced 

2.979 

365 

6357 

Note:  Calibration  was  on  field  data  on  biomass  and  water  transparency  from  Lake  Orange,  FL,  1977  (Bowes  et  al.  1979)  and 
climatological  data  from  Gainesville,  FL,  1980.  Sprouting  was  timed  after  Hailer  et  al.  (1976). 


Table  3d.  Relationship  between  development  phase  (DVS)  of  M.  spicatum,  day  of  year,  at  a  reference 
temperature  of  30  °C  and  3  °C  day-degree  sum  in  a  temperate  climate  (DVRVT=  0.022;  DVRRT=  0.015). 


Developmental  Phase 

Day  Number 

3  °C  Day-Degree 
Sum 

Description 

DVS  Value 

First  Julian  day  number  ->  sprouting,  initiation 
elongation,  and  leaf  expansion  C0H0RT1 

0  -»  0.375 

0  114 

1  ->  191 

Sprouting,  initiation  elongation,  and  leaf  expan¬ 
sion  ->  floral  initiation,  anthesis,  and  induction  of 
senescence  C0H0RT1 

0.376  ->  1.000 

115  162 

192  ->  900 

Floral  initiation,  anthesis,  and  induction  of  se¬ 
nescence  ->  senescence  C0H0RT1 

1.001  1.630 

163  ->  212 

901  ^  2012 

Senescence  ->  senesced  C0H0RT1 

1.631  ->  2.000 

213  -»  245 

2013  2669 

Sprouting,  initiation  elongation,  and  leaf  expan¬ 
sion  ->  floral  initiation,  anthesis,  and  induction  of 
senescence  C0H0RT2 

1.001  ->  1.630 

163  212 

901  ->  2012 

Floral  initiation,  anthesis,  and  induction  of  se¬ 
nescence  ->  anthesis  and  senescence  C0FI0RT2 

1.631  ->  2.000 

213  245 

2013  ^  2669 

Senescence  ->  senesced  C0FI0RT2 

2.001  ->  2.570 

246  -»  365 

2670  ->  3508 

Senesced  COFIORT  1  and  2 

2.570 

365 

3508 

Note:  Calibration  was  on  field  data  on  biomass  (Adams  and  McCracken  1974),  and  on  water  transparency,  temperature,  and 
irradiance  from  Lake  Wingra,  Wl,  1970  (Lee  and  Kluesener  1972). 
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Table  3e.  Relationship  between  development  phase  (DVS)  of  M.  spicatum,  day  of  year,  at  a  reference 
temperature  of  30  °C  and  3°C  day-degree  sum  in  a  tropical  climate  (DVRVT=  0.022;  DVRRT=  0.015). 


Developmental  Phase 

Day  Number 

3  °C  Day-Degree 
Sum 

Description 

DVS  Value 

First  Julian  day  number  A  sprouting,  initiation  elon¬ 
gation,  and  leaf  expansion  COHORT  1 

0  A  0.375 

0  A  25 

1  A  431 

Sprouting,  initiation  elongation,  and  leaf  expansion 

A  floral  initiation,  anthesis,  and  induction  of  senes¬ 
cence  COHORT  1 

0.376  ->  1.000 

115  A  61 

432  A  1163 

Floral  initiation,  anthesis,  and  induction  of  senes¬ 
cence  ->  senescence  COHORT  1 

1.001  ->  1.630 

62  A  162 

1164  ->  3844 

Senescence  A  senesced  COHORT  1 

1.631  A  2.000 

163  A  188 

3845  A  4490 

Sprouting,  initiation  elongation,  and  leaf  expansion 

A  floral  initiation,  anthesis,  and  induction  of  senes¬ 
cence  COHORT  2 

1.001  A  1.630 

62  A  162 

1164  ->  3844 

Floral  initiation,  anthesis,  and  induction  of  senes¬ 
cence  ->  anthesis  and  senescence  COHORT  2 

1.631  A  2.000 

163  A  188 

3845  A  4490 

Senescence  A  senesced  COHORT  2 

2.001  A  2.570 

164  A  233 

4491  A  5492 

Sprouting,  initiation  elongation,  and  leaf  expansion 

A  floral  initiation,  anthesis,  and  induction  of  senes¬ 
cence  COHORT  3 

2.001  A  2.447 

164  A  223 

4491  A  5273 

Floral  initiation,  anthesis,  and  induction  of  senes¬ 
cence  A  anthesis  and  senescence  COHORT  3 

2.448  A  3.500 

224  A  307 

5274  A  7125 

Senescence  A  senesced  COHORT  3 

3.501  ->  4.141 

308  A  365 

7126  A  8254 

Senesced  COHORT  1,  2,  and  3 

4.141 

365 

8254 

Note:  Calibration  was  on  field  data  on  biomass  from  Kashmir  lakes,  India,  1970s  (Zutschi  and  Vass  1973),  and  on  climatologi¬ 
cal  data  from  Patancheru,  India,  1978. 


•  H.  verticillata:  0.012  d1  prior  to  the  flowering  period  and  0.012  d_1 
subsequently. 

•  M.  spicatum :  0.022  d1  prior  to  the  flowering  period  and  0.015  d1  sub¬ 
sequently. 

Each  simulation  starts  at  the  first  Julian  day  number,  i.e.,  January  1,  when 
the  development  phase  has  the  value  0.0  (Table  3).  The  quantities  of 
leaves,  stems,  and  roots  are  set  equal  to  o.  For  the  tuber-producing  species 
(Tables  2a,  2b,  2c),  the  simulation  starts  using  an  observed  tuber  bank 
density,  with  a  certain,  chosen  (see  “Wintering,  sprouting  in  spring,  and 
growth  of  sprouts  to  water  surface”)  individual  tuber  weight  as  initial 
values. 


For  M.  spicatum  (Tables  2d,  3e),  the  simulation  starts  using  an  observed 
rhizome  weight  at  the  beginning  of  the  year.  If  simulation  of  a  SAV  at  a 
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different  site  is  desired,  the  simulation  may  start  also  with  wintering 
plants  present;  first,  however,  initial  quantities  of  plant  organs  must  be 
calculated.  Initiation  of  growth  activity  for  a  tuber-producing  SAV  in  a 
temperate  climate  occurs  by  sprouting  of  the  tubers,  i.e.,  occurs  at  a  DVS 
between  0.2  and  0.3  (Tables  2a,  2b,  2c).  Sprouts  of  the  first  plant  cohort 
develop  through  remobilization  of  carbohydrates  from  tubers.  The  sprouts 
elongate  up  to  the  preset  maximum  plant  height  (or  to  the  water  surface  in 
cases  where  water  depth  is  less  than  the  maximum  preset  plant  height), 
and  subsequently  follow  a  species-characteristic  vertical  distribution  over 
the  water  column.  Anthesis  is  initiated  at  DVS  1.000  and  finishes  at  DVS 
2.000  just  before  new  tubers  are  initiated.  Tubers  can  be  formed  directly 
at  initiation  in  V.  americana  and  P.  pectinatus,  or  after  a  certain  lag 
period  in  H.  verticillata.  Tuber  formation,  downward  translocation,  and 
senescence  set  in  at  DVS  2.001  (2.201  in  H.  verticillata )  and  continue  until 
the  end  of  the  year  (Tables  3a,  3b,  3c). 

Initiation  of  growth  activity  for  a  M.  spicatum  community  in  a  temperate 
climate  (Table  2d)  occurs  by  sprouting  of  the  rhizome/root  crown  system, 
i.e.,  at  DVS  0.376.  Sprouts  of  plant  COHORT  1  develop  through  remobili¬ 
zation  of  carbohydrates  from  the  rhizome/root  crown  system.  The  sprouts 
elongate  up  to  the  water  surface  and  form  a  canopy  in  the  upper  water  lay¬ 
ers.  Anthesis  of  COHORT  1  is  initiated  at  DVS  1.000  and  finishes  at  DVS 
1.630,  just  before  downward  carbohydrate  translocation  and  senescence 
are  initiated.  Translocation  and  senescence  of  COHORT  1  set  in  at  DVS 
1.631  and  continue  until  DVS  2.000.  Sprouting  of  COHORT  2  starts  when 
translocation  and  senescence  of  COHORT  1  have  set  in.  This  timing  is 
based  on  the  assumption  that  at  that  time,  apical  dominance  by  the  exist¬ 
ing,  senescing  shoots  is  broken  and,  consequently,  new  shoots  develop. 
Sprouting  of  COHORT  2  occurs  from  growing  points  on  the  rhizome/root 
crown  system.  Anthesis  of  COHORT  2  is  initiated  at  DVS  1.631  and  fin¬ 
ishes  at  DVS  2.000.  Translocation  and  senescence  of  COHORT  2  set  in 
at  DVS  2.001  and  continue  until  the  end  of  the  year.  For  a  vegetation  in 
the  tropics  (Table  3e),  it  proved  impossible  to  generate  the  high  levels  of 
shoot  and  rhizome/root  crown  biomass  reported  (Zutschi  and  Vass  1973) 
with  two  plant  cohorts  active,  since  the  second  plant  cohort  had  already 
senesced  in  May.  However,  proper  biomass  levels  and  timing  were  at¬ 
tained  with  three  plant  cohorts  active,  the  third  cohort  being  switched  on 
at  latitudes  less  than  33  °N.  It  is  possible  that  a  particular  plant  process, 
such  as  sprouting,  is  sensitive  to  day  length,  and  that  this  process  decides 
for  the  population  to  activate  another  cohort.  However,  since  the  authors 
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are  not  aware  of  publications  on  this  topic  for  M.  spicatum,  the  switch  has 
been  set  at  the  cutoff  latitude  for  tropical  areas. 

Plants  in  tropical  regions  behave  similarly  in  terms  of  DVS  to  those  in  tem¬ 
perate  regions,  except  that  tropical  plants  require  on  average  a  1.6  x  higher 
3  °C  day-degree  sum  to  complete  their  individual  life  cycle  than  temperate 
cohorts. 

Plant  density  and  maximum  biomass 

V.  americana.  For  V.  americana,  typical  plant  density  of  30  plants  per 
square  meter  (NPL)  was  calculated  by  dividing  the  maximum  standing 
crop  of  an  established,  monotypic  vegetation  (50.1  g  DW  nr2)  by  the  aver¬ 
age  weight  of  an  individual  rosette  with  neighbor  plants  (1.65  g  DW  nr2)  in 
Chenango  Lake,  NY  (Titus  and  Stephens  1983).  Maximum  biomass  ranged 
from  253  g  DW  nr2  at  a  0.73-m  rooting  depth  in  Lake  Biwa,  Japan,  in  1969 
(lat  350  30'  N;  Ikusima  1970);  to  344  g  DW  nr2  at  a  1.2-1.5-m  rooting  depth 
in  University  Bay,  WI,  in  1976  (lat  430  08'  N;  Titus  and  Adams  1979b);  to 
496  g  DW  nr2  in  1.5-m  deep  experimental  ponds,  in  Orange  County,  FL 
(lat  28°  30'  N;  Haller  1974).  For  subtropical  to  tropical  areas  in  Florida,  a 
relatively  constant,  high  biomass  has  been  suggested  (Godfrey  and  Wooten 
1997).  The  highest  published  value  on  maximum  biomass  (<496  g  DW  nr 
2)  was  selected  as  upper  limit  of  plant  biomass  in  VALLA. 

P.  pectinatus.  For  P.  pectinatus,  typical  plant  density  of  30  plants  per 
square  meter  (NPL)  was  computed  by  dividing  the  maximum  standing 
crop  of  an  established,  monotypic  vegetation  at  a  2.5-m  rooting  depth 
(82.27  g  DW  nr2;  Best  and  Boyd  2003c)  by  the  highest  average  weight  of 
an  individual  plant  with  neighbor  plants  (2.76  g  DW  plant1)  in  shallow 
waters  studied  in  NL  and  The  Camargue,  France  (Best  and  Boyd  2003c; 
Van  Wijk  1989).  Other  literature  reviewed  did  not  provide  sufficiently  de¬ 
tailed  information  to  enable  calculations  of  plant  density.  In  a  temperate 
climate,  the  highest  biomass  was  found  in  brackish  and  saline  lakes,  i.e., 
1,313  g  DW  nr2  in  a  brackish  pool  near  Yerseke,  NL  (lat  510  30'  N;  Van 
Wijk  et  al.  1988);  1,568  g  DW  nr2  in  saline  Lake  Mariut,  Egypt  (lat  310 
io'  N;  Aleem  and  Samaan  1969);  and  1,952  g  DW  nr2  in  the  saline 
Swartvlei,  between  Port  Elizabeth  and  Cape  Town,  South  Africa  (lat  340 
o'  S;  Howard-Williams  1978).  A  somewhat  lower  maximum  of  712  g 
DW  nr2  occurred  in  freshwater  Badfish  Creek,  WI  (lat  43°4'  N;  Madsen 
and  Adams  1988).  Maximum  biomass  in  freshwater  lakes  in  the  tropics  is 
similar  to  that  in  a  temperate  climate,  i.e.,  370-445  g  DW  nr2  (lat  170  27'  N; 
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Sinha  1970;  Sahai  and  Sinha  1973).  The  highest  published  value  on  maxi¬ 
mum  biomass  (<1,952  g  DW  nr2)  was  selected  as  upper  limit  of  plant  bio¬ 
mass  in  POTAM. 

H.  verticillata.  For  H.  verticillata,  typical  plant  density  of  35  plants  per 
square  meter(NPL)  was  considered  as  representative.  A  range  of  34-53 
plants  nr2  was  calculated  by  dividing  the  standing  crop  of  an  established 
vegetation  (163  g  DW  nr2)  in  Lake  Orange  in  August  (1977;  Bowes  et  al. 
1979)  by  the  weight  range  of  individual  plants  (3.08-4.76  g  DW  plant1; 
Barko  and  Smart  1981).  Maximum  biomass  varied  from  163  g  DW  nr2  in 
Lake  Orange,  FL,  in  1977  (Bowes  et  al.,  1979);  to  500  g  DW  nr2  in  Lake 
Ramgarh,  India,  in  1967-68  (Sahai  and  Sinha  1973);  to  890  g  DW  nr2  in 
Lake  Trafford  (shoots;  Bowes  et  al.  1979). 

The  highest  published  value  on  maximum  biomass  (i.e.,  890  g  shoot  bio¬ 
mass  augmented  with  10  g  root  biomass  nr2)  of  900  g  DW  nr2)  was  se¬ 
lected  as  upper  limit  of  plant  biomass  in  HYDRIL. 

M.  spicatum.  For  M.  spicatum,  typical  plant  density  of  11  plants  per 
square  meter  (NPL)  was  the  average  within  a  range  of  3  to  32  “plant 
clumps”  per  square  meter,  consisting  of  a  variable  number  of  stems, 
measured  in  a  vegetation  in  Fish  Lake,  WI,  in  1990-92  (Budd  et  al.  1995). 
Maximum  biomass  was  2,283  g  DW  nr2  in  the  same  lake  and  year,  and 
similar  values  have  been  reported  for  the  more  southern  Lake  Gunters- 
ville,  AL,  in  1972  (Stanley  et  al.  1976).  The  latter  value  on  maximum  bio¬ 
mass  (<2,283  g  DW  nr2)  was  selected  as  upper  limit  of  plant  biomass  in 
MILFO. 

Wintering,  sprouting  in  spring,  and  growth  of  sprouts  to  water  surface 

V.  americana.  Tuber  bank  density  in  this  species  varies  over  a  large  range, 
from  o  in  early  summer  to  450  nr2  in  autumn.  The  large  range  found  may 
be  attributed  partly  to  the  patchy  spatial  distribution  of  the  community 
over  the  water  body,  limited  number  of  replicate  samples  taken  (Spencer 
et  al.  1994),  and  between-site  variation  in  rooting  depth  of  the  vegetation. 
Densities  published  are  (1)  101  nr2  in  Lake  Mendota,  WI  (Titus  and  Adams 
1979b),  (2)  170  nr2  in  Pool  9  of  the  Upper  Mississippi  River,  WI  (Donner- 
meyer  1982);  (3)  233  nr2  in  the  Lower  Detroit  River,  MI  (Korschgen  and 
Green  1988);  (4)  330-450  nr2  in  Chenango  Lake,  NY  (Titus  and  Stephens 
1983);  and  (5)  115-1,140  nr2  in  Lake  Mattamuskeet,  NC  (Loworn  and 
Gillingham  1996).  A  tuber  bank  density  of  233  nr2  with  a  tuber  size  of 
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0.09  g  DW  found  in  well-established  tuber  banks  (Korschgen  and  Green 
1988)  was  used  for  calibration. 

Tuber  sizes  published  are,  in  g  DW  tuber1,  (1)  0.04  to  0.18  (various  north- 
american  water  bodies;  Korschgen  and  Green  1988;  Korschgen  et  al. 

1997);  (2)  0.055  (average,  Chenango  Lake,  NY,  1.4-m  rooting  depth;  Titus 
and  Stephens  1983);  (3)  0.07  (average,  Lake  George,  NY,  2-  to  4-m  rooting 
depth,  J.  D.  Madsen  unpublished,  Vicksburg,  MS);  and  (4)  0.18  (maxi¬ 
mum,  just  after  tuber  completion,  Pool  9  Upper  Mississippi  River,  WI, 
1.1-m  rooting  depth;  Donnermeyer  1982).  Individual  tuber  weight  de¬ 
creased  almost  linearly  with  decreasing  light  level  (from  0.102  to  <0.01  g 
dry  weight  tuber1),  and  parabolically  with  concurrently  formed  tuber 
number  per  plant  (from  0.120  to  0.015;  Korschgen  et  al.  1997).  The  latter 
phenomenon  forms  the  basis  of  the  hypothesis  by  the  authors  of  this  re¬ 
port  that  V.  americana  follows  an  optimization  strategy  aimed  at  produc¬ 
ing  the  largest  tubers  for  the  light  level  experienced  because  large  tubers 
have  a  greater  survival  value  than  smaller  ones.  Thus,  an  established  vege¬ 
tation  experiencing  a  given  light  level  will  aim  at  producing  one  tuber  size 
class  only,  i.e.,  with  a  size  that  allows  new  plants  to  survive  at  that  site. 
Consequently,  the  differences  in  tuber  size  found  in  tuber  banks  may  be 
attributed  to,  besides  differences  in  environmental  conditions,  differences 
in  age  between  tuber  size  classes  with  tubers  in  the  oldest  class  weighing 
less  because  of  remobilization  processes,  in  the  last-completed  class  being 
full-size,  and  in  the  youngest  class  weighing  less  because  not  yet  finished. 
The  combination  of  5.5  concurrently  initiated  tubers  and  tuber  size  of 
0.090  g  DW  tuber1  was  selected  as  nominal.  The  typical  relationship  be¬ 
tween  tuber  size  and  concurrently  formed  tuber  number  per  plant  has  to 
be  used  to  calculate  the  proper  initial  values  for  simulations  started  from  a 
different  tuber  size  than  nominal  (Table  4). 

Tuber  death  rates  have  not  been  published.  For  use  in  VALLA,  the  value 
for  the  relative  tuber  death  rate  was  determined  by  applying  the  same  dif¬ 
ferential  equation  as  commonly  used  for  simple  exponential  growth  to  de¬ 
scribe  continuous  exponential  decrease  in  tuber  number,  with  a  negative 
specific  decrease  rate  (Thornley  and  Johnson  1990b;  Hunt  1982).  A  rela¬ 
tive  tuber  death  rate  ( RDTU)  of  0.018  d1  (on  number  basis)  was  found  for 
the  vegetation  in  Chenango  Lake,  NY  (Titus  and  Stephens  1983).  Tuber 
loss  from  grazing  is  expected  to  be  relatively  low  in  this  plant  species  be¬ 
cause  of  its  relatively  low  tuber  bank  density  (compared  to  H.  verticillata ) 
which  may  discourage  foraging  by  waterfowl  because  this  species  requires 
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a  relatively  long  search  time  (Loworn  1989;  Loworn  and  Gillingham 
1996). 

Table  4.  Relationship  between  tuber  size  and  concurrently  initiated  tuber  number  in 
V.  americana  and  P.  pectinatus. 


V.  americana 

P.  pectinatus 

Tuber  Size 
(g  DW  tuber1) 

Concurrently  Initiated 
Tuber  Number 
(N  plant1) 

Tuber  Size 
(g  DW  tuber1) 

Concurrently 

Initiated  Tuber 
Number 
(N  plant1) 

0.015 

0.6 

0.002 

0.5 

0.035 

2.0 

0.050 

3.0 

0.070 

4.0 

0.086 

5.0 

0.055 

5.0 

0.090 

5.5 

0.070 

6.0 

0.120 

5.8 

0.155 

12.0 

Sprouting  frequency  of  tubers  in  the  field  has  not  been  published.  Poten¬ 
tial  sprouting  frequency  is  high  in  a  temperate  climate,  being  <80  percent. 
Sprouting  frequency  in  an  established  community  is  probably  not  impor¬ 
tant,  unless  it  is  very  low,  as  long  as  the  typical  density  of  30  plants  per 
square  meter  is  somehow  reached,  since  density  tends  to  play  a  lesser  role 
in  biomass  production  compared  to  space  availability.  Sprouting  frequency 
is  not  affected  by  day  length,  but  is  higher  in  light  than  in  darkness. 
Sprouting  is  prevented  by  temperatures  <  5  °C,  optimal  between  15  and 
25  °C,  and  maximal  at  20  and  25  °C  (J.  D.  Madsen  unpub.,  Vicksburg,  MS), 
and  occurs  early  in  the  growth  season,  i.e.,  mid-April  to  end-of-May  in 
Pool  9  of  the  Upper  Mississippi  River,  WI,  in  1980-81  (Donnermeyer 
1982).  In  VALLA,  sprouting  is  governed  by  DVS  and,  thus,  indirectly  by 
temperature. 

Minimum  weight  required  for  a  tuber  to  sprout  is  0.003  g  DW  tuber1 
(Donnermeyer  and  Smart  1985;  J.  D.  Madsen  unpub.,  Vicksburg,  MS). 

In  VALLA  the  minimum  tuber  weight  has  been  set  at  0.003  g  dry  weight 
tuber1. 
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The  derivations  of  the  parameters  CRIFAC  and  SURPER  in  VALLA  are  de¬ 
scribed  in  this  paragraph;  those  of  ROC  and  RCHST  are  described  under 
H.  verticillata.  The  elongation  potential  of  plantlets  emerging  from  tubers 
is  limited,  i.e.  0.44  m  (range  0.04  to  0.18  g  DW  tuber1;  Korschgen  and 
Green  1988;  Korschgen  et  al.  1997).  Thus,  plants  can  rise  one  0.1-m  depth 
layer  in  the  water  column  after  filling  that  layer  with  a  minimum  of  0.0091 
or  maximum  of  0.041  g  DW.  Elongation  is  modeled  by  filling  each  succes¬ 
sive  water  layer  from  the  sediment  to  a  1.2-m  water  column  with  the 
minimum  required  shoot  biomass  (0.0091  g  DW  plant1,  termed  CRIFAC ). 
Respiration  of  dormant  tubers  is  low.  Measured  respiration  in  dormant 
P.  pectinatus  tubers  was  0.003623  g  C02  g  DW1  h1  at  20  °C  (Best  and 
Boyd  2003c).  Respiration  is  expected  to  decrease  by  a  factor  of  2  with  a 
10  °C-decrease  in  temperature  (Q10).  Thus,  daily  respiration  of  a  dormant 
tuber  at  10  °C  would  be  0.043476  g  C02  g  DW1  d1  (multiplication  by 
Vi  and  24),  and  yield  0.00391  g  C02  tuber1  d1  of  the  nominal  0.09-g-DW 
tuber  size.  A  tuber  of  this  size  would  survive  23  days  (SURPER),  not  taking 
conversion  costs  of  C02  into  dry  weight  into  consideration  (0.00391  g  C02 
tuber1  d1  x  23  days  =0.09  g  dry  weight)  without  net  photosynthesis  taking 
place.  Larger  tubers  have  longer,  and  smaller  tubers  shorter  survival 
periods. 

P.  pectinatus.  Tuber  bank  density  in  this  species  ranges  from  o  in  early 
summer  to  3,975  nr2  in  autumn  (Van  Wijk  1989).  Potential  causes  for  the 
large  range  found  have  been  outlined  under  V.  americana.  The  following 
densities  have  been  published:  (1)  45-115  nr2  on  a  wave-exposed  shallow 
site  within  the  brackish  Baltic  Sea,  Sweden  (Kautsky  1987);  (2)  270-385 
m  2  in  the  shallow,  fresh,  Lake  Veluwe,  NL  (Van  Dijk  et  al.  1992);  (3)  1,330 
and  3,975  nr2  in  an  oligohaline,  sheltered  ditch  Salin  De  Badon  and  shal¬ 
low  pool  Les  Garcines,  respectively,  both  in  the  Camargue,  France  (Van 
Wijk  1988).  A  tuber  bank  density  of  240  nr2  with  a  tuber  size  of  0.083  g 
DW  found  in  well-established  tuber  banks,  in  the  Western  Canal  near 
Zandvoort,  NL  (Best  and  Boyd  2003c)  was  used  for  calibration. 

Tuber  sizes  published  are,  in  g  DW  tuber1,  (1)  0.005  to  0.155  (average 
0.013  g,  laboratory  study  in  CA;  Spencer  and  Anderson  1987);  (2)  0.017 
(derived  from  0.056  g  fresh  weight,  and  DW:FW  ratio  of  0.299+0.034; 
freshwater  canals,  CA;  Spencer  1987);  (3)  0.083  (Best  and  Boyd  2003c). 
Tubers  may  be  formed  in  classes,  as  suggested  by  the  formation  of  10  tuber 
size  classes  with  an  average  fresh  weight  of  0.061  +  0.006  g  tuber1  (range 
0.02  to  0.3  g)  and  18  tuber  size  classes  with  an  average  fresh  weight  of 
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0.017  +  0.003  g  tuber1  (range  0.02  to  0.46  g)  in  the  fresh  Lake  Veluwe 
and  a  brackish  Texel  ditch,  respectively,  both  in  NL  (Hootsmans  1994). 
Tuber  size  decreased  almost  linearly  with  tuber  number  concurrently 
formed  per  plant  (from  0.155  to  0.002;  Spencer  and  Anderson  1987, 

Table  3).  The  combination  of  8.0  concurrently  initiated  tubers  and  tuber 
size  of  0.083  g  DW  tuber1  was  selected  as  nominal.  The  typical  relation¬ 
ship  between  tuber  size  and  concurrently  formed  tuber  number  per  plant 
is  used  to  calculate  the  proper  initial  values  for  simulations  started  from  a 
different  tuber  size  than  nominal  (Table  4). 

Tuber  death  rates  have  not  been  published.  For  use  in  POTAM,  the  value 
for  the  relative  tuber  death  rate  ( RDTU )  was  found  using  the  same  ap¬ 
proach  as  followed  for  VALLA.  A  RDTU  of  0.026  d1  (on  number  basis) 
was  found  for  the  SAV  in  a  shallow  ditch  in  Yerseke,  NL  (Van  Wijk  1988). 
Tuber  loss  from  grazing  is  expected  to  be  in  the  same  order  of  magnitude 
as  in  V.  americana. 

Sprouting  frequency  of  tubers  in  the  field  has  not  been  published.  Poten¬ 
tial  sprouting  frequency  is  high  in  a  temperate  climate,  being  <80  percent. 
It  is  increased  by  cold  stratification  (Van  Wijk  1989).  Sprouting  frequency 
decreases  with  tuber  burial  depth,  for  instance  (a)  small  tubers  in  the 
range  of  0.003  to  0.012  g  DW  tuber1  proved  less  likely  to  emerge  when 
planted  at  5  and  10  cm  than  at  2.5  cm  depth,  whereas  (b)  only  30  percent 
of  large  tubers  of  0.017  g  DW  sprouted  at  20  cm  depth  (Spencer  1987). 
Abundant  sprouting  was  observed  between  5.5  °C  and  10  °C  (Western  Ca¬ 
nal,  NL;  Best  and  Boyd  2003c)  and  25  °C  (Spencer  and  Anderson  1987, 

CA;  Van  Wijk  1989,  France),  but  sprouting  was  lower  at  25  °C  than  at 
20  °C  (Van  Wijk  1989).  Sprouting  occurs  early  in  the  growth  season,  i.e.,  in 
the  beginning  of  April  in  the  Western  Canal,  NL,  in  1987  (Best  and  Boyd 
2003c).  In  POTAM,  sprouting  is  governed  by  DVS  and,  thus,  indirectly  by 
temperature.  Minimum  weight  required  for  a  tuber  to  sprout  in  POTAM 
has  been  set  at  10  percent  of  the  initial  tuber  size. 

The  derivations  of  the  parameters  CRIFAC  and  SURPER  in  POTAM  are 
described  in  this  paragraph;  those  of  ROC  and  RCHST  are  described  under 
H.  verticillata.  The  elongation  potential  of  plantlets  emerging  from  tubers 
is  limited,  i.e.,  0.00714  m  (for  tubers  ranging  from  0.0054  to  0.155  g  DW 
tuber1  (Spencer  1987;  Spencer  and  Anderson  1987).  Thus,  plants  can  rise 
one  0.1-m  depth  layer  in  the  water  column  when  they  can  fill  that  layer 
with  a  minimum  of  0.0076  or  maximum  of  0.2170  g  DW.  Elongation  is 
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modeled  by  filling  each  successive  water  layer  from  the  sediment  to  a 
1.3-m  water  column  with  the  minimum  shoot  biomass  required  (0.0076  g 
plant  dry  weight  plant-1,  termed  CRIFAC).  A  tuber  survival  period 
( SURPER )  of  23  days  for  the  nominal  0.083-g-DW  tuber  size  was  calcu¬ 
lated,  following  the  same  approach  as  used  for  VALLA.  A  slightly  higher 
SURPER  of  27  days  was  selected  for  use  in  POTAM  since  this  value  al¬ 
lowed  the  simulated  SAV  to  survive  and  produce  biomass  similar  to  popu¬ 
lations  in  the  field,  while  with  a  SURPER  of  23  days  the  SAV  became  ex¬ 
tinct.  Larger  tubers  have  longer,  and  smaller  tubers  shorter  survival 
periods. 

H.  verticillata.  Tuber  bank  density  in  this  species  ranges  from  o  in  early 
summer  to  1,000  m-2  in  autumn  (Haller  et  al.  1976;  Bowes  et  al.  1979).  A 
tuber  bank  density  of  500  nr2  with  a  tuber  size  of  0.1  g  DW  found  in  well- 
established  tuber  banks,  in  Lake  Orange,  FL  (Bowes  et  al.  1979)  was  used 
for  calibration. 

Tubers  of  various  size  classes  have  been  found  in  the  field,  usually  around 
0.1  g  DW  tuber-1  (Van  et  al.  1978a).  Larger  tubers  have  been  found  also 
(Bowes  et  al.  1977).  The  largest  tubers  of  0.5  g  DW  tuber-1  were  reported  in 
plants  cultivated  in  the  laboratory  and  fertilized  (McFarland  and  Barko 
1990).  The  number  of  tubers  concurrently  formed  per  plant  was  derived 
from  tubers  collected  from  the  H.  verticillata  tuber  bank  in  Lake  Orange, 
FL  (Bowes  et  al.  1979),  and  ranged  from  7  to  11.  The  combination  of  7.0 
concurrently  initiated  tubers  ( NINTUB )  and  tuber  size  of  0.1  g  DW  tuber-1 
(INTUB)  was  selected  as  nominal.  A  typical  relationship  between  individ¬ 
ual  tuber  weight  and  concurrently  formed  tuber  number  per  plant,  which 
was  identified  in  both  other  tuber-forming  plants,  was  not  found  in  this 
plant  species. 

Tuber  death  rates  of  have  not  been  published.  For  use  in  HYDRIL,  the 
value  for  the  relative  tuber  death  rate  ( RDTU)  of  0.36  d-1  (on  number  ba¬ 
sis)  was  found  by  linear  extrapolation  from  a  density  of  500  m-2  (represen¬ 
tative  for  Lake  Orange,  FL)  to  o  nr2  in  3  years  (field  observation  North 
New  River  Canal,  FL;  D.L.  Sutton,  personal  communication  1995).  This 
value  is  similar  to  that  reported  for  a  potato  crop  in  Florida  (Ingram  and 
McCloud  1984).  Tuber  loss  from  grazing  is  expected  to  be  substantial  be¬ 
cause  of  the  high  tuber  bank  density,  and  relatively  high  compared  to 
other  tuber-forming  species. 
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Sprouting  frequency  of  tubers  in  the  field  has  not  been  published.  Poten¬ 
tial  sprouting  frequency  is  high  under  experimental  conditions,  i.e.,  up  to 
too  percent  (Bowes  et  al.  1977;  Spencer  and  Anderson  1986;  Steward  and 
Van  1986;  Sutton  and  Portier  1992).  Sprouting  may  be  triggered  by  a  cer¬ 
tain  red-far  red  ratio  in  the  light  reaching  the  plants’  photosensitive  pig¬ 
ments  (Spencer  and  Anderson  1986);  however,  the  latter  phenomenon  can 
only  play  a  role  when  wintering  shoot  biomass  is  present.  Phytochrome  is 
not  expected  to  trigger  tuber  sprouting  in  the  spring  for  H.  verticillata  in 
Lake  Orange,  since  it  winters  in  the  form  of  tubers  in  the  sediment  without 
shoot  biomass.  Sprouting  occurs  at  a  day  length  greater  than  13  h,  as  for 
potatoes  (Hahn  and  Hozyo  1983),  within  a  temperature  range  between 
18  °C  and  33  °C  (Haller  et  al.  1976),  and  may  start  early  in  the  growth  sea¬ 
son,  i.e.,  mid-March  in  Florida  lakes  (Haller  et  al.  1976).  Minimum  weight 
required  for  a  tuber  to  sprout  in  POTAM  has  been  set  at  10  percent  of  the 
initial  tuber  size. 

The  derivations  of  the  parameters  ROC  and  RCHST  in  HYDRIL  are  de¬ 
scribed  in  this  paragraph.  CRIFAC  and  SURPER  are  not  used  in  HYDRIL. 
Sprouting  tubers  deplete  tuber  carbohydrate  reserves.  Tuber  carbohy¬ 
drates  are  converted  into  plant  material  following  a  relative  tuber-to-plant 
conversion  rate  ( ROC)  of  0.0576  g  CH20  g  tuber  DW1  d1.  This  value  was 
derived  from  the  light  compensation  point  ( LCPt )  of  12  pE  nr2  s1  in  young, 
35-d  old,  H.  verticillata  shoots  (Bowes  et  al.  1977),  assuming  that  the 
originally  high  LCPo  of  too  decreases  exponentially  to  the  far  lower  LCPt 
due  to  conversion  of  tuber  into  plant  material,  according  to  the  following 
equation. 


LCPt  =  LCP0  x  e 


ROC(  35x24) 


(6) 


The  elongation  potential  of  plantlets  emerging  from  tubers  is  large,  i.e., 

12  m  g1  DW  on  average  (range  8.8  to  18.8  m  g-1  DW  (tuber  size  range  of 
0.011-0.056;  derived  from  Bowes  et  al.  1977)).  Elongation  is  modeled  by 
allowing  each  tuber  to  elongate  from  the  sediment  as  far  as  possible  or  up 
to  the  water  surface,  according  to  the  relation  coefficient  tuber  weight- 
stem  length  ( RCSHST)  of  12  m  g1  DW. 

M.  spicatum.  In  MILFO,  initial  rhizome/root  crown  biomass  has  been  set 
at  50  g  AFDW  nr2,  i.e.,  equal  to  the  lowest  shoot  biomass  found  in  Lake 
Wingra,  WI,  in  1970  (Adams  and  McCracken  1974),  and  equal  to  the 
belowground  biomass  measured  at  the  beginning  of  the  year  at  a  1.5-m 
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rooting  depth  in  the  same  lake  in  1977  (Smith  and  Adams  1986).  Rhizome/ 
root  crown  biomass  of  M.  spicatum  fluctuated  between  12  and  400  g 
DW  nr2  in  a  temperate  climate  (Madsen  1997),  and  between  32  and 
48  g  DW  nr2  in  a  tropical  climate  (Zutschi  and  Vass  1973). 

Sprouting  or  regrowth  of  the  rhizome/root  crown  in  the  field  has  not  been 
published.  Sprouting  potential  is  usually  high.  Sprouting  frequency  in  an 
established  community  is  probably  not  important,  unless  it  is  very  low,  as 
long  as  the  typical  plant  density  of  11  plants  nr2  is  somehow  reached  be¬ 
cause  plant  density  tends  to  play  a  lesser  role  in  biomass  production  com¬ 
pared  to  space  availability.  In  southern  areas  such  as  Texas  (lat  330  N), 
sprouting  occurs  already  in  March  (Madsen  1997).  In  northern  areas,  the 
timing  of  sprouting  may  be  similar,  but  no  observations  confirming  this 
have  been  published,  probably  because  at  that  time  water  temperatures 
are  still  very  low  impeding  field  work.  Minimum  weight  required  for  a  rhi¬ 
zome  to  sprout  ( CRRIZ)  is  12.0  g  dry  weight  nr2,  i.e.,  the  lowest  published 
value  (Madsen  1997). 

The  derivation  of  the  parameter  RDRIZ  in  MILFO  is  described  in  this 
paragraph;  those  of  ROC  and  RCHST  are  described  under  H.  verticillata. 
The  relative  rhizome  death  rate  was  not  derived  from  field  data  because 
the  data  availability  was  insufficient.  Instead,  a  relative  rhizome  death  rate 
has  been  derived  from  observations  on  terrestrial  rhizome  systems  where 
annual  turnover  rates  were  found  to  be  approximately  four  times  less  than 
those  of  aboveground  plant  biomass  in  the  growth  season,  but  dropping  by 
a  factor  of  1/100  in  inactive  periods  (Vogt  et  al.  1991).  Following  this  ap¬ 
proach,  a  tentative  relative  rhizome  death  rate  ( RDRIZ)  of  0.00042  d1  was 
calculated  (on  weight  basis),  being  1/100  of  the  M.  spicatum  plant  death 
rate.  Effects  of  grazing  by  waterfowl  on  the  milfoil  rhizome/root  crown 
systems  are  unknown,  but  they  are  expected  to  be  far  lower  than  on  H. 
verticillata  tubers. 

Light,  photosynthesis,  maintenance,  growth,  and  assimilate  partitioning 

V.  americana.  VALLA  is  calibrated  using  field  data  on  biomass  and  water 
transparency  (L)  from  Chenango  Lake,  NY,  1978  (Titus  and  Stephens 
1983)  and  climatological  data  from  Binghamton  (air  temperatures)  and 
Ithaca  (irradiance),  NY,  1978.  For  the  plant  community-specific  light  ex¬ 
tinction  coefficient  ( K)  in  VALLA,  the  value  of  0.0235  m2  g1  DW  was  used, 
measured  in  the  canopy  of  a  SAV  in  Lake  Mendota,  WI  (Titus  and  Adams 
1979a).  Other,  lower,  community-specific  extinction  coefficients  of  0.0051, 
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0.013-0.019,  and  0.018  m2  g1  DWhave  been  published  by  Blanch  et  al. 
(1998),  Titus  and  Adams  (1979a,  “vertical  K”)  and  Ikusima  (1970). 

The  reduction  in  light  interception  by  the  plant  through  epiphyte  shading 
in  VALLA  is  derived  from  measured  epiphyte  shading  up  to  a  0.43-fraction 
in  V.  americana  populations  in  the  Upper  Mississippi  River,  WI  (Best  et 
al.  2005;  Figure  6).  Light  interception  by  the  plant  is  multiplied  by  a  factor 
of  ( l-EPISHD ),  that  ranges  from  1  at  DVS  o  to  0.57  at  DVS  1.0  during  flow¬ 
ering,  to  1  at  DVS  2.3  in  a  temperate  climate  (and  DVS  >2.3  in  a  warmer 
climate).  This  fraction  may  be  typical  for  this  particular  plant  species  and 
site.  For  SAV  growing  in  aquatic  systems  with  less  epiphyte  shading  than 
in  the  Mississippi  River,  the  maximum  epiphyte  shading  may  have  to  be 
reduced  in  the  model  input  file. 


V.  americana 


Figure  6.  Relationship  between  development  phase  and  relative  epiphytic  light  interception 
used  for  calibration  of  VALLA. 

For  photosynthetic  activity  at  light  saturation  ( AMX)  in  VALLA,  the  value 
of  0.0165  g  C02  g1  DWtr1  was  used,  measured  in  V.  americana  in  the 
laboratory  (Titus  1977).  Light-  and  carbon-saturated  photosynthetic  rates 
can  be  far  higher  (Titus  and  Stone  1982;  Titus  et  al.  1990),  suggesting  that 
photosynthetic  activity  in  lakes  such  as  Lakes  Chenango  and  Mendota, 
where  DIC  concentrations  are  in  the  range  of  0.8  to  3.5  mmol  and  pH  7.6  - 
9.4,  can  be  carbon-limited.  For  photosynthetic  light-use  efficiency  ( EE )  in 
VALLA,  a  value  of  lixicr6  g  C02  J1  and  typical  for  C3  plants  was  used 
(Penning  de  Vries  and  Van  Laar  1982a).  Gross  assimilation  rate  at  light 
saturation  shows  a  distinct  seasonal  pattern  and  tends  to  decrease  with 
aging  (Titus  19 77).  Although  a  relative  function  describing  the  effect  of 
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aging  on  gross  assimilation  ( REDFT)  has  been  included  in  VALLA,  it  is  not 
active  in  the  nominal  version  (set  to  l)  because  in  simulations  this  effect 
was  negligible  compared  to  that  of  changes  in  irradiance.  Daily  changes  in 
pH  and  oxygen  concentrations  may  affect  AMX.  In  VALLA,  a  reduction 
factor,  REDAM,  can  be  used  to  take  effects  of  daily  changes  in  pH  and 
oxygen  concentration  into  account  by  reducing  the  AMX  by  a  factor  be¬ 
tween  o  and  l.  AMX  and  actual  photosynthesis  are  similar  at  pH  8.5  found 
in  V.  americana  communities  in  Chenango  Lake  (Titus  and  Stephens 
1983;  Titus  and  Stone  1982),  and,  therefore,  REDAM  currently  has  the 
value  of  1.  AMX  in  this  species  is  not  affected  by  changes  in  oxygen  con¬ 
centration.  In  VALLA,  the  effect  of  daytime  temperature  on  AMX  is  ac¬ 
counted  for  via  a  fitted,  relative  function,  AMTMPT,  which  is  based  on  the 
measured  photosynthetic  response  of  V.  americana  to  temperature  and 
has  its  optimum  at  35  °C  (Titus  and  Adams  1979a).  A  similar  relationship 
between  temperature  and  chlorophyll  concentration  was  found  by  Barko 
and  Filbin  (1982). 

The  relationship  between  photosynthesis  and  current  velocity  in  VALLA  is 
derived  from  measured  typical  plant  photosynthesis  in  the  related  Elodea 
nuttallii  at  low  current  velocity  (Best  and  Boyd  2003a),  and  typical  bio¬ 
mass  production  of  V.  americana  in  the  Red  Cedar  River,  WI,  at  high  cur¬ 
rent  velocity  (Best  et  al.  2005).  AMX  is  multiplied  by  the  relative  photo¬ 
synthetic  rate  (REDAM 1),  which  ranges  between  1  at  a  current  velocity  of 
7  cm  s1  and  o  at  a  current  velocity  of  81  cm  s_1  (Figure  7). 


V.  americana 


Figure  7.  Relationship  between  current  velocity  and  relative  photosynthetic  rate  used  for 
calibration  of  VALLA. 
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In  VALLA,  assimilates  are  partitioned  typically  over  the  plant  parts,  with 
71.8  percent  as  leaves,  15.9  percent  as  stems,  and  12.3  percent  as  roots  (Ti¬ 
tus  and  Stephens  1983;  Haller  1974).  Vertical  distribution  of  shoot  bio¬ 
mass  is  governed  by  a  dry  matter  partitioning  coefficient  function  ( DMPC ) 
following  the  typical  meadow  type  shape.  Distribution  proceeds  as  follows. 
First  plant  biomass  is  distributed  at  18.4  percent  over  the  three  depth  lay¬ 
ers  above  the  sediment  surface  and  at  11.4  percent  over  layers  4  and  5 
above  these.  Subsequently,  the  remaining  shoot  biomass  is  distributed 
over  maximally  seven  depth  layers  above  this  0.5  m  up  to  a  total  share  of 
9.7  percent.  When  five  depth  layers  or  less  are  present,  first  the  lowest  lay¬ 
ers  are  filled  according  to  the  typical  distribution  pattern,  and  subse¬ 
quently  the  remaining  shoot  biomass  is  added  by  equal  distribution  over 
the  found  layers.  Vertical  biomass  distribution  pattern  is  recalculated  and 
redistributed  by  VALLA  when  a  rooting  depth  other  than  nominal  (1.4  m) 
is  chosen  by  the  user.  In  the  model,  roots  always  contribute  12.3  percent 
to  total  plant  biomass. 

P.  pectinatus.  POTAM  is  calibrated  using  field  data  on  biomass  and  water 
transparency  (L)  from  the  Western  Canal  near  Zandvoort,  NL,  1987  (Best 
and  Boyd  2003c)  and  climatological  data  from  De  Bilt,  NL,  1987.  For  the 
plant  community-specific  light  extinction  coefficient  in  POTAM  ( K ),  the 
value  of  0.095  m2  g  DW1,  measured  in  the  canopy  of  a  SAV  in  the  same 
canal,  was  used  (Best  and  Boyd  2003c).  Other  community-specific  extinc¬ 
tion  coefficients  of  0.0183,  0.020,  and  1.79  m2  g  DW1  have  been  published 
by  Sher-Kaul  et  al.  (1995),  Westlake  (1965),  and  Best  et  al.  (2005). 

The  reduction  in  light  interception  by  the  plant  through  epiphyte  shading 
in  POTAM  is  derived  from  measured  epiphyte  shading  up  to  complete 
light  extinction,  i.e.,  a  l.oo-fraction,  in  P.  pectinatus  populations  in  the 
Upper  Mississippi  River,  WI  (Best  et  al.  2005;  Figure  8).  Light  intercep¬ 
tion  by  the  plant  is  multiplied  by  a  factor  of  ( l-EPISHD ),  that  ranges  from 
1  at  DVS  o  to  o  at  DVS  1.0  during  flowering,  to  1  at  DVS  2.3  in  a  temperate 
climate  (and  DVS  >2.3  in  a  warmer  climate).  This  fraction  may  be  typical 
for  this  particular  plant  species  and  site.  For  SAV  growing  in  aquatic  sys¬ 
tems  with  less  epiphyte  shading  than  the  Mississippi  River,  the  maximum 
epiphyte  shading  may  have  to  be  reduced  in  the  model  input  file. 
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P.  pectinatus 


Figure  8.  Relationship  between  development  phase  and  relative  epiphytic  light  interception 
used  for  calibration  of  POTAM,  HYDRIL,  and  MILFO. 

For  photosynthetic  activity  at  light  saturation  ( AMX )  in  POTAM,  the  value 
of  0.019  g  C02  g1  DW  h-1  was  used,  measured  in  P.  pectinatus  shoot  apices 
in  May  in  the  shallow,  fresh,  River  Susa,  Denmark  (Van  der  Bijl  et  al.  1989; 
AMX).  Other,  lower,  AMX  values  of  0.0099  and  0.012  g  C02  g  DW1  h1 
were  determined  under  laboratory  conditions  (Hootsmans  and  Vermaat 
1994;  Hootsmans  1994;  Madsen  and  Adams  1988).  Light-  and  carbon- 
saturated  photosynthetic  rates  of  P.  pectinatus  can  be  far  higher  (Best  and 
Boyd  2003c),  suggesting  that  photosynthetic  activity  in  waters  with  low 
DIC  concentrations  (<3.5  mmol  L1)  or  high  pH  (>8)  can  be  carbon- 
limited.  For  photosynthetic  light-use  efficiency  (EE)  in  POTAM,  the  same 
value  as  for  VALLA  was  used.  REDFT  and  REDAM  have  been  modeled  as 
in  VALLA.  In  POTAM,  the  effect  of  daytime  temperature  on  AMX  is  ac¬ 
counted  for  via  a  fitted,  relative  function,  AMTMPT,  which  is  based  on  the 
measured  photosynthetic  response  of  P.  pectinatus  to  temperature  and 
has  its  optimum  at  30.0  °C  (Best  and  Boyd  2003c). 

The  relationship  between  photosynthesis  and  current  velocity  in  POTAM 
is  derived  from  measured  typical  plant  photosynthesis  at  low  current  ve¬ 
locity  (Best  and  Boyd  2003c),  and  biomass  production  in  the  shallow  Ca¬ 
nadian  rivers  at  high  current  velocity  (Chambers  et  al.  1991).  AMX  is  mul¬ 
tiplied  by  the  relative  photosynthetic  rate  ( REDAMi ),  which  ranges 
between  1  at  a  current  velocity  of  7  cm  s1  and  o  at  a  current  velocity  of 
94  cm  s_1  (Figure  9). 
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P.  pectin  at  us 


Figure  9.  Relationship  between  current  velocity  and  relative  photosynthetic  rate  used  for 
calibration  of  POTAM,  HYDRIL,  and  MILFO. 

In  POTAM,  assimilates  are  partitioned  typically  over  the  plant  parts,  with 
73.1  percent  as  leaves,  18.3  percent  as  stems  and  8.6  percent  as  roots  (Best 
and  Boyd  2003c;  Sher-Kaul  et  al.  1995).  Vertical  distribution  of  shoot  bio¬ 
mass  is  governed  by  a  dry  matter  partitioning  coefficient  function  ( DMPC) 
following  the  typical  canopy  type  shape.  Distribution  proceeds  as  follows. 
First  plant  biomass  is  distributed  for  4.4  percent  over  the  upper  depth 
layer  1,  for  4.3  percent  over  layer  2,  for  23.1  percent  over  layer  3,  for 
25.4  percent  over  layer  4,  and  for  21.3  percent  over  layer  5  (78.4  percent  of 
total  plant  biomass;  Best  and  Boyd  2003c).  Subsequently,  the  remaining 
shoot  biomass  is  distributed  equally  over  the  depth  layers  below  this  0.5  m 
up  to  a  total  share  of  13  percent.  When  the  biomass  in  a  layer  is  lower  than 
the  critical  biomass,  the  total  number  of  layers  is  reduced  by  one  and  the 
distribution  procedure  is  redone.  When  five  depth  layers  or  less  are  pre¬ 
sent,  first  the  upper  layers  are  filled  according  to  the  typical  distribution 
pattern,  and  subsequently  the  remaining  shoot  biomass  is  added  by  equal 
distribution  over  the  found  layers.  Vertical  biomass  distribution  pattern  is 
recalculated  and  redistributed  by  POTAM  when  a  rooting  depth  other  than 
nominal  (1.3  m)  is  chosen  by  the  user.  In  POTAM,  roots  always  contribute 
8.6  percent  to  total  plant  biomass. 

H.  verticillatCL.  HYDRIL  is  calibrated  using  field  data  on  biomass  and  wa¬ 
ter  transparency  (L)  from  Lake  Orange,  FL,  1977  (Bowes  et  al.  1979)  and 
climatological  data  from  Gainesville,  FL,  1980.  For  the  plant  community- 
specific  light  extinction  coefficient  ( K)  in  HYDRIL,  the  value  of  0.01  m2  g1 
DW,  measured  in  the  canopy  of  a  H.  verticillata  community  in  Lake 
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Trafford,  FL,  was  used  (Bowes  et  al.  1979).  The  same  value  was  reported 
for  a  H.  verticillata  Caspary  community  in  Lake  Biwa,  Japan  (Ikusima 
1970). 

The  reduction  in  light  interception  by  the  plant  through  epiphyte  shading, 
and  the  relationship  between  photosynthesis  and  current  velocity  have 
been  modeled  as  in  POTAM  based  on  calibration  data  pertaining  to 
VALLA  because  calibration  data  on  H.  verticillata  were  lacking.  Once  per¬ 
tinent  calibration  data  are  available,  these  can  be  inserted  in  the  model  in¬ 
put  file. 

For  photosynthetic  activity  at  light  saturation  ( AMX)  in  HYDRIL,  the 
value  of  0.0158  g  C02  g1  DW  h1  was  used,  measured  in  shoot  apices  in 
Lake  Orange,  FL,  in  May  (Bowes  et  al.  1979;  Van  et  al.  1978a).  Light-  and 
carbon-saturated  photosynthetic  rates  can  be  far  higher  (Van  et  al.  1976), 
suggesting  that  photosynthetic  activity  in  waters  with  low  DIC  concentra¬ 
tions  (<3.5  mmol  L1)  or  high  pH  (>8)  can  be  carbon  limited.  For  photo¬ 
synthetic  light-use  efficiency  {EE)  in  HYDRIL,  the  same  value  as  for 
VALLA  was  used. 

REDFT has  been  modeled  as  in  VALLA.  In  HYDRIL,  a  reduction  factor, 
REDAM,  can  be  used  to  take  effects  of  daily  changes  in  pH  and  oxygen 
concentration  into  account  by  reducing  the  AMX  by  a  factor  between  o 
and  1  for  the  whole  day.  Daily  changes  in  pH  and  oxygen  concentrations 
affect  AMX  in  H.  verticillata  significantly.  REDAM  currently  has  the  value 
of  0.581  because  changes  in  pH  in  the  SAV  in  Florida  lakes  inhibited  pho¬ 
tosynthesis  by  a  factor  of  0.83  and  changes  in  02  concentration  by  a  factor 
of  0.7  (Van  et  al.  1976;  Bowes  et  al.  1979).  In  HYDRIL,  the  effect  of  day¬ 
time  temperature  on  AMX  is  accounted  for  via  a  fitted,  relative  function, 
AMTMPT,  which  is  based  on  the  measured  photosynthetic  response  of  H. 
verticillata  to  temperature  and  has  its  optimum  at  24.0  °C  (Barko  and 
Smart  1981).  A  higher  optimum  of  36.5  °C,  found  by  W.T.  Haller,  IFAS, 
Gainesville,  FL  (unpublished  results  1984)  and  by  Van  et  al.  (1976),  gener¬ 
ated  unrealistically  high  biomass  in  simulations,  and  was,  therefore,  not 
used  for  model  calibration. 

In  HYDRIL,  assimilates  are  partitioned  typically  over  the  plant  parts, 
with  34  percent  as  leaves,  60  percent  as  stems,  and  6  percent  as  roots 
(Haller  and  Sutton  1975;  Van  et  al.  1978b;  Van  der  Zweerde  1981).  Vertical 
distribution  of  shoot  biomass  is  governed  by  a  dry  matter  partitioning 
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coefficient  function  (DMPC)  following  the  typical  canopy  type  shape.  Dis¬ 
tribution  proceeds  as  follows.  First,  plant  biomass  is  distributed  for  21  per¬ 
cent  over  the  two  upper  depth  layers  1  and  2,  for  10  percent  over  the  two 
successive  layers  3  and  4,  and  for  9  percent  over  the  two  successive  layers 
5  and  6  (80  percent  of  total  plant  biomass;  Ambasht  and  Ram  1976;  Bowes 
et  al.  1979).  Subsequently,  the  remaining  shoot  biomass  is  allocated 
equally  to  the  depth  layers  below  this  0.6  m  up  to  a  total  share  of  14  per¬ 
cent.  Vertical  biomass  distribution  pattern  is  recalculated  and  redistrib¬ 
uted  by  HYDRIL  when  a  rooting  depth  other  than  nominal  (1.0  m)  is  cho¬ 
sen  by  the  user.  In  HYDRIL,  roots  always  contribute  6  percent  to  total 
plant  biomass. 

M.  spicatum.  MILFO  is  calibrated  as  follows:  (1)  For  a  SAV  in  a  temperate 
climate,  using  field  data  on  biomass  (Adams  and  McCracken  1974),  water 
transparency  (L),  temperature  and  irradiance  from  Lake  Wingra,  WI,  1970 
(Lee  and  Kluesener  1972);  (2)  For  a  SAV  in  a  tropical  climate,  using  field 
data  on  biomass  from  Kashmir  lakes,  India  (Zutschi  and  Vass  1973)  and 
climatological  data  from  Patancheru,  India,  1978.  For  the  plant  commu¬ 
nity-specific  light  extinction  coefficient  (/<")  in  MILFO,  the  value  of  0.006 
m2  g'1  DW,  measured  in  the  canopy  of  a  SAV  in  Lake  Wingra,  WI,  was  used 
(Titus  and  Adams  1979a).  Another,  higher,  extinction  coefficient  of  0.01 
m2  g1  DW  has  been  published  by  Ikusima  (1970). 

The  reduction  in  light  interception  by  the  plant  through  epiphyte  shading, 
and  the  relationship  between  photosynthesis  and  current  velocity  have 
been  modeled  as  in  POTAM  based  on  calibration  data  pertaining  to 
POTAM  because  calibration  data  on  M.  spicatum  were  lacking.  Once  per¬ 
tinent  calibration  data  are  available,  these  can  be  inserted  in  the  model  in¬ 
put  file. 

For  photosynthetic  activity  at  light  saturation  ( AMX)  in  MILFO,  the  value 
of  0.0165  g  C02  g1  DW  h1  was  used,  measured  in  M.  spicatum  in  Lake 
Wingra,  WI,  in  May  (Adams  and  McCracken  1974;  pH  8  and  a  total  alka¬ 
linity  of  190  mg  L1;  AMX).  For  photosynthetic  light-use  efficiency  (EE)  in 
POTAM,  the  same  value  as  for  VALLA  was  used.  REDFT  has  been  mod¬ 
eled  as  in  VALLA.  Daily  changes  in  pH  and  oxygen  concentrations  may  af¬ 
fect  AMX.  A  reduction  factor,  REDAM,  can  be  used  to  take  effects  of  daily 
changes  in  pH  and  oxygen  concentration  into  account  by  reducing  the 
AMX  by  a  factor  between  o  and  1.  In  MILFO,  actual  photosynthesis  is  re¬ 
duced  by  50  percent  compared  to  AMX  at  pH  8.8  found  in  M.  spicatum 
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SAV  in  Lake  Wingra  (Titus  and  Stone  1982;  Adams  and  McCracken  1974), 
and,  therefore,  REDAM  currently  has  the  value  of  0.5.  AMX  in  this  species 
is  negligibly  affected  by  changes  in  oxygen  concentration  (Van  et  al.  1976). 
In  MILFO,  the  effect  of  daytime  temperature  on  AMX  is  accounted  for  via 
a  fitted,  relative  function,  AMTMPT,  which  is  based  on  the  measured  pho¬ 
tosynthetic  response  of  M.  spicatum  to  temperature  and  has  its  optimum 
at  35.0  °C  (Titus  and  Adams  1979a;  confirmed  by  Stanley  and  Nailor  1972). 

In  MILFO,  assimilates  are  partitioned  typically  over  the  plant  parts,  with 
47  percent  as  leaves,  47  percent  as  stems  and  6  percent  as  roots  (Adams  et 
al.  1974).  Vertical  distribution  of  shoot  biomass  is  governed  by  a  dry  mat¬ 
ter  partitioning  coefficient  function  ( DMPC)  following  the  typical  canopy 
type  shape.  Distribution  proceeds  as  follows.  First,  plant  biomass  is  dis¬ 
tributed  for  10  percent  over  the  upper  depth  layer  1,  for  16  percent  over 
layer  2,  for  17  percent  over  layer  3,  for  10  percent  over  layer  4,  and  for  8 
percent  over  layer  5  (61  percent  of  total  plant  biomass;  Adams  et  al.  1974). 
Shoot  biomass  distribution  over  the  lower  water  layers  is  equal  up  to  a  to¬ 
tal  biomass  share  of  5  percent.  The  remaining  shoot  biomass  of  34  percent 
is  distributed  equally  over  all  layers.  Vertical  biomass  distribution  pattern 
is  recalculated  and  redistributed  by  MILFO  when  a  rooting  depth  other 
than  nominal  (1.5  m)  is  chosen  by  the  user.  In  MILFO,  roots  always  con¬ 
tribute  6  percent  to  total  plant  biomass. 

Flowering,  translocation,  senescence,  and  formation  of  wintering  organs 
V.  americana.  Flowering  starts  late  in  June  (Titus  and  Stephens  1983). 

Translocation  has  not  been  measured  directly  in  SAV.  In  V.  americana 
translocation  products  end  up  largely  in  tubers,  and,  thus,  in  VALLA 
translocation  rate  is  based  on  relative  tuber  growth  rate  (see  below).  A 
translocation  rate  of  24.7  percent  of  net  photosynthesis  generated  realistic 
biomass  production  in  simulations  and  was  used  for  calibration.  The  total 
transport  “cost”  of  downward  translocation  was  estimated  at  a  factor  of 
1.05  (10/9  ><36/38;  CVT)  because  transport  of  glucose  costs  dry  matter 
(decrease  by  a  factor  of  36/38),  whereas  conversion  of  starch  into  glucose 
increases  it  (by  a  factor  of  10/9).  Translocation  estimates  in  other  plant 
species  range  from  19  percent  of  net  production  in  seagrasses  (Wetzel  and 
Neckles  1996),  29  percent  in  cassava  (Gijzen  1985)  to  35  percent  in  certain 
potato  varieties  (Kooman  1995). 
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Senescence  rates  of  SAV  have  not  been  published.  The  value  for  the  rela¬ 
tive  death  rate  in  VALLA  was  determined  using  the  same  approach  as  fol¬ 
lowed  for  the  tuber  death  rate.  Thus,  a  relative  death  rate  (RDRT,  RDST) 
of  0.021  d1  was  calculated  for  the  SAV  in  Chenango  Lake,  NY  (Titus  and 
Stephens  1983). 

In  VALLA,  tubers  are  formed  after  flowering,  under  short  day  conditions 
(<14.7  h)  and  within  a  temperature  range  between  5  °C  and  25  °C,  based  on 
field  data  from  Titus  and  Stephens  (1983)  and  Donnermeyer  (1982).  The 
intensity  of  translocation  is  governed  by  the  actual  relative  tuber  growth 
rate,  RTRL,  consuming  24.7  percent  of  net  production  and  corrected  for 
the  total  transport  “cost”  CVT  and  temperature.  In  VALLA,  the  maximum 
relative  tuber  growth  rate  (RTR)  was  determined  by  applying  the  same  dif¬ 
ferential  equation  as  commonly  used  for  simple  exponential  growth 
(Thornley  and  Johnson  1990b;  Hunt  1982)  to  tuber  data  collected  in  the 
field  (Titus  and  Stephens  1983).  Thus,  a  RTR  of  0.247  d1  at  a  reference 
temperature  of  20  °C  was  computed.  Tubers  grow  until  the  preselected 
critical  weight  of  the  new  tuber  class  (TWCTUB)  is  reached.  The  value  of 
TWCTUB  is  calculated  by  multiplication  of  the  individual  tuber  weight 
with  the  concomitantly  initiated  tuber  number  and  plant  density  (5.5  x 

O. 09  x  30). 

P.  pectinatus.  Flowering  starts  in  mid-July  (Best  and  Boyd  2003c). 

Translocation  rate  in  POTAM  could  not  be  based  on  tuber  growth  rate  be¬ 
cause  the  availability  of  tuber  growth  data  was  insufficient.  A  translocation 
rate  of  19  percent  of  net  photosynthesis,  typical  for  the  related  seagrasses 
(Wetzel  and  Neckles  1996),  generated  realistic  biomass  production  in 
simulations  and  was  used  for  calibration.  The  total  transport  “cost”  of 
downward  translocation  was  estimated  at  a  factor  of  1.05  (CVT). 

The  value  for  the  relative  death  rate  in  POTAM  was  determined  using  the 
same  approach  as  followed  for  the  tuber  death  rate.  Thus,  a  relative  death 
rate  (RDRT,  RDST)  of  0.047  d1  was  calculated  for  the  SAV  in  the  Western 
Canal  near  Zandvoort,  NL  (Best  and  Boyd  2003c). 

In  POTAM,  tubers  are  formed  after  flowering  under  short-day  conditions 
(<16  h)  and  within  a  temperature  range  between  5  °C  and  25  °C  (Spencer 
and  Anderson  1987;  Van  Wijk  1988).  In  POTAM,  the  maximum  relative 
tuber  growth  rate  (RTR)  at  was  set  equal  to  the  translocation  rate  at 
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0.19  d1  at  a  reference  temperature  of  20  °C.  Tubers  grow  until  the  pre¬ 
selected  critical  weight  of  the  new  tuber  class  (TWCTUB)  is  reached. 
TWCTUB  is  calculated  by  multiplication  of  the  individual  tuber  weight 
with  the  concomitantly  initiated  tuber  number,  and  plant  density  (8  x 
0.083  x  30). 

H.  verticillata.  Flowering  starts  at  the  end  of  May  (Haller  et  al.  1976). 

Translocation  rate  in  HYDRIL  could  not  be  based  on  tuber  growth  rate  be¬ 
cause  the  availability  of  tuber  growth  data  was  insufficient.  A  translocation 
rate  of  40  percent  of  net  production,  slightly  higher  than  measured  in  po¬ 
tato  (Kooman  1995),  generated  realistic  biomass  production  in  HYDRIL 
and  was  used  for  calibration.  In  this  case,  only  the  dry  matter  cost  of  glu¬ 
cose  transport  of  glucose  was  included  in  the  model,  i.e.,  a  factor  of  1.1 
(CVT). 

The  value  for  the  relative  death  rate  in  HYDRIL  was  determined  using  the 
same  approach  as  followed  for  the  tuber  death  rate.  A  relative  death  rate 
(RDRT,  RDST)  of  0.033  d1  was  calculated  for  the  SAV  in  Lake  Orange,  FL 
(Bowes  et  al.  1979). 

In  HYDRIL,  tubers  are  formed  after  flowering  under  short-day  conditions 
(>10  h  and  <16  h)  and  within  a  temperature  range  between  14  °C  and  33  °C 
(Van  et  al.  1978b).  In  HYDRIL,  the  maximum  relative  tuber  growth  rate 
(RTR)  at  was  set  equal  to  the  translocation  rate  at  0.40  d1  at  a  reference 
temperature  of  20  °C.  Tubers  grow  until  the  preselected  critical  weight  of 
the  new  tuber  class  (TWCTUB)  is  reached.  TWCTUB  is  calculated  by  mul¬ 
tiplication  of  the  individual  tuber  weight  with  the  concomitantly  initiated 
tuber  number  and  plant  density  (7  x  0.1  x  35). 

M.  spicatum.  Flowering  starts  in  the  first  half  of  June  in  temperate  condi¬ 
tions  (Adams  and  McCracken  1974),  and  may  last  the  whole  growth  season 
in  the  tropics  (Zutschi  and  Vass  1973). 

Translocation  rate  (TRAFAC)  in  MILFO  could  not  be  based  on  rhizome/ 
root  crown  growth  rate  because  the  availability  of  growth  data  was  insuffi¬ 
cient.  A  translocation  rate  of  35  percent  of  net  production,  equal  to  that 
measured  in  potato  (Kooman  1995),  generated  realistic  biomass  produc¬ 
tion  in  MILFO  and  was  used  for  calibration.  Costs  of  glucose  transport  and 
conversion  of  starch  into  glucose  (CVT)  have  been  modeled  as  in  VALLA. 
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In  MILFO,  the  value  for  the  relative  death  rate  was  determined  using  the 
same  approach  as  followed  for  the  death  rates  of  other  modeled  SAV.  A 
relative  death  rate  (RDRT,  RDST)  of  0.042  d1  was  determined  for  the  SAV 
in  Lake  Wingra,  WI  (Adams  and  McCracken  1974). 
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4  Simulations 

VALLA 

VALLA  calibration:  Simulated  and  measured  behavior  of  a  V.  americana 
community  in  Chenango  Lake,  NY 

The  seasonal  changes  in  biomass  of  plants  and  tubers,  as  simulated  by 
VALLA,  are  shown  in  Figure  to.  Simulated  plant  biomass  compared  well 
with  average  plant  biomass  measured  in  Chenango  Lake,  NY  (Titus  and 
Stephens  1983).  Plant  biomass  reached  its  maximum  at  the  same  time, 
and  peak  biomass  was  somewhat  higher  in  the  simulated  than  in  the 
measured  SAV,  notably  56.1  versus  50.1  g  DW  nr2.  The  latter  may  be  due 
to  the  relatively  large  tuber  size  /concurrently  initiated  tuber  number 
combination  (0.09  g  DW  tuber1,  5.5  tubers  plant'1)  used  to  initiate  this 
nominal  run.  Measured  tuber  size  was  0.055  g  DW  tuber1,  and  another 
model  run  starting  with  the  measured  tuber  size  generated  a  peak  biomass 
of  48  g  DW  nr2  (see  also  Figure  12).  Simulated  transport  of  carbohydrates 
was  substantial  in  the  beginning  of  the  growth  season  when  upward  car¬ 
bohydrate  remobilization  from  the  tubers  supports  initial  sprouting,  but 
far  higher  after  flowering  when  downward  carbohydrate  translocation 
from  plant  organs  supports  the  filling  of  the  tubers.  Carbohydrate  trans¬ 
port  could  be  in  the  same  range  as  net  assimilation  at  the  beginning  and 
end  of  the  growth  season.  Maintenance  respiration  was  usually  considera¬ 
bly  lower  than  assimilation  but  could  be  in  the  same  range  of  translocation 
just  after  flowering  (data  not  shown). 

Running  VALLA  with  (24-h  average)  air  temperatures  with  a  lag -period  of 
1  day  instead  of  with  measured  water  temperatures  as  forcing  variables 
yielded  higher  assimilation  (Figure  11)  and  plant  biomass  values  than 
found  in  the  nominal  case.  This  difference  was  attributed  to  the  fact  that 
water  temperatures  in  the  lake  were  relatively  low  compared  to  air  tem¬ 
peratures  because  of  the  large  inflow  of  groundwater  (Titus  and  Stephens 
1983).  In  the  authors’  experience,  a  lag  period  of  7  days  between  model 
daily  air  and  measured  temperatures  usually  describes  this  relationship 
well  for  shallow,  up  to  5-  or  6-m  deep  water  bodies  with  a  relatively  small 
inflow  of  groundwater. 


Biomass  (g  DW  m 
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Figure  10.  Simulated  biomass  of  plants  (A),  dormant  and  new  tuber  numbers  (B),  and  measured  plant 

biomass  (C)  of  V.  americana  in  Chenango  Lake,  NY.  Nominal  run.  Field  data  1987  from  Titus  and 
Stephens  (1983);  climatological  data  1987,  Binghamton,  NY  (lat  42°  15'  N,  long  75°  50'  W);  water 
depth  1.4  m;  light  extinction  coefficient  0.43  nr1. 


Figure  11.  Simulated  photosynthetic  rates  of  V.  americana  in  Chenango  Lake,  NY,  with  water 
or  air  temperatures  as  input  (initial  plant  parameter  values  as  in  nominal  run). 
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V.  americana 


Time  (d)  Time  (d)  Time  (d) 

Figure  12.  Simulated  biomass  of  plants  and  tubers  of  V.  americana  in  Chenango  Lake,  NY,  started  from 

different  initial  biomass,  but  run  in  the  same  environmental  and  climatological,  nominal,  conditions. 
(A)  Plant  biomass  0;  tuber  size  0.09  g  DW;  tuber  bank  233  nr2;  (B)  Plant  biomass  50  g  DW  nr2; 
tuber  size  0.09  g  DW;  tuber  bank  density  233  rrr2;  (C)  Plant  biomass  0;  tuber  size  0.055  g  DW; 
tuber  bank  density  233  nr2. 

Running  VALLA  for  the  same  lake  and  year,  but  with  both  plants  and  tu¬ 
bers  initially  present,  showed  that  peak  plant  biomass  was  then  far  higher 
and  more  tuber  classes  were  finished  (4  instead  of  3  in  the  nominal  case; 

Figure  12).  This  large  difference  in  peak  biomass  is  due  to  the  ability  of  the 
plant  community  to  fully  capture  the  high  spring  irradiance  at  this  latitude 
of  430  N,  which  they  can  not  without  wintering  shoots.  Thus,  wintering 
shoots  would  provide  a  distinct  advantage  for  this  plant  species.  However, 
wintering  shoots  have  only  been  reported  to  occur  in  a  tropical  climate.  A 
simulation  started  from  the  measured  tuber  size/chosen  concurrently  ini¬ 
tiated  tuber  number  combination  yielded  peak  plant  biomass  that  was  al¬ 
most  equal  in  simulated  and  measured  plant  community  (Figure  12C;  Fig¬ 
ure  10),  and  the  simulated  tuber  numbers  were  within  the  range  found  in  a 
V.  americana  community  in  the  same  lake  (Titus  and  Stephens  1983). 

VALLA  validation:  Simulated  and  measured  behavior  of  a  V.  americana 
community  at  other  latitudes 

To  explore  whether  VALLA  could  be  used  to  simulate  behavior  of  a  V. 
americana  community  at  other  sites  besides  the  nominal  one,  runs  were 
made  for  one  more  western  site,  Lake  Mendota,  WI,  and  another,  tropical 
more  southern  site,  ponds  at  Fort  Lauderdale,  FL. 


ERDC/EL  TR-07-14 


62 


A  simulation  was  performed  of  SAV  in  Lake  Mendota,  WI,  starting  from 
site-specific  community,  water  depth,  transparency,  and  climatological 
data.  Plant  community-specific  data  included  the  following:  initial  plant 
biomass  absent,  tuber  size  0.03  g  DW,  concurrently  initiated  tuber  num¬ 
ber  1.5  plant1,  and  tuber  bank  density  233  nr2.  Site-specific  environmental 
data  included  the  following:  water  depth  1.2  m  and  light  extinction  coeffi¬ 
cient  0.7  nr1.  In  these  conditions  simulated  plant  biomass  remained  low, 
maximally  25  g  DW  nr2,  and  no  tuber  classes  were  finished.  Because  the 
simulated  maximum  plant  biomass  was  far  lower  than  published  (average 
peak  biomass  344,  and  minimum  266  g  DW  nr2;  Titus  and  Adams  1979a), 
and  V.  americana  populations  in  Lake  Mendota  have  been  described  as 
stable,  several  other  simulations  were  done  to  explore  community-  and 
site-specific  characteristics  favoring  such  a  sustainable  population.  In¬ 
creasing  the  tuber  size  from  0.03  to  0.09  g  DW  increased  peak  biomass 
from  25  to  105  g  DW  nr2  and  the  number  of  finished  tuber  classes  to  5.  In¬ 
creasing  water  transparency  by  decreasing  the  light  extinction  coefficient 
from  0.7  to  0.4  nr1  increased  peak  biomass  from  25  to  130  g  DW  nr2  and 
the  number  of  finished  tuber  classes  to  7.  Introduction  of  the  largest  tuber 
size  published,  0.11  g  DW,  increased  peak  biomass  to  150  g  DW  nr2  but 
allowed  only  5  tuber  classes  to  be  finished.  Even  in  very  clear  water  (ex¬ 
tinction  coefficient  0.4  nr1),  a  very  shallow,  0.2-m  water  depth,  and  with 
the  most  profitable  tuber  size,  the  SAV  would  only  be  able  to  produce 
maximally  195  g  DW  nr2  in  Lake  Mendota  (Figure  13B).  The  latter  result 
led  the  authors  to  believe  that  the  high  plant  biomass  range  of  266  to  344  g 
DW  nr2  published  (Titus  and  Adams  1979a)  is  either  an  overestimate,  or 
that  spring  growth  starts  partly  from  wintering  shoots. 

A  simulation  was  performed  of  a  V.  americana  community  in  earthen 
ponds  near  Fort  Lauderdale,  FL,  starting  from  site-specific  community, 
water  depth,  transparency,  and  climatological  data.  Plant  community- 
specific  data  included  the  following:  initial  plant  biomass  50  g  DW  nr2,  tu¬ 
ber  size  0.09  g  DW,  concurrently  initiated  tuber  number  5.5  plant’1,  tuber 
bank  density  233  nr2.  Site-specific  environmental  data  included  the  fol¬ 
lowing:  water  depth  1.5  m  and  light  extinction  coefficient  0.4  nr1.  In  these 
conditions,  simulated  plant  biomass  was  high,  maximally  403  g  DW  nr2, 
and  three  tuber  classes  were  finished  at  the  very  end  of  the  year.  Tuber 
weights  and  numbers  of  the  Fort  Lauderdale  community  were  not  pub¬ 
lished,  so  comparison  between  simulated  and  measured  tuber  data  is  not 
possible.  Simulated  maximum  plant  biomass  in  this  case  was  within  the 
measured  peak  biomass  range  of  298  to  496  g  DW  nr2  (Figure  13C). 
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Maximum  plant  biomass  was  enhanced  to  450  g  DW  nr2  by  decreasing  the 
K  value  from  0.0235  to  0.019  m2  g  DW'1,  and  even  surpassed  the  maxi¬ 
mum  measured  value  by  decreasing  the  K  value  further  to  0.0051  m2  g 
DW1. 


Comparison  of  simulated  biomass  produced  in  the  various  climatological 
conditions  (Figure  13)  indicates  that  in  a  temperate  climate  generally  less 
biomass  is  produced  and  investment  in  vegetative  reproduction  is  rela¬ 
tively  higher  than  in  a  tropical  climate. 


V.  americana 


Time  (d)  Time  (d)  Time  (d) 


Figure  13.  Simulated  biomass  of  plants  and  tubers  of  V.  americana  at  sites  differing  in  latitude.  (A)  Chenango 
Lake,  NY  (lat  42°  15'  N,  long  75°  50'  W);  tuber  size  0.09  g  DW,  tuber  bank  density  233  rrr2;  water 
depth  1.4  m;  light  extinction  coefficient  0.43  nr1;  climatological  data  1987;  measured  data  1978, 
Titus  and  Stephens  (1983).  (B)  Lake  Mendota,  Wl  (lat  43°  08'  N,  long  89°  20'  W);  tuber  size  0.11  g 
DW,  tuber  bank  density  233  nr2;  water  depth  0.2  m;  light  extinction  coefficient  0.4  nr1;  climato¬ 
logical  data  1971-75;  measured  data  1976,  Titus  and  Adams  (1979a).  (C)  Fort  Lauderdale  ponds, 
FL  (lat  27°  58'  N,  long  82°  32'  W);  plant  biomass  50  g  DW  nr2;  tuber  size  0.09  g  DW,  tuber  bank 
density  233  rrr2;  K  value  0.018  m2  g1  DW;  water  depth  1.5  m;  light  extinction  coefficient  0.4  nr1; 
climatological  data  1975-84;  measured  data  1973-74,  Haller  (1974). 


VALLA  example  applications 

Application  1:  Historical  and  simulated  behavior  of  a  V.  americana  community 
in  a  riverine  environment  subject  to  flooding 

V.  americana  is  an  important  riverine  macrophyte  that  provides  food  and 
habitat  resources  for  waterfowl,  fishes,  and  invertebrates  in  the  Upper 
Mississippi  River  (UMR),  WI.  Navigation  pools  along  the  Mississippi  Fly¬ 
way  have  historically  been  used  by  migrating  waterfowl  as  staging  areas  in 
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part  because  of  abundant  populations  of  V.  americana  (Bellrose  et  al. 

1983;  Korschgen  and  Green  1988).  Recent  declines  in  these  populations 
have  caused  concern  and  have  been  attributed  to  eutrophication,  competi¬ 
tion  by  other  macrophytes  such  as  M.  spicatum  L.  and  Nelumbo  lutea 
Willd.,  drought,  and  flooding.  It  is  desirable  to  improve  management  pro¬ 
grams  aimed  at  enhancing  V.  americana  populations,  but  these  programs 
require  an  improved  understanding  of  the  population  dynamics  of  the  lat¬ 
ter  and  factors  affecting  these.  Detailed  biomass  dynamics  of  V.  ameri¬ 
cana  have  been  recorded  in  the  early  1980s,  when  populations  were  still 
substantial  (Pool  9;  Donnermeyer  1982).  Subsequently,  populations  have 
been  included  in  regular  surveys  of  several  navigation  pools,  but  only  by 
species  presence  or  absence,  not  characterized  by  plant  and  tuber  biomass. 

Simulations  using  VALLA  were  carried  out  to  evaluate  the  effects  of  daily 
changes  in  water  level  during  different  hydrological  years  on  a  typical  V. 
americana  in  Pool  8  of  the  UMR.  Stage  data  collected  at  the  dam  of  Pool  8 
were  used  to  document  water  level  fluctuations  over  a  10-year  period,  i.e. 
from  1985  to  1994.  In  this  period,  1985  is  considered  as  a  normal  hydro- 
logical  year,  1986  as  a  normal  flood  year  with  floods  in  spring  and  autumn, 
1993  as  an  abnormal  flood  year  with  one  flood  in  summer,  and  1986  as  a 
drought  year  (Figure  14).  Disappearance  of  large  portions  of  V.  americana 
populations  in  1993  and  1988  have  been  reported  (Rogers  1994;  Spink  and 
Rogers  1996).  The  simulations  were  done  starting  from  a  nominal  SAV, 
site-specific  environmental  data,  and  site-  and  year-specific  climatological 
data.  Site-specific  environmental  data  included  the  following:  water  depth 
daily  varying  as  would  be  experienced  by  a  community  at  0.5-m  rooting 
depth;  light  extinction  coefficient  ranging  from  2.619  to  3.173  nr1  during 
the  period  of  May  to  October,  and  set  to  2.0  nr1  the  rest  of  the  year  (con¬ 
verted  via  Giesens  relationship  from  Secchi  disk  readings,  cf  Giesen  et  al. 
(1990),  correlated  with  10-year  data  on  total  suspended  solids  concentra¬ 
tions).  Historical  data  indicate  that  in  1980  plant  biomass  peaked  in  mid- 
August  at  217  g  DW  nr2,  and  that  tuber  biomass  was  low  (maximally  14-16 
g  DW  nr2;  Donnermeyer  1982;  Figure  15A). 

Running  VALLA  with  nominal  initial  plant  biomass  and  tuber  bank  data  at 
a  constant  1.1-m  water  depth  (annual  average  in  the  historical  dataset) 
with  a  10-year  average  climate,  indicated  that  a  peak  biomass  of  only  60  g 
DW  nr2  was  formed,  but  no  tubers.  Water  fluctuations  were  substantial  in 
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Figure  14.  Water  level  fluctuations  over  a  10-year  period  measured  at  the  dam  of  Pool  8  of  the  Upper  Mississippi  River,  Wl  (long  91°  30'  E,  lat  43°  10'  N); 
data  J.  H.  Wlosinski,  La  Crosse,  Wl.  Flat  pool  is  considered  to  be  the  average  summer  (June,  July,  August)  value  during  normal  hydrological 
years  over  the  1985-94  period;  abnormal  hydrological  years  were  1986, 1988,  and  1993. 


ERDC/EL  TR-07-14 


66 


UMRS,  V.  americana 


1  100  200  300  400 


Time  (d) 


Figure  15.  Comparison  of  historical  and  simulated  data  on  biomass  of  plants  and  tubers  of  V.  americana  in 
the  Upper  Mississippi  River  System,  Wl.  (A)  Historical  data:  measured  plant  data  Pool  9, 1980, 
average  water  depth  1.1  m  (lat  43°  10'  N,  long  91°  30'  W;  Donnermeyer  1982).  Simulations: 
Nominal  initial  biomass  data;  light  extinction  coefficients  derived  from  10-year  average  background 
total  suspended  solids  values  measured  in  the  nearby  Pool  4  in  the  1980s  (Bartell  et  al.  1999); 
climatological  data,  Minnesota/St.Paul,  Wl  (long  93°  W,  lat  45°  N).  (B)  Water  level,  annual  average 
historical  data  set;  climate,  average  1985-94;  (C)  Water  level,  daily  values  0.5  m  depth  class; 
climate  1985;  (D)  Water  level,  daily  values  0.5  m  depth  class;  climate,  1986;  (E)  Water  level:  daily 
values  0.5-m  depth  class;  climate  1993;  (F)  Water  level,  daily  values  0.5-m  depth  class;  climate, 
1986. 


1980,  and  it  was,  therefore,  concluded  that  calculation  of  the  annual  aver¬ 
age  water  depth  from  Pool  9  historical  data  might  have  yielded  an  errone¬ 
ous,  probably  too  high,  average  depth  and,  therefore,  the  simulated  plant 
biomass  was  too  low.  The  next  simulation  at  a  constant,  shallower  0.5-m 
water  depth  indicated  that  in  this  case  more  biomass  was  produced  and 
that  two  tuber  classes  were  finished  (Figure  15B).  Surprisingly,  far  more 
biomass  could  be  produced  and  five  tuber  classes  were  finished  under  a 
normal  water  level  fluctuation  regime  in  Pool  8  pointing  to  a  tentative 
positive  influence  of  relatively  small  water  level  fluctuations  (Figure  15C). 
Normal  flooding  inhibited  biomass  and  tuber  production  somewhat, 
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allowing  only  four  tuber  classes  to  be  finished  but  the  population  to  persist 
(Figure  15D).  The  relatively  small  size  of  this  effect  was  attributed  to  the 
fortuitous  timing  of  the  high  water  levels  that  occurred  only  in  spring  and 
autumn,  still  allowing  the  plants  to  fully  benefit  from  the  high  summer  ir- 
radiance  at  normal  water  levels.  Abnormal  flooding,  however,  reduced  the 
finished  number  of  tuber  classes  by  a  factor  of  >2  in  the  0.5-m  depth  class 
(Figure  15  E),  and  completely  prevented  tuber  formation  in  the  l-m  depth 
class.  The  harshness  of  this  effect  was  attributed  to  the  fact  that  the  plants 
could  not  fully  benefit  from  the  high  summer  irradiance  due  to  the  high 
summer  water  levels.  The  effect  of  the  1988  drought  was  most  surprising 
and  detrimental.  In  this  year  substantial  plant  biomass  could  be  produced 
peaking  relatively  early  in  the  growth  season.  However,  tuber  formation 
was  severely  inhibited  because  (1)  water  levels  were  kept  relatively  high 
later  in  summer,  possibly  as  a  water  conservation  measure,  causing  in¬ 
creased  extinction  of  light  within  the  water  column,  and  (2)  temperatures 
were  relatively  high  causing  increased  respiration  and  senescence.  The  in¬ 
creased  light  extinction  in  the  water  column  may  even  have  been  larger  in 
situ  than  in  the  simulation  because  during  droughts  not  only  water  levels 
may  change  but  also  extinction  within  the  water  column  increases  by 
stimulation  of  algal  blooms.  The  seasonal  changes  in  the  light  extinction 
coefficient  were  kept  the  same  in  all  simulations. 

Application  2:  Importance  of  current  velocity  and  epiphyte  cover  for  a 
V.  americana  community  in  a  riverine  environment 

VALLA  was  used  to  explore  the  habitat  quality  at  a  site  in  the  UMR,  Turtle 
Island  (WI)  for  V.  americana.  This  site  is  characterized  by  shallow  water 
(depth  of  0.2  m  at  the  day  of  harvest),  current  velocities  tolerable  for  SAV 
(ranging  from  0.03  to  0.37  m  s_1),  and  high  turbidity  in  the  plant  growth 
season  (light  extinction  coefficients  ranging  from  3.00  to  4.23  nr1).  Re¬ 
sults  of  a  run  performed  for  this  site  with  a  temperate  climate  in  2001  in¬ 
dicated  that  simulated  shoot  biomass  at  the  day  of  harvest,  without  and 
with  corrections  for  epiphyte  shading  effects,  was  within  the  range  of  the 
measured  shoot  biomass,  but  with  corrections  for  current  velocity  and 
both  epiphyte  shading  and  current  velocity  it  was  below  the  measured 
range  (Figure  16A).  All  simulated  peak  shoot  biomass  values  were  within 
the  range  of  the  measured  values.  Peak  shoot  biomass,  without  corrections 
for  the  effects  of  current  velocity  and  epiphyte  shading,  was  similar  to 
mean  shoot  biomass  measured  at  the  day  of  harvest,  but  was  about  30 
days  delayed  in  time.  Peak  biomass  generated  by  runs  corrected  for  cur¬ 
rent  velocity  effects  was  50  percent  lower,  corrected  for  epiphyte  shading 
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10  percent  lower,  and  corrected  for  both  effects  64  percent  lower  than 
peak  biomass  without  corrections.  Thus,  simulated  plant  biomass  was 
greatly  decreased  by  onsite  current  velocity  and  epiphyte  cover.  Simulated 
end-of-year  tuber  number  was  always  substantial,  ranging  from  80  to 
around  350  tubers  nr2.  This  indicates  that  the  V.  americana  population 
would  persist  because  the  end-of-year  tuber  number  was  >1,  and  the  tuber 
size  was  large  enough  to  enable  sprouts  to  become  self-supporting  in  their 
carbon  gain  at  shallow  sites  (0.5  m  depth  at  relatively  turbid  water  com¬ 
mon  in  the  UMR;  Best  et  al.  2005). 


Figure  16.  Simulated  biomass  of  shoots  (A)  and  tubers  (B)  of  V.  americana  in  Upper  Mississippi  River  Pool  8 
at  Turtle  Island  (lat  43°  10'  N,  long  91°  30'  W).  Plant  biomass  0;  tuber  size  0.09  g  DW;  tuber  bank 
density  233  nr2;  site-specific  field  data  on  light  extinction  coefficient  (3.00-4.23  nr1),  water 
temperature  (5  °C  -30.5  °C),  current  velocity  (0.03-0.37  m  s1),  water  depth  (0.20  m),  2001; 
climatological  data  LaCrosse,  Wl,  2001  (long  and  lat  as  above);  measured  shoot  biomass  values 
are  indicated  by  the  mean  (solid  circle)  and  range  (between  arrows;  Bestetal.  2005). 


POTAM 

POTAM  calibration:  Simulated  and  measured  behavior  of  a  P.  pectinatus 
community  in  the  Western  Canal  near  Zandvoort,  NL 

The  seasonal  changes  in  biomass  of  plants  and  tubers,  as  simulated  by 
POTAM,  are  shown  in  Figure  17.  Simulated  plant  biomass  compared  well 
with  average  plant  biomass  measured  in  the  Western  Canal  near  Zand¬ 
voort,  NL  (Best  and  Boyd  2003c).  Plant  biomass  reached  its  maximum  13 
days  later,  and  peak  biomass  was  somewhat  higher  in  the  simulated  than 
in  the  measured  SAV,  notably  101.9  versus  78.5  g  DW  nr2.  The  higher 
simulated  than  measured  biomass  was  attributed  to  the  use  of  air  tem¬ 
peratures  for  this  nominal  run  instead  of  the  measured  water  tempera¬ 
tures.  Air  temperatures  with  a  lag-period  of  7  days  (default)  were  used  be¬ 
cause  the  temperature  of  the  water  surrounding  the  majority  of  the 
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Figure  17.  Simulated  biomass  of  plants  (A),  dormant  and  new  tuber  numbers  (B),  and  measured  plant 
biomass  (C)  of  P.  pectinatus  in  the  Western  Canal  near  Zandvoort,  NL.  Nominal  run.  Field 
data  1987  from  Best  and  Boyd  (2003a);  climatological  data  1987,  De  Bilt,  NL  (lat  52°  06'  N, 
long  05°  11'  E);  water  depth  1.3  m;  light  extinction  coefficient  1.07  m1. 


plant  shoots  in  summer  was  closer  to  the  temperature  of  the  air  than  to  the 
temperature  of  the  water  filling  the  canals  via  upward  seepage.  Another 
explanation  may  be  that  actual  maximum  plant  biomass  occurred  at  the 
end  of  August,  as  found  for  the  2.5-m  depth  class  (Best  and  Boyd  2003c), 
but  since  biomass  was  measured  only  once  a  month  the  actual  peak  was 
missed. 


Running  POTAM  for  the  same  canal  and  year,  but  with  both  plants  and 
tubers  initially  present,  showed  that  peak  plant  biomass  was  greatly  in¬ 
creased  but  no  tuber  class  was  finished.  The  large  increase  in  peak  biomass 
is  due  to  the  ability  of  the  plant  community  to  fully  capture  the  high  spring 
irradiance  at  this  latitude  (520  N),  which  they  can  not  without  wintering 
shoots.  However,  since  maintenance  respiration  increased  proportionally 
to  plant  biomass  and  light  extinction  due  to  self-shading  was  high,  carbo¬ 
hydrate  reserves  invested  in  tubers  were  insufficient  to  finalize  a  large¬ 
sized  tuber  class  and,  consequently,  such  a  plant  population  would  become 
extinct  the  next  year.  Thus,  wintering  shoots  would  not  provide  a  distinct 
advantage  for  this  plant  species  in  a  temperate  climate.  A  simulation 
started  with  initial  plant  biomass  but  from  a  smaller  tuber  size/concurrently 
initiated  tuber  number  combination  (a  0.07  g  DW  and  6  tubers  per  plant) 
yielded  a  far  higher  peak  plant  biomass  and  50  percent  reduced  tuber 
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numbers  (data  not  shown)  compared  to  the  values  measured  in  the  West¬ 
ern  Canal  (Best  and  Boyd  2003c). 

POTAM  validation:  Simulated  and  measured  behavior  of  a  P.  pectinatus 
community  in  Lake  Veluwe,  NL,  in  two  consecutive  years  with  greatly 
different  turbidities 

Lake  Veluwe,  located  in  the  center  of  NL,  is  a  man-made,  shallow  wind- 
exposed,  eutrophic  lake  in  which  remnants  of  P.  pectinatus  populations 
remained  at  the  end  of  the  1770s  (Bick  and  Van  Schaik  1980).  The  decline 
in  SAV  coverage  area  was  attributed  to  increased  turbidity  due  to  eutro¬ 
phication,  and  management  strategies  were  sought  to  reverse  this  decline. 
Management,  including  regular  flushing  of  the  lake  and  a  reduction  in  ex¬ 
ternal  phosphorus  loading,  resulted  in  a  more  transparent  water  column 
and  a  slow  increase  in  the  area  colonized  by  P.  pectinatus.  Because  this 
management  scenario  did  not  lead  to  a  full  recovery  of  the  SAV,  several 
studies  have  been  undertaken  to  quantify  the  relationship  between  light 
availability  within  the  water  column  and  the  production  of  plant  biomass 
and  tubers  in  this  lake.  Field  data  on  plant  and  tuber  biomass,  and  on  local 
water  transparency  generated  by  the  studies  mentioned  above  have  been 
used  to  validate  the  data  simulated  using  POTAM  (Van  Dijk  et  al.  1992; 
Van  Dijk  and  Achterberg  1992). 

POTAM  was  run  for  two  consecutive  years  greatly  differing  in  turbidity, 
i.e.,  1986  and  1987,  using  as  initial  values  the  measured  tuber  bank  densi¬ 
ties  and  predominant  tuber  sizes.  Model  results  indicated  that  a  higher 
plant  biomass  but  lower  tuber  class  number  (four  classes  in  1986  versus 
seven  classes  in  1987)  would  be  expected  in  1986  than  in  1987,  with  plant 
biomass  and  newly  produced  tuber  densities  similar  to  measured  ones 
(Figure  18).  The  timing  at  which  the  simulated  maxima  in  plant  biomass 
and  tuber  density  occurred  coincided  with  the  measured  values  in  1986, 
but  were  delayed  in  1987  (Figure  18).  This  example  illustrates  that  al¬ 
though  intuitively  it  would  be  expected  that  a  larger  peak  plant  biomass 
would  lead  to  a  higher  tuber  density,  it  is  possible  that  in  reality  less  tubers 
are  formed  since  the  latter  process  does  not  only  depend  on  plant  biomass 
but  also  on  the  critical  combination  of  day  length  and  temperature. 
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Figure  18.  Simulated  and  measured  biomass  of  P.  pectinatus  plants  and  tubers  in  Lake  Veluwe,  NL 
(lat  52°  17'  N,  long  05°  19'  E),  during  two  successive  years  differing  in  water  trans¬ 
parency.  Field  data  from  Van  Dijk  et  al.  (1992);  1986,  tuber  size  0.034  g  DW,  tuber  bank 
density  79  nr2,  water  depth  0.5  m;  1987,  tuber  size  0.024  g  DW,  tuber  bank  density 
100  nrr2,  water  depth  0.5  m.  Light  extinction  coefficients:  1986,  average  May-September 
2.58  nrr1;  1987,  average  4.30  m1  (Van  Dijk  and  Achterberg  1992).  Climatological  data 
1986,  1987,  De  Bilt,  NL  (lat  52°  06'  N,  long  05°  11'  E). 


POTAM  example  application:  Simulated  and  measured  behavior  of  a 
P.  pectinatus  community  at  other  latitudes 

To  investigate  whether  POTAM  could  be  used  to  simulate  the  behavior  of  a 
P.  pectinatus  community  at  other  sites  besides  the  nominal  one,  runs  were 
made  for  one  more  western  site,  the  Byrne  Canal,  CA,  and  one  more 
southern  site,  tropical  Lake  Ramgarh,  India. 

A  simulation  was  performed  of  a  P.  pectinatus  community  in  the  Byrne 
Canal,  CA,  starting  from  site-specific  community,  water  depth,  transpar¬ 
ency,  and  climatological  data.  For  this  site,  plant  community-specific  data 
included  the  following:  initial  plant  biomass  absent,  tuber  size  0.025  g 
DW,  concurrently  initiated  tuber  number  3  plant1,  and  tuber  bank  density 
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700  nr2.  Environmental  data  included  the  following:  water  depth  0.2  m 
and  light  extinction  coefficient  0.4  nr1.  In  these  conditions  simulated  plant 
biomass  showed  a  maximum  of  130  g  DW  nr2  using  the  nominal  self¬ 
shading  coefficient  (X)  of  0.095  m2  g  DW1,  but  a  maximum  close  to  the 
measured  biomass  of  175  g  DW  nr2  at  a  lower  self-shading  coefficient  of 
0.0183  m2  g  DW1  (as  reported  by  Sher-Kaul  et  al.  1995,  for  clear  Lake  Ge¬ 
neva,  Switzerland).  In  both  simulations  many  tuber  classes  were  finished, 
giving  rise  to  extremely  high  tuber  densities  -  just  as  found  in  situ  (Figure 
19B;  D.  F.  Spencer  1990,  unpub.  results).  Another  simulation  was  per¬ 
formed  of  a  P.  pectinatus  community  in  Lake  Ramgarh,  India,  starting 
from  a  nominal  community,  water  depth,  transparency,  but  climatological 
data  pertaining  to  Patancheru,  India.  For  this  site,  plant  community- 
specific  data  included  the  following:  no  initial  plant  biomass,  tuber  size 
0.083  g  DW,  concurrently  initiated  tuber  number  8  plant1,  and  a  tuber 
bank  density  of  240  nr2.  In  these  conditions  simulated  plant  biomass  was 
high,  maximally  375  g  DW  nr2,  and  only  two  tuber  classes  were  finished 
within  a  year.  Tuber  weights  and  numbers  of  the  Lake  Ramgarh  SAV  were 
not  published,  so  comparison  between  simulated  and  measured  tuber  data 
is  not  possible.  However,  simulated  maximum  plant  biomass  in  this  case 
was  in  the  same  order  of  magnitude  as  the  measured  peak  biomass  of 
445  g  DW  nr2  (Figure  19C).  Comparison  of  biomass  production  character¬ 
istics  in  the  various  climatological  conditions  (Figure  19)  indicates  that  in 
a  temperate  climate  generally  less  biomass  is  produced,  but  investment  in 
vegetative  reproduction  is  higher  than  in  a  tropical  climate.  The  criteria  for 
tuber  formation,  of  (a)  timed  just  after  flowering,  (b)  day  length  between 
12  and  16  hours,  and  (c)  temperature  between  5.5  °C  and  25  °C,  were  met 
for  long  continuous  periods  in  NL  and  CA,  but  only  occasionally  in  India 
because  the  temperature  usually  exceeded  the  critical  range.  In  NL,  tubers 
could  be  initiated  from  day  215  to  250  at  relatively  moderate  irradiance, 
allowing  only  one  tuber  class  to  be  finished;  in  CA  from  day  161  to  250  at 
far  higher  irradiance  giving  rise  to  many  tuber  classes;  in  India  from  day 
150  to  220  only  two  tuber  classes  might  be  finished.  The  extremely  limited 
window  for  tuber  formation  in  tropical  regions  was  confirmed  by  Pilon 
(1999),  who  reported  that  the  number  of  tubers  formed  per  plant  de¬ 
creases  from  around  o  at  a  latitude  of  25°N  to  8  at  68°  N,  but  predominant 
tuber  size  decreases  from  0.022  to  0.006  g  DW  tuber1.  Although  the  latter 
characteristics  were  suggested  at  first  to  have  a  genetic  basis,  it  was  re¬ 
ported  later  that  they  may  be  environmentally  regulated  (Pilon  and 
Santamaria  2002). 
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Figure  19.  Simulated  biomass  of  plants  and  tubers  of  P.  pectinatus  at  sites  differing  in  latitude. 
(A)  The  Western  Canal  near  Zandvoort,  NL  (lat  52°  06'  N,  long  05°  11'  E);  tuber  size 
0.083  g  DW,  tuber  bank  density  240  rrr2;  water  depth  1.3  m;  light  extinction 
coefficient  1.07  nr1;  climatological  data  1987;  measured  data  1987,  Best  and  Boyd 
(2003a).  (B)  Byrne  Canal,  CA  (lat  38°  32'  N,  long  121°  47'  W);  tuber  size  0.025  g  DW, 
tuber  bank  density  700  nr2;  K  value  0.0183  m2  g  DW1;  water  depth  0.2  m;  light 
extinction  coefficient  0.4  nr1;  climatological  data  1990;  measured  data  1990,  D.  W. 
Spencer,  USDA-ARS,  Davis,  CA,  unpublished).  (C)  Lake  Ramgarh,  India  (lat  26°  05'  N, 
long  83°  26'  E);  tuber  size  0.083  g  DW,  tuber  bank  density  240  nr2;  K  value 
0.0183  m2  g  DW1;  water  depth  1.3  m;  light  extinction  coefficient  1.07  nr1; 
climatological  data  Patancheru,  India,  1978  (lat  17°  27'  N,  long  78°  28'  E);  measured 
data  1968,  Sahai  and  Sinha  (1973). 


HYDRIL 

HYDRIL  calibration;  Simulated  and  measured  behavior  of  the  H.  verticillata 
community  in  Lake  Orange,  FL 

The  seasonal  changes  in  biomass  of  plants  and  tubers,  as  simulated  by 
HYDRIL,  are  shown  in  Figure  20.  Simulated  plant  biomass  compared  well 
with  measured  plant  biomass  found  in  the  lake  (Bowes  et  al.  1979).  Plant 
biomass  reached  its  maximum  somewhat  earlier  in  the  simulated  than  in 
the  measured  SAV.  The  lower  simulated  than  measured  biomass  was  at¬ 
tributed  to  the  rather  low  frequency  of  field  observations  (no  measure¬ 
ments  between  September  and  November).  Simulated  tuber  number  was 
well  within  the  tuber  number  range  of  approximately  500  nr2  found  in  the 
lake. 
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Figure  20.  Simulated  biomass  of  plants  (A),  dormant  and  new  tuber  numbers  (B),  and  measured 
plant  biomass  (C)  of  H.  verticillata  in  Lake  Orange,  FL.  Nominal  run.  Field  data  1977 
from  Bowes  et  al.  (1979);  climatological  data  1980,  Gainesville,  FL  (lat  29°  38'  N, 
long  82°  22'  W);  water  depth  1.0  m;  light  extinction  coefficient  0.83  m 1 . 


HYDRIL  validation:  Simulated  and  measured  behavior  of  the  H.  verticillata 
community  in  Lake  Trafford,  FL 

To  investigate  whether  HYDRIL  could  be  used  to  simulate  the  behavior  of 
a  H.  verticillata  community  at  other  sites,  a  run  was  made  using  initial 
values  of  another  Florida  lake  at  a  similar  latitude,  Lake  Trafford.  The 
simulated  biomass  compared  well  with  the  measured  biomass  (Figure  21). 
Comparison  of  the  tuber  behavior  is  not  useful  in  this  case  since  measured 
tuber  data  had  such  a  wide  scatter  that  no  conclusions  could  be  drawn. 
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Figure  21.  Simulated  biomass  of  plants  and  tubers  of  H.  verticillata  at  sites  differing  in 

latitude.  (A)  Lake  Orange,  FL  (lat  29°  38'  N,  long  82°  22'  W);  tuber  size  0.1  g  DW, 
tuber  bank  density  500  nr2;  water  depth  1.0  m;  light  extinction  coefficient 
0.83  m1;  climatological  data  Gainesville  1980  (long  and  lat  as  above);  measured 
data  1977,  Bowes  et  al.  (1979).  (B)  Lake  Trafford,  FL  (lat  26°  30'  N,  long  82° 

30'  W);  initial  plant  biomass  200  g  DW  nr2;  climatological  data  Tampa  (lat  27° 

58'  N,  long  82°  32'  W),  average  1980-85;  measured  data  1977,  Bowes  et  al. 
(1979).  (C)  Lake  Ramgarh,  India  (lat  26°  05'  N,  long  83°  26'  E;  initial  plant 
biomass  154  g  DW  nr2;  climatological  data  Patancheru,  India  (lat  17°  27'  N, 
long  78°  28'  E),  1978;  measured  data  1970s,  Sahai  and  Sinha  (1973). 

HYDRIL  example  applications 

Application  1:  Simulated  and  measured  behavior  of  the  H.  verticillata 
community  at  other  latitudes 

Simulation  of  a  H.  verticillata  community  in  Lake  Ramgarh,  India,  indi¬ 
cated  that  usually  lower  peak  biomass  than  in  FL  is  reached,  when  started 
merely  from  subterranean  tubers.  This  was  attributed  to  the  formation  of 
more  tubers  and  their  inherent  sink  function  for  assimilates  in  India.  A 
simulation  started  from  wintering  plant  biomass  in  addition  to  tubers  gen¬ 
erated  a  far  higher  peak  biomass  (Figure  21C).  The  latter  simulation 
showed  biomass  curves  that  agreed  well  with  the  onsite  measured  biomass 
(Figure  21C).  Using  the  model  calibrated  on  North  Florida  to  calculate  the 
timing  of  phenological  events  in  other  climatological  conditions  indicated 
postponement  of  anthesis  and  reduced  tuber  formation  in  autumn  only  in 
temperate  areas,  but  earlier  anthesis  and  extensive  tuber  formation  in  two 
periods  of  the  year  in  the  tropics.  Both  of  these  model  results  are  con¬ 
firmed  by  literature  data  from  India  and  Ireland,  respectively  (data  not 
shown;  Sahai  and  Sinha  1973;  Scannel  and  Webb  1976). 
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Application  2:  Simulated  effects  of  cutting  and  grazing 

Effects  of  man-made  control  activities,  such  as  plant  biomass  harvesting  at 
different  times  and  at  various  water  depths  to  prevent  excessive  biomass 
production  can  also  be  calculated  using  the  HYDRIL  model  (Table  5). 
Thus,  in  the  latter  case  the  model  can  be  used  as  a  tool  for  aquatic  plant 
management  agencies.  From  the  data  in  this  table  it  can  be  concluded  that 
harvesting  at  the  end  of  July  to  a  water  depth  of  0.8  m  requires  removal  of 
a  relatively  low  amount  of  biomass  but  yields  the  lowest  tuber  bank  den¬ 
sity  at  the  end  of  the  year.  This  situation  can  be  seen  as  favorable  to  con¬ 
trol  H.  verticillata.  In  contrast,  harvesting  later  in  the  year  requires  re¬ 
moval  of  relatively  more  plant  biomass  and  allows  for  a  relatively  higher 
tuber  bank  density.  Removing  only  the  top  layer  of  the  plant  community 
later  in  the  year  may  lead  even  to  increased  numbers  of  tubers  at  the  end 
of  the  year,  probably  due  to  a  higher  light  penetration  within  the  commu¬ 
nity  in  the  period  of  tuber  formation.  Effects  of  plant  biomass  control 
through  grazing  of  tubers  by  waterfowl  can  be  calculated  also  by  reducing 
the  initial  tuber  bank  density  in  the  simulations. 


Table  5.  Simulated  effects  of  mechanical  harvesting  date  and  depth  on  H.  verticillata  plant  biomass  and  tuber 
bank  density.  Results  of  a  1-year  simulation  starting  from  a  tuber  bank  density  of  500  tubers  nr2  without 
wintering  shoot  biomass  present,  in  Lake  Orange,  FL;  climatological  data  1980,  Gainesville,  FL. 


Harvest 

Day 

Depth 

(m) 

Live  Biomass 
Day  260 
(g  DW  m-2) 

Preharvest 
Biomass 
(g  DW  nr2) 

Postharvest 
Biomass 
(g  DW  nr2) 

Day  with  Zero 
Biomass 

Final  Tuber 
Biomass 
(g  DW  nr2;  (no)) 

212 

0.8 

0.0 

72.8 

7.65 

257 

51  (330) 

243 

0.8 

6.3 

127.5 

13.2 

273 

56  (490) 

273 

0.8 

160.3 

124.0 

12.3 

287 

65  (490) 

273 

0.1 

160.3 

124.0 

95.5 

>365 

82  (804) 

304 

0.8 

160.3 

65.2 

6.5 

316 

73  (647) 

MILFO 

MILFO  calibration:  Simulated  and  measured  behavior  of  a  M.  spicatum 
community  in  Lake  Wingra,  Wl 

The  seasonal  changes  in  biomass  of  plant  shoots  and  roots,  and  of  the  rhi¬ 
zome/root  crown  system,  as  simulated  by  MILFO,  are  shown  in  Figure  22. 
Simulated  shoot  biomass  compared  well  with  shoot  biomass  measured  in 
Lake  Wingra.  Peak  biomass  appeared  to  be  reached  somewhat  earlier  in 
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the  simulation  than  found  in  the  lake.  The  latter  was  attributed  to  the  low 
frequency  of  field  observations  (no  measurements  between  September  and 
November).  The  simulated  biomass  of  the  rhizome/root  crown  system 
showed  two  maxima  per  year.  Variation  was  within  the  range  found  in  a 
M.  spicatum  community  in  the  same  lake  in  later  years  (Smith  and  Adams, 
1986). 
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Figure  22.  Simulated  biomass  of  shoots,  roots  (A),  rhizomes/root  crowns  (B),  and  measured 
shoot  biomass  (C)  of  M.  spicatum  in  Lake  Wingra,  Wl.  Nominal  run.  Field  data  on 
shoot  biomass,  1970,  from  Adams  and  McCracken  (1974),  and  on  initial 
rhizomes/root  crown  biomass,  1970,  from  Smith  and  Adams  (1986); 
climatological  data  1970,  Lake  Wingra,  Wl  (lat  43°  08'  N,  long  89°  20'  W);  water 
depth  1.5  m;  light  extinction  coefficient  1.15-2.00  (Lee  and  Kluesener  1972). 


Running  MILFO  with  the  low  assimilate  requirements  suggested  as  typical 
for  SAV  by  Spencer  et  al.  (1997),  showed  that  peak  biomass  of  M.  spicatum 
shoots  would  increase  by  approximately  20  percent,  oscillations  in  bio¬ 
mass  of  the  rhizome/root  crown  system  would  be  greater,  and  final  bio¬ 
mass  of  the  rhizome/root  crown  system  would  be  somewhat  increased 
(Figure  23).  However,  as  indicated  under  “Maintenance,  growth,  and  as¬ 
similate  partitioning,”  a  construction  cost  of  0.99  to  1.11  for  M.  spicatum 
plant  tissues  is  on  the  low  side,  taking  the  usually  high  N  concentrations  of 
shoots  into  consideration. 
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Figure  23.  Simulated  biomass  with  assimilate  requirement  for  growth  (ASRQ)  value  nominal 
(A)  or  as  suggested  for  SAV  (B),  and  measured  shoot  biomass  (C)  of  a  M.  spicatum 
community  in  Lake  Wingra,  Wl  (SAV  ASRQ  data;  Spencer  et  al.  1997).  Initial 
biomass  and  climatological  data  as  in  nominal  run. 

MILFO  validation:  Simulated  and  measured  behavior  of  a  M.  spicatum 
community  at  a  more  southern  latitude 

To  investigate  whether  MILFO  could  be  used  to  simulate  behavior  of  a 
M.  spicatum  community  at  other  sites,  runs  were  made  for  a  site  at  a  more 
southern  latitude,  Lake  Guntersville,  AL.  Behavior  of  M.  spicatum  in  this 
lake  is  particularly  interesting  because  the  lake  is  long,  oriented  from 
northeast  to  southwest,  and  situated  at  a  latitude  of  340  N,  being  very  close 
to  tropical  (33°N).  The  biomass  dynamics  of  M.  spicatum  communities  in 
this  lake  have  been  described  as  highly  variable  in  time  and  space  (Grace 
and  Wetzel  1978;  Stanley  et  al.  1976);  unfortunately,  in  these  descriptions 
no  attention  was  paid  to  differences  in  latitude  of  the  various  sites  within 
the  lake  nor  to  local  differences  in  temperature  or  other  environmental 
factors.  It  was  mentioned,  however,  that  flowering  and  subsequent  slough¬ 
ing  periods  were  less  predictable  in  southern  locales  than  in  northern 
ones. 

Initial  plant  biomass  values  measured  at  a  site  in  this  lake  studied  in  1990 
were  very  low,  possibly  due  to  grass  carp  herbivory  the  previous  year 
(M.  S.  Stewart,  personal  communication,  Waterways  Experiment  Station, 
Vicksburg),  and  initial  rhizome/root  crown  mass  for  the  simulation  has, 
therefore,  been  set  at  the  critical  value  of  12  g  DW  nr2.  Rooting  depth  in 
the  simulation  was  kept  at  1.5  m,  although  actual  water  depth  may  have 
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varied  over  1.0  +  0.7  m  within  a  year  (Stanley  et  al.  1976).  Simulated  bio¬ 
mass  of  the  first  plant  cohort  remained  low  (Figure  24).  Only  one  apparent 
biomass  peak  could  be  distinguished,  which  originated  from  the  second 
plant  cohort.  Simulated  shoot  biomass  coincided  in  timing  with  measured 
shoot  biomass,  but  the  simulated  peak  was  lower  than  the  measured  one. 
The  latter  difference  was  attributed  to  temporal  decreases  in  water  depth 
in  the  lake,  while  water  depth  was  kept  constant  in  the  simulation,  leading 
to  an  underestimate  of  simulated  plant  biomass.  Relatively  small  changes 
in  water  depth  can  cause  large  changes  in  net  assimilation  and  biomass 
production  (See  “Sensitivity  Analysis  of  the  Models”  and  “Discussion”).  To 
investigate  which  consequences  a  warm  year  for  the  M.  spicatum  commu¬ 
nity  in  this  lake  might  have,  when  three  instead  of  two  plant  cohorts  are 
expected  to  be  active,  a  model  run  was  made  with  the  same  initial  biomass 
and  climatological  data  but  with  a  third  plant  cohort  activated.  It  turned 
out  that  such  a  SAV  was  able  to  generate  far  higher  shoot  biomass  values 
of  approximately  950  g  DW  nr2  in  one  year,  similar  in  timing  and  value  to 
maximum  biomass  values  reported  for  the  nearby  Melton  Hill  Lake 
(Stanley  et  al.  1976),  with  rhizome/root  crown  biomass  accumulating 
towards  the  end  of  the  year.  However,  similar  biomass  values  could  also 
be  reached  earlier  in  the  year  when  higher  initial  (nominal)  biomass 
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Time  (d)  Time  (d)  Time  (d) 

Figure  24.  Simulated  biomass  of  shoots,  roots  (A),  rhizomes/root  crowns  (B),  and  measured  shoot 
biomass  (C)  of  M.  spicatum  in  Lake  Guntersville,  AL.  Field  data  on  shoot  and  initial 
rhizomes/root  crown  biomass  from  M.  S.  Stewart  1990,  unpublished)  extrapolated  to 
Julian  day  number  1.  Climatological  data  1990,  Lake  Guntersville,  AL  (lat  34°  08'  N, 
long  85°  50'  W). 
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values  were  used  as  input,  and  only  two  cohorts  active;  in  the  latter  case, 
biomass  peaks  of  both  plant  cohorts  could  be  distinguished  and  rhizome/ 
root  crown  biomass  was  well  above  the  critical  level  but  not  accumulating. 

MILFO  example  application:  Simulated  and  measured  behavior  of  a  M. 
spicatum  community  at  other  latitudes 

Running  MILFO  with  nominal  biomass  values  and  climatological  data 
typical  for  sites  representative  for  temperate,  temperate  to  tropical,  and 
tropical  climates  (Figure  25)  indicated  the  following:  (1)  in  all  climates  one 
clear  biomass  peak  is  generated;  (2)  only  in  a  temperate  climate  the  bio¬ 
mass  peak  of  both  first  and  second  cohorts  can  be  distinguished  (flowering 
coinciding  with  every  biomass  maximum  is  always  a  suitable  indicator,  but 
it  is  often  not  noted  in  biomass  studies);  (3)  peak  biomass  is  expected  to  be 
highest  in  the  tropics  (that  highest  biomass  values  have  been  found  at 
northern  latitudes  may  be  because  most  biomass  studies  on  SAV  have 
been  performed  at  the  latter  latitudes  and  biomass  data  from  tropical  ar¬ 
eas  are  extremely  scarce);  and  (4)  end-of-year  accumulation  of  rhi¬ 
zome/root  crown  biomass  usually  occurs  in  tropical,  but  not  in  temperate 
climates  (that  is,  when  three  plant  cohorts  are  active). 
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Figure  25.  Simulated  biomass  of  shoots,  roots,  and  rhizomes  of  M.  spicatum  at  sites  differing  in 
latitude,  starting  from  a  nominal  initial  biomass.  (A)  Lake  Wingra,  Wl  (lat  43°  08'  N, 
long  89°  20'  W);  climatological  data  1970.  (B)  Lake  Guntersville,  AL  (lat  34°  08'  N, 
long  85°  50' W);  climatological  data  1990.  (C)  Patancheru,  India,  1978  (lat  17°  27'  N, 
long  78°  28'  E);  climatological  data  1978. 
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The  maximum  shoot  biomass  generated  under  tropical  conditions  was 
somewhat  higher  than  the  measured  range  of  288  to  640  g  DW  nr2  and 
the  rhizome/root  crown  biomass  varied  over  a  range  close  to  the  measured 
range  of  32  to  160  (Zutschi  and  Vass  1973).  The  higher  simulated  shoot 
biomass  in  this  case  might  partly  be  attributed  to  the  use  of  climatological 
data  from  Patancheru,  which  is  located  more  south,  and  thus  warmer,  than 
the  Kashmir  lakes  from  which  the  measured  biomass  values  originated 
(Patancheru  170  N,  Kashmir  Lake  320  N);  however,  more  northern  clima¬ 
tological  data  from  India  were  not  available.  It  was  also  investigated 
whether  M.  spicatum  may  benefit  from  adaptation  to  the  tropics  by  pro¬ 
ducing  thinner  leaves.  This  was  done  because  a  higher  leaf  surface  area  : 
dry  weight  ratio  (K  value)  has  been  found  for  milfoil  in  Japan  (0.01  m2  g 
DW1;  Ikusima  1970)  than  in  WI  (0.006  m2  g  DW1;  Titus  and  Adams 
1979a).  Timing  was  very  similar  and  simulated  plant  biomass  about  ten 
percent  higher  using  the  higher  K  value  (data  not  shown). 
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5  Sensitivity  Analysis  of  the  Models 

A  sensitivity  analysis  of  a  simulation  model  is  required  to  assess  the  pa¬ 
rameters  likely  to  strongly  affect  model  behavior.  The  current  analysis  is 
based  on  the  effect  of  a  change  in  one  parameter  when  all  other  parame¬ 
ters  are  kept  the  same.  The  parameter  under  study  was  changed  and  l-year 
simulations  were  conducted  under  nominal  environmental  conditions.  As 
reference  level,  the  nominal  parameter  values  were  chosen  as  presented  in 
Tables  2a,  2b,  2c,  and  2d.  The  results  were  compared  with  those  of  a 
nominal  run.  Each  parameter  was  once  increased  by  20  percent  and  once 
decreased  by  20  percent.  The  relative  sensitivity  (RS)  of  a  parameter  was 
defined  as  the  relative  change  in  the  variable  on  which  the  effect  was 
tested  divided  by  the  relative  change  in  the  parameter  (Ng  and  Loomis 
1984).  The  effects  of  nine  parameters  on  two  state  variables,  representing 
different  plant  biomass  compartments,  were  tested.  A  model  variable  is 
considered  sensitive  to  a  change  in  the  value  of  a  parameter  at  RS>0.5  and 
<-0.5. 


RS  = 


(yield t  -  yields  )  /yield,. 
(parami  - paramr) / paramr 


(7) 


where 

yieldi  =  value  at  parameter  value  i; 
yieldr  =  value  at  reference  parameter  value  r; 
parami  and  paramr  as  above 

Among  the  four  models,  VALLA  was  the  most  and  HYDRIL  the  least  sensi¬ 
tive  to  parameter  changes  (Table  6).  Among  the  models  for  tuber-formers, 
VALLA  was  far  more  sensitive  than  POTAM  and  HYDRIL.  In  the  models, 
maximum  plant  biomass  was  generally  more  sensitive  to  changes  in  the 
tested  parameter  values  than  end-of-year  tuber  number  in  tuber-forming 
plants  or  rhizome/root  crown  biomass  in  M.  spicatum  (Table  6).  Maxi¬ 
mum  plant  biomass  was  sensitive  to  changes  in  all  parameters  in  VALLA, 
but  only  to  changes  in  DVRRT,  NPL,  INTUB,  AMX  and  EE  in  HYDRIL. 
End-of-year  tuber  number/rhizome  biomass  was  sensitive  to  changes  in 
almost  all  parameters,  except  increased  NPL  and  ROC,  and  decreased 
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INTUB,  while  it  was  only  sensitive  to  changed  NPL  and  decreased  RTR  in 
HYDRIL. 

Table  6.  Relative  sensitivity  (RS)  of  two  state  variables  to  deviations  in  parameter  values  from  their  nominal 
values  as  presented  in  Table  2.  The  RS  of  a  parameter  is  the  relative  change  in  the  variable  on  which  the  effect 
was  tested  divided  by  the  relative  change  in  the  parameter.  A  model  variable  is  considered  sensitive  to  a 
change  in  the  value  of  a  parameter  at  RS>0.5  and  <-0.5.  Results  were  obtained  in  1-year  simulations  under 
nominal  conditions.  Sensitive  values  shaded. 


Parameter 

Name 

Parameter 

Value 

Relative  Sensitivity 

VALLA 

POTAM 

HYDRIL 

MILFO 

Max. 

Plant 

Biomass 

EOY2 

Tuber 

No. 

Max. 

Plant 

Biomass 

EOY2 

Tuber 

No 

Max. 

Plant 

Biomass 

EOY2 

Tuber 

No. 

Max. 

Plant 

Biomass 

EOY2 

Rhiz. 

Biomass 

DVRVT 

+20% 

0.56 

-2.5 

-1.36 

-3.36 

-0.20 

0 

-0.23 

-0.81 

-20% 

-6.04 

-1.39 

-0.91 

4.91 

-0.44 

0 

-0.26 

-0.99 

DVRRT 

+20% 

0.98 

-2.47 

-0.39 

-0.43 

-1.04 

0 

-0.69 

-0.89 

-20% 

-2.19 

0.24 

-0.45 

-3.12 

-2.19 

0 

-0.79 

0.66 

NPL 

+20% 

3.39 

-0.01 

0.28 

1.20 

0.62 

-0.62 

0.16 

-0.79 

-20% 

-0.82 

2.71 

0.35 

1.14 

0.73 

-0.54 

-0.16 

0.79 

INTUB 1 

+20% 

3.25 

-1.79 

0.25 

0 

0.50 

0 

0.20 

0.17 

-20% 

-0.92 

-0.03 

0.34 

0.19 

0.62 

0 

0.22 

0.18 

ROC 

+20% 

2.65 

-0.43 

0.09 

0 

0.02 

0 

0.19 

0.17 

-20% 

-1.37 

2.33 

0.14 

0 

0.03 

0 

0.21 

0.18 

AMX 

+20% 

5.00 

4.46 

1.72 

-1.58 

3.40 

0 

1.96 

2.00 

-20% 

3.02 

2.04 

1.94 

5.00 

3.16 

0 

1.97 

2.02 

EE 

+20% 

0.50 

-0.73 

0.25 

-0.83 

1.47 

0 

1.10 

1.14 

-20% 

0.56 

1.44 

0.33 

-3.10 

1.61 

0 

1.22 

1.25 

RTR 1 

+20% 

1.76 

-0.77 

-0.10 

-2.15 

-0.11 

0 

-0.13 

0.57 

-20% 

-2.62 

2.19 

-0.10 

5.00 

-0.10 

0.97 

-0.14 

0.72 

RDR 

+20% 

2.25 

0.71 

0 

0 

0 

0 

-0.62 

-1.01 

-20% 

-3.03 

0.22 

0 

-2.93 

0 

0 

-0.77 

-1.36 

Notes:  1  Synonymous  terms  in  MILFO:  IWGRIZ  and  TRAFAC;  2  EOY,  end-of-year. 


The  sensitivity  of  maximum  plant  biomass  and  end-of-year  tuber  number/ 
rhizome  biomass  to  changes  in  environmental  factors  was  assessed  follow¬ 
ing  the  same  approach  as  for  sensitivity  analysis  of  the  model  parameters. 
For  this  purpose,  parameter  changes  were  based  on  value  ranges  taken 
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from  literature,  which  sometimes  differed  more  than  20  percent  from  the 
nominal  parameter  values  presented  in  Table  2.  Also  in  this  analysis, 

VALLA  was  the  most  sensitive  among  the  models  (Table  7).  In  the  models, 
maximum  plant  biomass  was  generally  more  sensitive  to  changes  in  the 
tested  parameter  values  than  end-of-year  tuber  number  in  tuber-forming 
plants  or  rhizome/root  crown  biomass  in  M.  spicatum.  Changes  in  light 
reflection  coefficient  at  the  water  surface  ( RC )  had  no  significant  effect, 
whereas  the  effects  of  changes  in  the  other  tested  parameter  values  de¬ 
creased  in  the  order  of  LT>  DPTT. 

Table  7.  Environmental  factor  analysis,  expressed  as  relative  sensitivity  (RS)  of  two  state  variables  to 
deviations  in  parameter  values  from  their  nominal  values  as  presented  in  Table  2.  The  RS  of  a  parameter  is 
the  relative  change  in  the  variable  on  which  the  effect  was  tested  divided  by  the  relative  change  in  the 
parameter.  A  model  variable  is  considered  sensitive  to  a  change  in  the  value  of  a  parameter  at  RS>0.5  and 
<-0.5.  Results  were  obtained  in  1-year  simulations  under  nominal  conditions.  Sensitive  values  shaded. 


Parameter 

Name 

Parameter 

Value 

Relative  Sensitivity 

VALLA 

POTAM 

HYDRIL 

MILFO 

Max. 

Plant 

Biomass 

EOY2 

Tuber 

No. 

Max. 

Plant 

Biomass 

EOY2 

Tuber 

No. 

Max. 

Plant 

Biomass 

EOY2 

Tuber 

No. 

Max. 

Plant 

Biomass 

EOY2  Rhiz. 
Biomass 

Climate1 

Lat  42°  N 

- 

- 

Lat  27°  N 

-0.49 

-0.87 

Climate1 

Lat  52°  N 

- 

- 

Lat  38°  N 

-1.54 

1.43 

Climate1 

Lat  29°  N 

- 

- 

Lat  17°  N 

1.66 

0.47 

Climate1 

Lat  43°  N 

- 

- 

Lat  34°  N 

-0.01 

0.89 

RC3 

1.00 

-0.06 

-0.06 

-0.06 

-0.06 

-0.06 

-0.03 

-0.06 

-0.04 

0.00 

-0.43 

-0.05 

-0.02 

0.09 

-0.09 

0 

-0.07 

-0.07 

LT 

+20% 

2.09 

0.04 

-0.12 

0.18 

-0.44 

0 

0.92 

-1.01 

-20% 

-2.79 

0.66 

-0.08 

0.43 

0.44 

0 

-1.12 

-1.01 

DPTT 

+20% 

1.47 

-2.16 

-0.02 

0 

-0.07 

0 

-0.31 

-0.31 

-20% 

-2.43 

0.48 

0.03 

0 

-0.54 

0 

-0.34 

-0.33 

Notes: 1  Climates  at  Chenango  Lake,  NY  (lat  42°  N),  Ft.  Lauderdale,  FL  (lat  27°  N),  De  Bilt,  NL  (lat  52°  N),  Davis,  CA  (lat  38°  N), 
Gainesville,  FL  (lat  29°  N),  Patancheru,  India  (lat  17°  N),  Lake  Wingra,  Wl  (lat  43°  N)  and  Lake  Guntersville,  AL  (lat  34°  N). 

2  EOY,  end-of-year;  3  Light  reflection  coefficient  at  water  surface;  to  enable  calculation  of  the  relative  sensitivity,  a  very  low  value 
of  0.000001  was  used. 
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6  Discussion 

The  models  can  be  used  to  simulate  the  daily  changes  in  carbon-flow- 
based  plant  processes  and  biomass  under  various  site-specific  conditions 
and  management  scenarios  over  a  l-  to  5-year  period  of  four  SAV  species, 
in  climates  varying  from  temperate  to  tropical.  VALLA,  POTAM,  and 
HYDRIL  can  be  run  with  default  or  user-specified  input  values  on  plant 
biomass,  tuber  size/concurrently  initiated  tuber  number,  and  tuber  bank 
density,  and  MILFO  with  user-specified  input  values  on  plant  and 
rhizome/root  crown  biomass.  For  all  runs,  climate-  and  site-specific 
environmental  variable  data  can  be  entered  by  the  user,  and  scenarios 
interesting  to  the  user  can  be  explored. 

Developmental  cycle 

In  the  models,  phenology  is  tied  to  day-degree  sum.  This  enables  simula¬ 
tions  for  different  sites  and  climates  facilitating  model  operation  by  users 
who  do  not  possess  a  full  data  set  on  plant  characteristics  and  environ¬ 
mental  variables  for  the  water  body  for  which  they  desire  to  run  the  models. 
Earlier  or  later  flowering  biotypes  are  suited  to  different  environments. 

The  effect  of  flowering  date  can  be  tested  with  the  model  by  varying  the 
development  rate  of  the  vegetation.  Slower  rates  represent  later  and  faster 
rates  earlier  biotypes.  Development  rate  slower  or  faster  than  the  nominal 
rate  leads  to  lower  biomass.  Faster  development  leads  to  a  shorter  growing 
season  and  less  biomass,  incomplete  light  interception,  and  lower  carbo¬ 
hydrate  availability  for  organ  formation.  At  the  same  time,  however,  the 
rate  of  organ  formation  increases  but  the  duration  of  each  organ  formation 
shortens.  Intuitive  prediction  of  biotype  behavior  under  such  highly  vari¬ 
able  climatic  conditions  is,  therefore,  hazardous.  The  models  show  some 
promise  in  being  able  to  reproduce  some  of  these  complex  responses  of  the 
vegetation  and  may  be  useful  in  evaluating  long  term  implications  of  dif¬ 
ferences  in  development  rate.  Because  only  vegetative  reproduction  is  in¬ 
cluded  in  the  life  cycle  formulation,  neither  the  tentative  difference  in  im¬ 
portance  of  sexual  reproduction  between  climates  nor  the  tentative 
importance  of  seeds  long-term  survival  (>5  years)  in  highly  variable  envi¬ 
ronmental  conditions  such  as  drought  and  extreme  flooding  can  be  ex¬ 
plored  with  the  current  model  versions. 
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Carbon-flow-based  aspects 

The  models  are  based  on  the  carbon  flow  through  a  particular  plant  com¬ 
munity,  without  limitation  by  nitrogen  and  phosphorus.  Calibration  has 
been  done  on  field  data  pertaining  to  lakes  with  high  alkalinity  and  rela¬ 
tively  high  pH,  leaving  light  as  major  limiting  factor  followed  by  DIC.  The 
results  of  the  sensitivity  analysis  showed  that  the  models  are  very  sensitive 
to  changes  in  process  parameters  influencing  carbon  flow,  with  VALLA  be¬ 
ing  most  sensitive  and  HYDRIL  least.  The  higher  sensitivity  of  VALLA 
compared  to  the  other  models  may  be  a  consequence  of  its  “meadow-type” 
vertical  biomass  distribution  over  the  water  column  that  would  make  it 
more  susceptible  to  plant-process  and  environment-related  changes  in 
light  absorption  in  the  overlying  water  layers  than  the  “canopy-type”  dis¬ 
tribution  in  the  other  models. 

Unique  characteristics  of  the  models 

The  models  contain  the  following  unique  descriptions  of 

1.  Species-characteristic  vertical  distribution  of  shoot  biomass  in  the  water 
column  that  enables  the  calculation  of  the  fraction  of  irradiance  actually 
available  for  absorption  by  the  SAV. 

2.  Recalculation  procedures  of  this  vertical  distribution  with  daily  changes  in 
water  level. 

3.  Recalculation  procedures  of  this  vertical  distribution  with  daily  removal  of 
shoot  biomass  at  various  levels  within  the  water  column. 

4 .  Species-characteristic  epiphytic  light  interception  that  accommodates  ac¬ 
curate  estimates  of  the  fraction  of  irradiance  available  for  absorption  by 
the  SAV  at  sites  differing  in  epiphyte  cover. 

5.  Species-characteristic  effects  of  current  velocity  on  photosynthesis  that  ac¬ 
commodates  accurate  estimates  of  the  fraction  of  AMY  accomplished  by 
the  SAV  at  sites  differing  in  current  velocity. 

6.  Removal  of  periodic  shoot  and  tuber/root  crown  biomass  that  enables  es¬ 
timates  of  effects  of  mechanical  control  and  propagule  removal  by  grazing 
and  dredging  on  plant  biomass. 

7.  Relationships  of  plant  process  parameters  with  site-specific  climate  by 
linkage  with  formatted  weather  files  and  calculation  of  latitude-specific  ef¬ 
fects  on  day  length  that  enables  estimates  of  effects  of  different  climates. 
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Application  potential 

The  models  can  be  used  to  simulate  the  daily  changes  in  plant  processes 
and  biomass  under  various  site-specific  conditions  and  management  sce¬ 
narios  over  a  l-  to  5-year  period  of  four  SAV  species,  in  climates  varying 
from  temperate  to  tropical. 

Example  applications  have  illustrated  how  relatively  low  and  high  fre¬ 
quency  fluctuations  in  water  levels,  current  velocity,  and  epiphyte  cover 
might  affect  submersed  plant  populations,  without  taking  plant- 
characteristic  adaptation  into  consideration.  Ability  of  plants  to  adapt  to 
these  changes  may  be  important  characteristics  for  their  persistence  in 
rivers,  reservoirs,  and  estuaries.  Although  this  ability  is  a  rather  inten¬ 
sively  discussed  research  topic,  pertinent  ecological  data  are  currently 
largely  lacking. 

The  models  can  be  run  using  air  or  water  temperatures  as  inputs.  The  re¬ 
lationship  between  air  and  water  temperature  may  differ  from  the  model 
calibration,  because  temperatures  within  each  water  body  are  influenced 
by  catchment  morphometry,  wind  speed,  fetch,  mixing  processes,  and  up¬ 
ward  seepage.  In  the  authors’  experience,  a  lag  period  of  7  days  is  usually  a 
reasonable  estimate  for  shallow  lakes. 

The  models  have  been  developed  as  stand-alone  versions.  They  can  rela¬ 
tively  easily  be  modified  to  communicate  with  other  models,  because  they 
are  programmed  in  FORTRAN77  and  their  structure  is  simple.  For  in¬ 
stance,  VALLA  version  1.0  was  linked  via  a  Geographical  Information  Sys¬ 
tem  through  an  appropriate  interface  to  hydrodynamic  and  sediment 
transport  models,  and  this  model  combination  was  used  to  explore  habitat 
suitability  in  Peoria  Lake,  IL,  for  V.  americana  (Black  et  al.  2003;  Best  et 
al.  2004).  VALLA  and  POTAM  versions  1.0  were  linked  to  a  decision  sup¬ 
port  system  that  was  used  for  ecological  risk  assessment  on  the  UMR 
(Bartell  et  al.  2000,  2003).  Currently  tools  are  being  developed  that  would 
greatly  facilitate  model  application  by  providing  links  to  (1)  Websites 
where  standardized  weather  files  can  be  obtained,  and  (2)  websites  where 
gauging  and  water  quality  data  required  for  the  input  file  can  be  selected/ 
obtained  (U.S.  only). 

Predecessors  of  the  models  have  been  described  in  technical  reports  (Ver¬ 
sions  1.0;  Best  and  Boyd  1996, 1999a,  1999b,  2001a,  2001b,  2003a, 

2003b;  Boyd  and  Best  1996).  A  summary  of  the  versions  3.0  was  published 
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recently  (Best  and  Boyd  2007).  Model  versions,  user  manuals,  and  execu¬ 
table  model  versions  are  available  via  the  internet  at 
http://el.erdc.usace.mil/products.cfm?Topic=model&Tvpe=aauatic. 
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7  Conclusions 

A  dynamic  simulation  modeling  approach  to  describing  carbon-flow-based 
ecophysiological  processes  and  biomass  dynamics  of  freshwater  sub¬ 
mersed  aquatic  plant  species  has  been  developed. 

The  models  (VALLA,  POTAM,  HYDRIL,  and  MILFO)  describe  major,  car¬ 
bon  flow-based  ecophysiological  processes  and  biomass  dynamics  of  four 
common  freshwater  species  and  how  these  are  influenced  by  factors  such 
as  light,  temperature,  current  velocity,  DIC,  availability,  oxygen  concentra¬ 
tion,  and  human  influences  such  as  management  measures  (changes  in 
turbidity,  mechanical  harvesting,  grazing,  flooding). 

Model  validation  results  showed  a  good  fit  with  measured  biomass  tem¬ 
perate  and  tropical  climates.  Validation  for  tropical  climate  was  limited 
greatly  by  the  scarcity  in  available  data  pertaining  to  tropical  areas.  Sensi¬ 
tivity  analysis  showed  that  the  models  are  very  sensitive  to  changes  in 
process  parameters  influencing  carbon  flow,  with  VALLA  being  most  sen¬ 
sitive  and  HYDRIL  least. 
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